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Abstract: 

We review the techniques necessary for the calculation of virtual electroweak and soft 
photonic corrections at the one-loop level. In particular we describe renormalization, 
calculation of one-loop integrals and evaluation of one-loop Feynman amplitudes. We 
summarize many explicit results of general relevance. We give the Feynman rules and the 
explicit form of the counter terms of the electroweak standard model, we list analytical 
expressions for scalar one-loop integrals and reduction of tensor integrals, we present the 
decomposition of the invariant matrix element for processes with two external fermions 
and we give the analytic form of soft photonic corrections. These techniques are applied 
to physical processes with external W^-bosons. We present the full set of analytical for- 
mulae and the corresponding numerical results for the decay width of the ly-boson and 
the top quark. We discuss the cross section for the production of W^-bosons in e'^e~- 
annihilaton including all 0{a) radiative corrections and finite width effects. Improved 
Born approximations for these processes are given. 
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1 Introduction 



All known experimental facts about the electroweak interaction are in agreement with 
the Glashow-Salam- Weinberg (GSW) model [1, 2, 3, 4]. Therefore, this theory is called 
the standard model (SM) of electroweak physics. Despite its extraordinary experimental 
success it is by no means tested in its full scope. Many more experimental and theoretical 
efforts arc needed for its further confirmation. 

An important step in this direction is provided by the e^e~ colliders SLC and LEPlOO 
which started a new era of precision experiments. The first important results from these 
experiments were the determination of the number of light neutrinos and the precise 
measurement of the mass of the neutral weak gauge boson, the Z-boson [5]. Furthermore 
the total and partial widths of the Z-boson and various on-resonance asymmetries have 
been determined and will be measured with increasing accuracy. These experiments will 
uniquely allow to study in great detail all the properties of the Z-boson and its couplings 
to fermions. 

There are, however, ingredients of the electroweak SM, which are not directly accessible 
at SLC and LEPlOO. The most important one is probably the gauge boson self-interaction 
which is crucial for the nonabelian structure of the GSW model. It will be directly tested 
for the first time at LEP200, the upgraded version of LEP. There the center of mass 
energy will be high enough to produce pairs of charged weak gauge bosons, the W^-bosons, 
such that one can study the reaction e^e' — > W^W~ in great detail. It will allow the 
investigation of the nonabelian three-gauge boson interactions 'yW~^W~ and ZW^W~ at 
the classical level of the theory. Moreover, all the properties of the VT-boson, like its mass 
and its total and partial widths can be measured directly there. The statistics will not be 
as good as on the Z-peak. One expects of the order of 10^ l^-pairs and thus an accuracy 
at the percent level. The examination of several independent methods indicates that an 
error of about 0.1% for the VT-mass determination can be reached [6]. 

Theoretical predictions should have an accuracy comparable to or even better than 
the experimental errors. If the experimental precision is of the order of one percent the 
classical level of the theory is no longer sufficient. One is forced to take into account quan- 
tum corrections: the radiative corrections. In the case of the electroweak SM these can 
reach several percent. For the high precision experiments at LEPlOO even the first order 
corrections are inadequate, one has to take into account leading higher order corrections, 
too. 

Radiative corrections arc not only compelling for the precise comparison between 
the theoretical predictions and the experimental results, but offer the possibility to get 
informations about sectors of the theory that are not directly observable. While the direct 
investigation of certain objects may not be possible because the available energy is too 
small to produce them they may affect the radiative corrections noticeably. 

In the electroweak SM there arc at least two such objects. The top quark, the still 
undiscovered constituent of the third fermion generation, and the Higgs boson, the phys- 
ical remnant of the Higgs-Kibble mechanism of spontaneous symmetry breaking. Both 
particles seem to be too massive to be produced directly in the existing colliders. How- 
ever, the high precision experiments performed so far together with the precise knowledge 
of the radiative corrections of the electroweak SM already allow to derive limits on the 
mass of the top quark within the SM [7, 5] . Since the sensitivity of radiative corrections 
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to the mass of the Higgs boson is weaker, the restrictions on this parameter are at present 
only marginal [8]. The situation may improve with increasing experimental accuracy. 
While direct determinations of physical parameters are in general to a large extent model 
independent, the information extracted from radiative corrections depends on the entire 
structure of the underlying theory. 

Finally there is a third important issue concerning radiative corrections. It is likely 
that the electroweak SM, despite its experimental success, is only an effective theory, 
the low-energy approximation of a more general structure. This would manifest itself 
typically in small deviations from the SM predictions. Furthermore most of the presently 
discussed new physics is connected with scales bigger than the experimentally accessible 
energies. Therefore new phenomena will show up predominantly via indirect effects rather 
than via direct production of new particles. In order to disentangle these small effects 
one has to know once again the predictions of the SM accurately and thus needs radiative 
corrections. 

The actual evaluation of the radiative corrections is a tedious and time consuming task. 
It requires extensive calculations involving many different techniques, like renormalization, 
evaluation of loop integrals, Dirac algebra calculations, phase space integrations and so 
on. Fortunately the whole procedure can be organized into different independent steps. 
Furthermore many steps can be facilitated with the help of computer algebra [9, 10, 11, 12] . 

For the interesting processes at LEPlOO radiative corrections have been calculated by 
many authors [13]. Their structure is relatively simple since the masses of the external 
fermions can be neglected. Calculations for gauge boson production processes at LEP200 
are already more complicated because the masses of the external gauge bosons are non- 
negligible. Such calculations have been performed by several groups and we will give 
the most important results in the second part of this review. The whole complexity of 
one-loop corrections will show up when considering reactions where all external particles 
are massive like e.g. gauge boson scattering processes which may be investigated at the 
LHC or SSC. The calculation of radiative corrections to these processes has just started. 

In the first part of this review we collect the relevant formulae and techniques necessary 
for the calculation of electroweak one-loop radiative corrections. Although we discuss 
everything in the context of the SM the presented material is - apart from the explicit 
form of the renormalization constants - applicable to extended models as well. In the 
second part these methods are applied to physical processes with external l^-bosons. This 
part not only gives examples for the calculation of one-loop electroweak corrections, but 
also provides a survey on the status of radiative corrections for the production and the 
decay of VF-pairs in e'^e~ annihilation. The corresponding experiments will be carried 
through in a few years at LEP200. 

The general techniques described in this paper are restricted to the virtual part of the 
electroweak corrections and soft photon bremsstrahlung. We do not consider the methods 
appropriate for hard photon bremsstrahlung. This can be efliciently treated using spinor 
techniques [14] and Monte Carlo simulations [15]. Furthermore we do not touch the 
methods developed for calculating higher order QCD corrections. 

This paper is organized as follows: 

In chapter 2 we specify the Lagrangian of the electroweak SM. Chapter 3 outlines 
the on-shell renormalization for the physical sector of the electroweak SM and provides 
explicit expressions for the counter terms. All relevant formulae for the calculation of one- 
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loop Feynman integrals are collected in chapter 4. In chapter 5 we introduce the standard 
matrix elements, a concept which allows to represent the results for one-loop diagrams in 
a systematic and simple way. In chapter 6 we show how everything is put together in the 
actual calculation of one- loop amplitudes and provide first simple examples. The relevant 
formulae for the calculation of the soft photon corrections are summarized in chapter 7. 
Chapter 8 serves to define our input parameters and the way of resumming higher order 
corrections. 

The remaining chapters are devoted to applications. In chapter 9 we give results for 
the width of the H^-boson, in chapter 10 for the width of the top quark. Finally the 
radiative corrections to the production of VF-pairs in e^e~ annihilation are discussed in 
chapter 11. 

The appendices contain the Feynman rules of the electroweak SM, the explicit expres- 
sions for the self energies of the physical particles and the vertex functions as well as the 
bremsstrahlung integrals relevant for the W^-boson and top quark decay width. 
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2 The Glashow-Salam- Weinberg Model 



The Glashow-Salam- Weinberg (GSW) model of the electroweak interaction has been 
proposed by Glashow [1], Weinberg [2], and Salam [3] for leptons and extended to the 
hadronic degrees of freedom by Glashow, Iliopoulos and Maiani [4] . It is the presently most 
comprehensive formulation of a theory of the unified electroweak interaction: theoretically 
consistent and in agreement with all experimentally known phenomena of electroweak ori- 
gin. For energies that are small compared to the electroweak scale it reproduces quantum 
electrodynamics and the Fermi model, which already accomplished a good description of 
the electromagnetic and weak interactions at low energies. It is minimal in the sense that 
it contains the smallest number of degrees of freedom necessary to describe the known 
experimental facts. 

The electroweak standard model (SM) is a nonabelian gauge theory based on the non- 
simple group SU{2)w X U{1)y- From experiment we know that three out of the four 
associated gauge bosons have to be massive. This is implemented via the Higgs-Kibble 
mechanism [16]. By introducing a scalar field with nonvanishing vacuum expectation 
value the SU{2)]y x f/(l)y gauge symmetry is spontaneously broken in such a way that 
invariance under the electromagnetic subgroup U{l)em is preserved. The SM is chiral 
since right- and left-handed fcrmions transform according to different representations of 
the gauge group. Consequently fermion masses are forbidden in the symmetric theory. 
They are generated through spontaneous symmetry breaking from the Yukawa couplings. 
Diagonalization of the fermion mass matrices introduces the quark mixing matrix in the 
quark sector. This can give rise to CP-violation. Fermions appear in generations. The 
model does not fix their number, but from experiment we know that there are exactly 
three with light neutrinos [5] . 

The SM is a consistent quantum field theory. It is renormalizable, as was proven by 
't Hooft [17], and free of anomalies. Therefore it allows to calculate unique quantum 
corrections. Given a finite set of input parameters measurable quantities can be predicted 
order by order in perturbation theory. 

The classical Lagrangian Cc of the SM is composed of a Yang-Mills, a Higgs and a 
fermion part 

Cc — Cym + Ch + C,F ■ (2.1) 
Each of them is separately gauge invariant. They are specified as follows: 

2.1 The Yang-Mills-part 

The gauge fields are four vector fields transforming according to the adjoint represen- 
tation of the gauge group SU{2)w x U{1)y. The isotriplet W^, a = 1,2,3 is associated 
with the generators I^, of the weak isospin group SU{2)w, the isosinglet with the 
weak hypercharge Yw of the group U{1)y- The pure gauge field Lagrangian reads 

Cym = -\ {d,W: - d.Wl H- g^e^^'WlWlf - \ {d,B, - d.B,f , (2.2) 
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where e^^'^ are the totally antisymmetric structure constants of SU{2). Since the gauge 
group is non-simple there are two gauge coupling constants, the SU{2)w gauge coupling 
g2 and the U{1)y gauge coupling gi. The covariant derivative is given by 

D, = d^- ig2l^W; + ig.^B,. (2.3) 

The electric charge operator Q is composed of the weak isospin generator 7^ and the 
weak hypercharge according to the Gell-Mann Nishijima relation 

Q = /^ + ^- (2.4) 



2.2 The Higgs part 

The minimal Higgs sector consists of a single complex scalar SU{2)w doublet field 
with hypercharge IV — 1 

*'^) -(:;':,')■ 

It is coupled to the gauge fields with the covariant derivative (2.3) and has a self coupling 
resulting in the Lagrangian 

£^ = (L'^$)^(OT)-T/($). (2.6) 

The Higgs potential 

V($) = ^ ($t$^^ _^2^t$ (2.7) 

is constructed in such a way that it gives rise to spontaneous symmetry breaking. This 
means that the parameters A and are chosen such that the potential V($) takes its 
minimum for a nonvanishing Higgs field, i.e. the vacuum expectation value ($) of the 
Higgs field is nonzero. 



2.3 Fermionic Part 



The left-handed fermions of each lepton (L) and quark [Q) generation are grouped 
into SU{2)w doublets (we suppress the colour index) 

L'- - = [^^ Q']^ uM, = {^. (2.8) 

the right-handed fermions into singlets 

I'f = uJ+l'j, u'f = uj+u'j, d'f = uj+d'j, (2.9) 

where u;± — is the projector on right- and left-handed fields, respectively, j is the 
generation index and u, I, u and d stand for neutrinos, charged leptons, up-type quarks 
and down-type quarks, respectively. The weak hypercharge of the right- and left-handed 
multiplets is chosen such that the known electromagnetic charges of the fermions are 
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reproduced by the Gell-Mann-Nishijima relation (2.4). There are no right-handed neutri- 
nos. These could be easily added, but they would induce nonvanishing neutrino masses, 
which have not been observed experimentally so far. 
The fermionic part of the Lagrangian reads 

i 

+ E (jf^rD.l'f + ^fil^D,v!f + dlfi^^D.d'f) (2.10) 

i 

- E (p'G'rj^f^ + Q'G'l^f^ + Q-Gyf^ + h.c) . 

ij 

Note that in the covariant derivative acting on right-handed fermions the term in- 
volving g2 is absent, since they are SU{2)w singlets. The primed fermion fields are by 
definition eigenstates of the electroweak gauge interaction, i.e. the covariant derivatives 
are diagonal in this basis with respect to the generation indices. Glj, G^- and Gfj are 

the Yukawa coupling matrices, $ = (0°*, — 0~)^ is the charge conjugated Higgs field and 
(f)~ = {4>~^)* ■ The SU{2)^r X U{1)y symmetry forbids explicit mass terms for the fermions. 
The masses of the fermions are generated through the Yukawa couplings via spontaneous 
symmetry breaking. 



2.4 Physical fields and parameters 

The theory is constructed such that the classical ground state of the scalar field satisfies 

I(^)r = ^ = y7^0. (2.11) 

In perturbation theory one has to expand around the ground state. Its phase is chosen 
such that the electromagnetic gauge invariance U{l)em is preserved and the Higgs field is 
written as 

where the components 0"*", H and x have zero vacuum expectation values, (f)'^, (f)~ and x 
are unphysical degrees of freedom and can be eliminated by a suitable gauge transforma- 
tion. The gauge in which they are absent is called unitary. The field H is the physical 
Higgs field with mass 

Mh = V2fx. (2.13) 

Inserting (2.12) into Cc the vacuum expectation value v introduces couplings with mass 
dimension and mass terms for the gauge bosons and fermions. 
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The physical gauge boson and fermion fields are obtained by diagonalizing the corre- 
sponding mass matrices 




where 



/ y-'^w cw J J ' (2.14) 

Ji — '-'ik J ki 

fR — Tjf>RfiR 
Ji — '^ik J k ' 

Cw = COS 9w — —? , s^ = sin^vK, (2-15) 

with the weak mixing angle 9w and / stands for u, I, u or d. The resulting masses are 

X 1 

Mw = -g2v, Mz = -^jgl + gl 

(2.16) 

The neutrinos remain massless since the absence of the right-handed neutrinos forbids the 
Yukawa couplings which would generate their masses. With (2.16) we find for the weak 
mixing angle 

cw = (2.17) 

Identifying the couphng of the photon field to the electron with the electrical charge 
e = -\/47ra yields 

e = , (2.18) 

or 

91 = —. 92 = —. (2.19) 
Cw sw 

The diagonalization of the fermion mass matrices introduces a matrix into the quark- W- 
boson couplings, the unitary quark mixing matrix 

Vij = Ur^'^Utf. (2.20) 

There is no corresponding matrix in the lepton sector. Since there is no neutrino mass 
matrix, U"'^ is completely arbitrary and can be chosen such that it cancels f/''^ in the 
lepton- W-boson couplings. The same would also be true for the quark sector if all up- 
type or down-type quarks would be degenerate in masses. For degenerate masses one can 
choose [/■^ = U^^ arbitrary without destroying the diagonality of the corresponding mass 
matrix and thus eliminate Vij. 

The above relations (2.13, 2.16, 2.18, 2.20) allow to replace the original set of 
parameters 

g„ 52, A, /.^ G\ G", (2.21) 
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by the parameters 

e, Mw, Mz, Mh, mf,i, V^j (2.22) 

which have a direct physical meaning. Thus we can express the Lagrangian (2.1) in terms 
of physical parameters and fields. 

Inserting (2.12) into Lc generates a term linear in the Higgs field H which we denote 
by tH{x) with 

t = v{ii'-^v'). (2.23) 

The tadpole t vanishes at lowest order due to the choice of v. We use t instead of v in 
the following. Choosing v as the correct vacuum expectation value of the Higgs field $ is 
equivalent to the vanishing of t. 

2.5 Quantization 

Quantization of Cc and higher order calculations require the specification of a gauge. 
We choose a renormalizable 't Hooft gauge with the following linear gauge fixings 

= {e.r'^d^Z, - MziiDh, (2.24) 

FT = (e7)-^aM^, 

leading to the following gauge fixing Lagrangian 
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fix 



{F^f + {F^ f + 2F+F-\ . (2.25) 



Cfix involves the unphysical components of the gauge fields. In order to compensate 
their effects one introduces Faddeev Popov ghosts u°'{x), u"{x) {a — ±,7,Z) with the 
Lagrangian 

Cfp^u-{x)-^^u^{x). (2.26) 
is the variation of the gauge fixing operators under infinitesimal gauge trans- 



se^ix) 

formations characterized by 9^{x). 

The 't Hooft Feynman gauge = 1 is particularly simple. At lowest order the poles 
of the ghost fields, unphysical Higgs fields and longitudinal gauge fields coincide with 
the poles of the corresponding transverse gauge fields. Furthermore no gauge-field-Higgs 
mixing occurs. 

With Cfix and jCfp the complete renormahzable Lagrangian for the electroweak SM 
reads 

^Gsw = Cc + Cfix + Cfp ■ (2.27) 
The corresponding Feynman rules are given in App. A. 
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3 Renormalization 



The Lagrangian (2.1) of the minimal SU{2)]y x f/(l)y model involves a certain num- 
ber of free parameters (2.22) which have to be determined experimentally. These are 
chosen such that they have an intuitive physical meaning at tree level (physical masses, 
couplings), i.e. they are directly related to experimental quantities. This direct relation 
is destroyed through higher order corrections. Moreover the parameters of the original 
Lagrangian, the so-called bare parameters, differ from the corresponding physical quanti- 
ties by UV-divergent contributions. However, in renormalizable theories these divergencies 
cancel in relations between physical quantities, thus allowing meaningful predictions. The 
renormalizability of nonabelian gauge theories with spontaneous symmetry breaking and 
thus of the SM was proven by 't Hooft [17]. 

One possibility to evaluate predictions of a renormalizable model is the following: 

• Calculate physical quantities in terms of the bare parameters. 

• Use as many of the resulting relations as bare parameters exist to express these in 
terms of physical observables. 

• Insert the resulting expressions into the remaining relations. 

Thus one arrives at predictions for physical observables in terms of other physical quan- 
tities, which have to be determined from experiment. In these predictions all UV- 
divergencies cancel in any order of perturbation theory. The predictions obtained from 
different input parameters differ in finite orders of perturbation theory by higher order 
contributions. This treatment of renormalization has been pioneered by Passarino, Velt- 
man and Consoli [18] and is the basis of the so-called 'star' scheme of Kennedy and Lynn 
[19]. 

We use the counterterm approach. Here the UV-divergent bare parameters are ex- 
pressed by finite renormalized parameters and divergent renormalization constants (coun- 
tcrtcrms). In addition the bare fields may be replaced by renormalized fields. The coun- 
tcrtcrms are fixed through renormalization conditions. These can be chosen arbitrarily, 
but determine the relation between renormalized and physical parameters. Further eval- 
uation proceeds like described above. The results depend in finite orders of perturbation 
theory not only on the choice of the input parameters but also on the choice of the 
renormalized parameters. Clearly the physical results are unambiguous up to the orders 
which have been taken into account completely. The renormalization procedure can be 
summarized as follows: 

• Choose a set of independent parameters (e.g. (2.22) in the SM). 

• Separate the bare parameters (and fields) into renormalized parameters (fields) and 
renormalization constants (see Sect. 3.1). 

• Choose renormalization conditions to fix the counterterms (see Sect. 3.2). 

• Express physical quantities in terms of the renormalized parameters. 

• Choose input data in order to fix the values of the renormalized parameters. 

• Evaluate predictions for physical quantities as functions of the input data. 
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The first three items in this hst specify a renormahzation scheme. 

Putting the counterterms equal to zero, the renormahzed parameters equal the bare 
parameters and we recover the first approach. 

However, wc can choose the counterterms such that the finite renormahzed parameters 
are equal to physical parameters in all orders of perturbation theory. This is the so- 
called on-shell renormahzation scheme. In the SM one uses the masses of the physical 
particles Mw, Mz, MH,mf, the charge of the electron e and the quark mixing matrix 
Vij as renormahzed parameters. This scheme was proposed by Ross and Taylor [20] 
and is widely used in the electroweak theory. The advantage of the on-shell scheme 
is that all parameters have a clear physical meaning and can be measured directly in 
suitable experiments^. Furthermore the Thomson cross section from which e is obtained 
is exact to all orders of perturbation theory. However, not all of the particle masses are 
known experimentally with good accuracy. Therefore other schemes may sometimes be 
advantageous. 

Renormahzation of the parameters is sufficient to obtain finite S-matrix elements, 
but it leaves Green functions divergent. This is due to the fact that radiative correc- 
tions change the normalization of the fields by an infinite amount. In order to get finite 
propagators and vertex functions the fields have to be renormahzed, too. Furthermore 
radiative corrections provide nondiagonal corrections to the mass matrices so that the 
bare fields are no longer mass eigenstates. In order to rediagonalize the mass matrices 
one has to introduce matrix valued field renormahzation constants. These allow to define 
the renormahzed fields in such a way that they arc the correct physical mass eigenstates 
in all orders of perturbation theory. If one does not renormalize the fields in this way, 
one needs a nontrivial wave function renormahzation for the external particles. This is 
required in going from Green functions to S-matrix elements in order to obtain a properly 
normalized S-matrix. 

The results for physical S-matrix elements arc independent of the specific choice of field 
renormahzation. There exist many different treatments in the literature [21, 22, 23, 24, 25]. 
Calculations without field renormahzation were performed by [26]. 



3.1 Renormalization constants and counterterms 



In the following we specify the on-shell renormalization scheme for the electroweak SM 
quantitatively. As independent parameters we choose the physical parameters specified in 
(2.22). The renormahzed quantities and renormalization constants are defined as follows 
(we denote bare quantities by an index 0) 



eo 


= Zee 


= + ' 




= M^ 


+ SM^, 




= M| 


+ SMl 




= M| 


+ SM],, 




= mf,i 





-'^This is not the case for the quark masses, due to the presence of the strong interaction. In practice 
these are replaced by suitable experimental input parameters (see Sect. 8.1). 
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Ui and U2 are unitary matrices since Vij^o and V^j are both unitary. 

Radiative corrections affect the Higgs potential in such a way that its minimum is 
shifted. In order to correct for this shift one introduces a counterterm to the vacuum 
expectation value of the Higgs field, which is determined such that the renormalized v is 
given by the actual minimum of the effective Higgs potential. Since we have replaced v by 
t (2.23) wc must introduce a counterterm 6t. This is fixed such that it cancels all tadpole 
diagrams, i.e. that the effective potential contains no term linear in the Higgs field H. 

The counterterms defined above are sufficient to render all S-matrix elements finite. 
In order to have finite Green functions we must renormalize the fields, too. As explained 
above we need field renormalization matrices in order to be able to define renormalized 
fields which are mass eigenstates 



W, 



Z^'^W^ = (1 + ^SZw)W 



± 




( 71/2 71/2 

^ZZ ^ZA 

7V2 7I/2 

\^AZ ^AA 




7l/2,/,L fL 
^ij Jj 



1 + ^SZzz ISZzA \fz\ 
jSZaz 1 + \^Zaa I 

1 + \bZH)E, 



(3.2) 



fR _ 71/2,/,^ /■i? _ / r I IX^/.-Rwi? 
Ji,0 - "^ij Jj - K^ij + 2^^U )jj ■ 

We do not discuss the renormalization constants of the unphysical ghost and Higgs 
fields. They do not affect Green functions of physical particles and are not relevant for 
the calculation of physical one-loop amplitudes. Furthermore the renormalization of the 
unphysical sector decouples from the one of the physical sector. It is governed by the 
Slavnov- Taylor identities. A discussion of this subject can be found e.g. in [24, 25]. 

In writing Z — 1 + 5Z ior the multiplicative renormalization constants (matrices) 
we can split the bare Lagrangian jCq into the basic Lagrangian C and the counterterm 
Lagrangian SC 

Cq = C + sc. (3.3) 

C has the same form as Cq but depends on renormalized parameters and fields instead 
of unrenormalized ones. 6C yields the counterterms. The corresponding Feynman rules 
are listed in App. A. They give rise to counterterm diagrams which have to be added to 
the loop graphs. Since we are only interested in one-loop corrections, we neglect terms of 
order {SZy everywhere. 



3.2 Renormalization conditions 



The renormalization constants introduced in the previous section are fixed by imposing 
renormalization conditions. These decompose into two sets. The conditions which define 
the renormalized parameters and the ones which define the renormalized fields. While 
the choice of the first affects physical predictions to finite orders of perturbation theory. 
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the second are only relevant for Green functions and drop out when calculating S-matrix 
elements. Nevertheless their use is very convenient in the on-shell scheme. They not only 
allow to eliminate the explicit wave function renormalization of the external particles, 
but also simplify the explicit form of the renormalization conditions for the physical 
parameters considerably. 

In the on-shell scheme all renormalization conditions are formulated for on mass shell 
external fields. The field renormalization constants, the mass renormalization constant 
and the renormalization constant of the quark mixing matrix are fixed using the one- 
particle irreducible two-point functions. For the charge renormalization we need one 
three-point function. For this we choose the ee7-vertex function. In the following renor- 
malized quantities are denoted by the same symbols as the corresponding unrenormalized 
quantities, but with the superscript". 

As discussed above the first renormalization condition involves the tadpole T, the 
Higgs field one-point amputated renormalized Green function 



(3.4) 
and simply states 

f ^T + St^O. (3.5) 

As a consequence of this condition no tadpoles need to be considered in actual calculations. 

Next we need the renormalized one-particle irreducible two-point functions. These are 
defined as follows (we are using the 't Hooft-Feynman gauge) 





w,, r77?\ ^^ 
k 



I* „ „ // ///U „ „ „ " _ rW 



= -^9,.{k^ - M^) - . - ^] S^(A;^) - ^^tY{P) 




k 



fj(^) 



where a,b^A, Z, M\ = 0, (3.6) 



f WZ) -V"{k)^^{k--Ml) + ^t"{k% 




12 



p 




fi 



ff,(p) 



The corresponding propagators are obtained as the inverse of these two-point functions. 
Note that we have to invert matrices for the neutral gauge bosons and for the fermions. 

The renormahzed mass parameters of the physical particles are fixed by the require- 
ment that they are equal to the physical masses, i.e. to the real parts of the poles of the 
corresponding propagators which are equivalent to the zeros of the one-particle irreducible 
two-point functions. In case of mass matrices these conditions have to be fulfilled by the 
corresponding eigenvalues resulting in complicated expressions. These can be considerably 
simplified by requiring simultaneously the on-shell conditions for the field renormalization 
matrices. These state that the renormahzed one-particle irreducible two-point functions 
are diagonal if the external lines are on their mass shell. This determines the nondiagonal 
elements of the field renormalization matrices. The diagonal elements are fixed such that 
the renormahzed fields are properly normalized, i.e. that the residues of the renormal- 
ized propagators are equal to one. This choice of field renormalization implies that the 
renormalization conditions for the mass parameters (in all orders of perturbation theory) 
involve only the corresponding diagonal self energies. Thus we arrive at the following 
renormalization conditions for the two-point functions for on-shell external physical fields 



= 0, 



0, 



0, 



rt^{k)e''{k) 



fc2=0 



0, 



r'^:^{k)e%k) 



0, 



hm — —pr- 

k^^M^ k^ - 



Ref;X(^K(fc) = -i£^(A;), 



lim 



1 



1 



ReT',.'{k)e-^{k) = -ie,{k), lun^-RcT^^^{k)e^{k) = -ie,{k), 



Rer"{k) 



k'^=Mlj 



Rer{j{p)uj{p) 



0, 



f,3 



0, 



lim 



^ + rrifi ~ - f 



J}^2^ k^ - Ml 
Reu,{p')%{p') 



Re f^(A;) 



0, 



hm^ ^Z,(p')Reff,(p')^i^ 



m 



iui{p'). 
(3-1) 

e(A:), u{p) and u{p') are the polarization vectors and spinors of the external fields. Re 
takes the real part of the loop integrals appearing in the self energies but not of the quark 
mixing matrix elements appearing there. Since we restrict ourselves to the one-loop order 
we apply it only to those quantities which depend on the quark mixing matrix at one loop. 
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In higher orders Re must be replaced by Re everywhere. Re and Re are only relevant 

above thresholds and have no effect for the two-point functions of on-shell stable particles. 
If the quark mixing matrix is real Re can be replaced by Re. This holds in particular for 
a unit quark mixing matrix which is often used. 

From the above equations we obtain the conditions for the self energy functions. 



ReE5r(M^) = 0, 
ReE|^(M|) = 0, 
E^^(O) = 0, 



ReE^^(M|) = 0, 
E^^(O) = 0, 2 



(3.8) 



Re 



Re 



ReE^(M|) = 0, 



0, 
0, 



Re 



Re 



aE^^(A; 



dk^ 

9E^(A;2) 



fe2=0 



0, 



0, 



(3.9) 



"^/jR-e E{f{mjj) + m/jRe Efj 



m 



0, 



^/jR'eEfj^(m^^^) + mf,iReE(f{mjj) = 0, 



ReE^'^(m^,J + ReEr(m^, 



f,Li 



(3.10) 



d 



+ '^rnj— {Re±{fip') + ReE^'^(/) + 2ReS^'^(p^ 



= 0. 



Note that the (unphysical) longitudinal part of the gauge boson self energies drops out 
for on-shell external gauge bosons. 

Our choice for the renormalization condition of the quark mixing matrix Vij can be 
motivated as follows. To lowest order Vij is given by (see eq. 2.20) 



Vo,ij = UTUl 



ik ^ kj 1 



(3.11) 



where the matrices U^'^ transform the weak interaction eigenstates /q to the lowest order 
mass eigenstates /o 

Uij'^f^^ = /to- (3-12) 

In the on-shell renormalization scheme the higher order mass eigenstates are related to 
the bare mass eigenstates through the field renormalization constants of the fermions 



2-l/2,/,L 



I J 



(3.13) 



We define the renormalized quark mixing matrix in analogy to the unrenormalized one 
through the rotation from the weak interaction eigenstates to the renormalized mass 



^This condition is automatically fulfilled due to a Ward identity. 
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eigenstates. In the one-loop approximation the rotation contained in the fermion wave 
function renormahzaton 1 + ^SZ^ is simply given by the anti-Hermitean part 5Z^^ of 

SZi]^''-liSZ!-'-dZ!f). (3.14) 
Thus we are lead to define the renormalized quark mixing matrix as 

v., = {6,, + I6z^^''''')u;:;^u^^:^{6r., + isz'^f") 



(3.15; 



It has been shown that this condition correctly cancels all one-loop divergencies and that 
Vij = Vo^ij in the limit of degenerate up- or down- type quark masses [27]. 

Finally the electrical charge is defined as the full ee7-coupling for on-shell external 
particles in the Thomson limit. This means that all corrections to this vertex vanish 
on-shell and for zero momentum transfer^ 




^ ieu{p)-f^u{p). (3.16) 



The momenta p, p' flow in the direction of the fermion arrows. Due to our choice for 
the field renormalization the corrections in the external legs vanish and we obtain the 
condition 

u{p)'^l^^iP,P)'^{p) 2_ ^=i'^u{p)'^f,u{p), (3.17) 
for the (amputated) vertex function 



K\p.p') 




(3.18) 



e ,p 



3.3 Explicit form of renormalization constants 

The renormalized quantities defined in Sect. 3.2 consist of the unrenormalized ones 
and the countcrtcrms as specified by the Fcynman rules in App. A. The renormalization 
conditions allow to express the counterterms by the unrenormalized self energies at special 
external momenta. This is evident for all renormalization constants apart from the one 

^Due to the wave function renormalization of the external particles the self energy corrections in the 
external legs contribute only with a factor 1 /2 to the S-matrix elements. 
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for the electrical charge. In this case, however, we can use a Ward identity to eliminate 
the vertex function. 

From conditions (3.5, 3.8, 3.9) we obtain for the gauge boson and Higgs sector 



5t = -T, 

5M2, = ReE^(M2,), 

5Ml =ReE|^(M|), 



5Zw — —Re 
5Zzz = -Re 



5Zaz = -2Re 



Ml 



6Z. 



ZA 



SZaa = 



fe2=0 



SZ, 



H 



-Re 



aE^(A;2) 



dk^ 



dk^ 

5Mjj = ReE^(M^), 

In the fermion sector (3.10) yields 

5mf, = ^Rc (e^ (mj,) + Er(mJ,) + 2E{f^f/)) , 



(3.19) 



Re 



+K. + ^.)sf K.)|' ^^^^ 



-Re 



m i . — m 



+2m/,im/jEf/(m5^^.) I , i ^ j. 



(3.20) 



= -ReE^'^(mJ,J - mJ^.^Re [E^'^(p2) + T^ip') + 2Er(/); 
= -Re^'^'l^,) - ^.^I^e [e^'^(/) + Tffip") + 2Er (p2) 



2 2 ' 



p^=m 



The use of Re ensures reality of the renormalized Lagrangian. Furthermore it yields 

5Zl^ = 5Zij {ml ^ m]) , (3.21) 

and in particular 



(3.22) 
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In the lepton sector we have Vij — 5ij. Consequently all lepton self energies are 
diagonal and the off-diagonal lepton wave function renormalization constants are zero. 
The same holds for the quark sector if one replaces the quark mixing matrix by a unit 
matrix as is usually done in calculations of radiative corrections for high energy processes. 

The renormalization constant for the quark mixing matrix Vij can be directly read off 
from (3.15) 

1 



sv,. 



(3.23) 



Inserting the fermion field renormalization constants (3.20) yields 



tj 2 



2 



Vu 



kj 



ik 



(3.24) 



^l,k^lk{^l,k) + ml,i^lk{^l,d 
Hmlk + mli) {^l%rnlk) + «j) 

^h^kti^h) + ^Ik^kti^lk) 

+ mlj) (Sjf (m|,) + S^f (mi J)' 

It remains to fix the charge renormalization constant 6Ze. This is determined from the 
ee7-vertex. To be more general we investigate the //7-vertex for arbitrary fermions /. 
The renormalized vertex function reads 



For on-shell external fermions it can be decomposed as {k — p' — p) 

{p' -P)n 



(3.25) 



2m/ 



2m/ 



(3.26) 

Expressing the renormalized quantities by the unrenormalized ones and the countcrtcrms 
and inserting this in the analogue of the renormalization condition (3.17) for arbitrary 
fermions we find, using the Gordon identities, 



= u{p)klll{p,p)u{p) 
= u{p)-f^u{p) [-Qf{6Ze + 8Z(^ + \5Zaa) + (0) + A^(0) + Vf\5ZzA 
-«(]?) 7m75m(p) [-Qf^Z(i^ + A;^(0) + af\5ZzA\ , 

where 



(3.27) 



(3.28) 
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and Vf,af are the vector and axialvector couplings of the Z-boson to the fermion /, given 
exphcitly in (A. 15). This yields in fact two conditions, namely 

= -QfiSZe + SZlt" + \5Zaa) + A^(0) + A^(0) + Vf]^5ZzA, (3.29) 

= -Qf5Z(i^ + A^(0) + af]^5ZzA. (3.30) 

The first one (3.29) for / = e fixes the charge renormahzation constant. The second 
(3.30) is automatically fulfilled due to a Ward identity which can be derived from the 
gauge invariance of the theory. The same Ward identity moreover yields 

A((0) + A^(0) - Q^5Z(i'' + afUzzA = 0. (3.31) 



Inserting this in (3.29) we finally find (using Vf — Uf — —QfSw/cw) 

Sw S^^(O) 



1 swi.^ idi:rik' 

oZf,^ = —-oZaa -oZzA ~ 



2 c^2 2 dk 



2 



k2^Q Cw Mz 



(3.32) 



This result is independent of the fermion species, reflecting electric charge universality. 
Clearly it does not depend on a specific choice of field renormahzation. Consequently the 
analogue of (3.17) holds for arbitrary fermions /. 

In the on-shell scheme the weak mixing angle is a derived quantity. Following Sirlin 
[26] we define it as 



sm^9w = sl^ = l--^, (3.33) 

Mz 

using the renormalized gauge boson masses. This definition is independent of a specific 
process and valid to all orders of perturbation theory. 

Since the dependent parameters sw and cw frequently appear, it is useful to introduce 
the corresponding counterterms 

cw,o ^cw + Scw, swfi ^ Sw + Ssw- (3.34) 

Because of (3.33) these are directly related to the counterterms to the gauge boson masses. 
To one-loop order we obtain 

Sew 1 fSM^ SMl\ ^ 1 ~ /^S^(M2,) S|^(M|)\ 



cw 2\M^ Ml 2 \ M"^ Ml y 

(3.35) 

Sw sl^ Cw 2sl. Ml )' 

We have now determined all renormahzation constants in terms of unrenormalized self 
energies. In the next sections we will describe the methods to calculate these self energies 
and more general diagrams at the one-loop level. 
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4 One-loop integrals 



Perturbative calculations at one-loop order involve integrals over the loop momentum. 
In this chapter we discuss their classification and techniques for their calculation. The 
methods described here are to a large extent based on the work of Passarino and Veltman 
[18], 't Hooft and Veltman [28], and Melrose [29]. 

4.1 Definitions 

The one-loop integrals in D dimensions are classified according to the number N of 
propagator factors in the denominator and the number P of integration momenta in the 
numerator. For P + D — 2N > these integrals arc UV- divergent. These divergencies are 
regularized by calculating the integrals in general dimensions D ^ A (dimensional regu- 
larization). The UV-divergencies drop out in renormalized quantities. For renormalizable 
theories we have P < N and thus a finite number of divergent integrals. 

We define the general one-loop tensor integral (see Fig. 4.1) as 



Tai...ap {Pi,---, Pn-1, mo, ... , miv-l) = 



with the denominator factors 



(27111) 



4-D 



I 



N-l 



(4.1) 



Do^ — ml + is, 



Di^ (q + Pi)^ -mf + ie, i^l,...,N-l, (4.2) 
originating from the propagators in the Feynman diagram. Furthermore we introduce 



Pio = Pi and pij ^Pi- Pj. 



(4.3) 



Evidently the tensor integrals are invariant under arbitrary permutations of the prop- 
agators Di, i ^ and totally symmetric in the Lorentz indices /ik- is an infinitesimal 
imaginary part which is needed to regulate singularities of the integrand. Its specific 
choice ensures causality. After integration it determines the correct imaginary parts of 
the logarithms and dilogarithms. The parameter n has mass dimension and serves to keep 



P21 



q + P2 



q + Pi 



q + PN- 



Pnn-1 



q + pn-2 

PN-IN-2 



Figure 4.1: Conventions for the N-point integral. 
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the dimension of the integrals fixed for varying D. Conventionally is denoted by the 
A'"th character of the alphabet, i.e. = A, T"^ = B, . . . , and the scalar integrals carry 
an index 0. 

Lorentz covariance of the integrals allows to decompose the tensor integrals into ten- 
sors constructed from the external momenta pi, and the metric tensor 5^,^ with totally 
symmetric coefficient functions T^^ We formally introduce an artificial momentum po 
in order to write the terms containing Qjj^^ in a compact way 

N-l 

^^...;up(Pi>--->Piv-i,mo, ...,mjv-i) = T^i...ipPhui---PiPi^P- (4-4) 

ii,...,ip=Q 

Prom this formula the correct g^j^i, terms are recovered by omitting all terms containing an 
odd number of po's and replacing products of even numbers of po's by the corresponding 
totally symmetric tensor constructed from the g^y, e.g. 



P0/ilP0/^2P0/i3P0/i4 ^ 9mfi2 9fJ.3fJ.4 ~^ 9filfJ.39fJ.2fi4 ~^ 9f^lf^4 9f^2f^3^ 
The explicit Lorentz decompositions for the lowest order integrals read 

Bn — PifiBi, 

-B^i/ = 9fji.vBoo + pi^piyBii, 

2 

C/x = PlfiCi + P2/xC2 = ^PifjCi, 

i=l 

Cf^u = 9fii^Coo + Pl^PluCii + P2fj.P2vC22 + {pinP2v + P2fi.Plu)Cl2 
2 

PifiPjvCij.1 



(4.5) 



(4.6) 



C^xi/p = {9u.vPlp + 9upPl,j. + 9u.pP\u)CoOl + {9fiuP2p + 9vpP2n + 9fipP2u)Co02 (4.7) 

+ PlfiPluPlpClu + P2nP2uP2pC222 

+ {PlfiPluP2p+PlnP2vP\p+P2nPluPlp)Cii2 

+ {.P2pP2uPlp + P2pPlyP2p + PlnP2uP2p)Cl22 
2 2 
= Y^i.9fiuPip + gvpPifi + 9fipPiu)CoOi + X] 

PifxPi vPkp Cij k 1 

1=1 i,j,k=l 
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3 

1=1 

3 

3 3 
Df^up = 'y^XOuvPip + QvpPip, + 9ppPiv)DoQi + ^ PipPjvPkpDijki 

(4.8) 

Dp,vpa — {g pug pa- + gppgva + gpagvp)DoQOQ 
3 

+ H {gpuPipPja + gupPipPja + gppPiuPja 

+ gpaPiuPjp + gvaPipPjp + gpaPipPjiy)DoOij 

3 

+ X] PipPjvPkpPlaDijkl- 

i,j,k,l=l 

Since the four dimensional space is spanned by four Lorentz vectors the terms involving 
g^i, should be omitted for iV > 5 and at most four Lorentz vectors should be used in the 
decomposition (4.4). Consequently the Lorentz decomposition for a general tensor integral 
with > 5 in four dimensions can be written as 

4 

Tli^,...pp{Pi,---,PN-i,mo,...,mN-i) ^ T^\...ipPiii^i- ■ -Pipi^P, (4-9) 

ii,...,ip=l 

where pi, . . . ,p4 is any set of four linear independent Lorentz vectors out of Pi, ■ ■ ■ ,Pn-i- 
The symmetry of the tensor integrals under exchange of the propagators yields relations 
between the scalar coefficient functions. Exchanging the arguments (pj,mi) {pj,mj) 
together with the corresponding indices i ^ j leaves the scalar coefficient functions in- 
variant 



T .t^... j...j...{Pi: • • • • • • • • • ,Piv-i,mo, ...,mi,...,mj,.. . ,miv-i) 

= T!^ ... j...j...iPu ...,Pj,...,Pi,.. .,PN-i,mo, ...,mj,...,mi,.. . ,mjv-i), 



(4.10) 



e.g. 

Ci{p-i,p2,mo,mi,m2) = 6*2(^2,^1, "^o, "^2, "7,i), 

Coo{pi,P2,mo,mi,m2) = Cooip2,Pi,mo,m2,mi), (4.11) 

Ci2ipi,P2,mo,mi,m2) = Ci2ip2,Pi,mQ,m2,mi). 

All one-loop tensor integrals can be reduced to the scalar ones T^^ . This is done in 
Sect. 4.2. General analytical results for the scalar integrals Aq, Bq, Cq are Dq are listed 
in Sect. 4.3. The scalar integrals for N > 4 can be expressed in terms of Dq^s in four 
dimensions. The relevant formulae are given in Sect. 4.4. They apply as well to the 
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tensor integrals with N < A in the kinematical regions, where the usual tensor integral 
reduction breaks down, because the Gram determinants appearing there are zero. The 
UV-divergent parts of the one- loop integrals are explicitly given in Sect. 4.5 

4.2 Reduction of tensor integrals to scalar integrals 

Using the Lorentz decomposition of the tensor integrals (4.4) the invariant functions 
Tf^ j can be iteratively reduced to the scalar integrals [18]. We derive the relevant 
formulae for the general tensor integral. 

The product of the integration momentum with an external momentum can be 
expressed in terms of the denominators 



QPk ^ ^[Dk - Do - fk], fk^pl-ml + ml. 



(4.12) 



Multiplying (4.1) with pk and substituting (4.12) yields 



R 



'N,k 

■/Ul-./Up-l 



rpN HP 

ix\...ixpirk 



1 (27r/x)^-^ 

2 i7r2 



/A 

• • ■ 



• • • 



Do... Dk-iDk+i ...D 



Di...D 



N-l 



Dn...D 



N-l 



pN-1 
-/ii.../ip_ 



N-l 

N 

//i...//P_i 



(4.13) 



where the argument k of the tensor integrals in the last line indicates that the propagator 
Dk was cancelled. Note that T^^~^^^__^{0) has an external momentum in its first propaga- 
tor. Therefore a shift of the integration momentum has to be performed in this integral 
in order to bring it to the form (4.1). All tensor integrals on the right-hand side of eq. 
(4.13) have one Lorentz index less than the original tensor integral. In two of them also 
one propagator is eliminated. 

For P > 2 we obtain one more relation by contracting (4.1) with g^^'^ and using 



= Do + ml. 



(4.14) 



This gives 



R 



N,00 

'/il.../ip_2 



7^^^ fiiJ-p-iiJ'P 

/il...^pi' 



(27r/i) 



4-D 



qi^i 



Di...D 



gMP-2 I 2^/^ 

+ 171 



N 



lXl...ilp-2 



° Do...D 



N 



(4.15) 



Inserting the Lorentz decomposition (4.4) for the tensor integrals T into (4.13) and (4.15) 
we obtain set of linear equations for the corresponding coefficient functions. This set 
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decomposes naturally into disjoint sets of A^" — 1 equations for each tensor integral. If the 
inverse of the matrix 



X 



N-l 



Pi 
P2P1 



P1P2 
Pi 



\Pn-iPi PN-1P2 



PiPn-1 
P2PN-1 

Pn-1 J 



(4.16) 



exists, these can be solved for the invariant functions T^^ yielding them in terms of 
invariant functions of tensor integrals with fewer indices (see eqs. (4.18) and (4.19) below). 
In this way all tensor integrals are expressed iteratively in terms of scalar integrals Tq 
with L < A^. 

If the matrix Xjv_i becomes singular, the reduction algorithm breaks down. If this 
is due to the linear dependence of the momenta we can leave out the linear dependent 
vectors of the set pi, . . . ,pn-i in the Lorentz decomposition resulting in a smaller matrix 
Xm- If Xm is nonsingular the reduction algorithm works again. This happens usually at 
the edge of phase space where some of the momenta pi become coUinear. 

If the determinant of X7V-1, the Gram determinant, is zero but the momenta are not 
linear dependent^ one has to use a different reduction algorithm [29, 30]. This will be 
discussed in Sect. 4.4. 

Here we give the results for the reduction of arbitrary N-point integrals depending 
on M < iV — 1 linear independent Lorentz vectors in D dimensions for nonsingular Xm- 
Inserting the Lorentz decomposition of Tn and R^''' as well as 



M 



TDN,k 

^ Hi... lip -I 



R 



Af,00 



rpN IJ'P 

fj,i...^ipl^k 



rpN llp_ilJ,p 
-^/il.../ipi' 



n,...,ip_i=0 
M 

ii,...,ip_2=0 



N,k 
ii...ip. 



.^Ph 



fJ-i 



■p. 



^p-ifJ-p-i 1 



(4.17) 



N,00 



Pilfj.1 ' ' ' P'. 



«P-2A'P-2' 



into the first hues of (4.13) and (4.15) these equations can be solved for the T{^_i^'- 

1 



Tr 



N 

00ii...ip_2 



D+P-2- M 



M 



R 



N,00 



'll...lp_2 



P-1 



.■lp-2 



k=l 



■^ii...ip_i / J ^ir OOii...ir-iir+i---ip-i 
r=l 



(4.18) 



Note that the numerator of the prefactor in the first equation is always positvc in the 
relevant cases P >2 and D > M. Using the third lines of (4.13) and (4.15) the i?'s can 
be expressed in terms of T^^ 



N rpN 

ip-1^ il...ip_2' 



, and Tj^../ , with q < P as follows 



R: 



N,00 

'h...iq M...M 



2rpN 



Zrrf 



P-2-q 



+ (-1) 



1 



fci=i 



^This can happen, because of the indefinite metric of space time. 
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+ 



P-2-q 



M-l 
ki,k2=l 



V ^/ fel,...,fep_2- 



Af-l 



.(0) 



■n...i„ M...M 



P-l-q 




(4.19) 



P-l-g 



P-l-q 



p-l-q 



-(-1) 



P-l-q 



p-l-q 



M-l 

/ ^ -"^ ii...iqki (^) 
fei=l 



+ 



X] ^ii...igfcifc2(0) + • ■ ■ 



\ ^/ fcl,...,fep_l_q = l 



Af-1 

h...iqkl...kp-l-q 



where ii, . . . , L 7^ M and 



6'(A; I ii, . . . ,ip-i) 



1 ir 7^ k, r — 1, . . . , P — 1, 
else. 



(4.20) 



The indices i refer to the i-th momentum of the corresponding iV-point function but 
to the (i — l)-th momentum of the N — 1-point function T^~^{k) if i > k. Again the 
arguments of the T's indicate the cancelled propagators. The tilde in T(0) means that a 
shift of the integration variable q ^ q — Pm has been performed in order to obtain the 
standard form of these integrals. This shift generates the terms in the square brackets 
of (4.19). It is also the reason for the unsymmetric appearance of the index M in the 
above equations. A different shift would result in similar results. An explicit example 
illustrating the use of these reduction formulae is given in App. C. 

The recursion formulae above determine the coefficients Tj^...jp regardless of their sym- 
metries. Consequently coefficients whose indices are not all equal are obtained in different 
ways. This allows for checks on the analytical results as well as on numerical stability. 

If the number M of linear independent momenta equals the dimension D of space-time 
then the terms containing g^^, in the Lorentz decomposition have to be omitted, since g^i, 
can be built up from the D momenta. In this case the coefficients are obtained 



from the second equations in (4.18) and (4.19) with T^^^ 



■lP-2 



= 0. 



4.3 Scalar one-loop integrals for N <A 



With the methods described in the last section all one-loop integrals can be reduced to 
the scalar ones Tj^provided the matrices Xm are nonsingular. General analytical results 
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for Ao,Bq, Co and Dq were derived in [28]. Algorithms for the numerical calculation of the 
scalar one-loop integrals based on these results have been presented in [31]. Here we give 
a new formula [32] for Dq involving only 16 dilogarithms compared to 24 of the solution 
of [28]. For completeness we first list the results for Aq, Bq and Cq. 

4.3.1 Scalar one-point function 

The scalar one-point function reads 

Ao{m) = -m2(^)^r(l - ^) = m'{A - log — + 1) + 0{D - 4), (4.21) 
with the UV-divergence contained in 

2 



-7E + log47r (4.22) 



and 7e is Euler's constant. The terms of order 0{D — 4) are only relevant for two- or 
higher-loop calculations. 

4.3.2 Scalar two-point function 

The two-point function is given by 

B.(p.o, m„, = A - /' log \p\«-' - -A + ^e] ^o(D-i) 

Jo IjL^ 

momi ml -ml mi m^mi 1 

= A 2 - log \ — - log 5 — ( r) logr 

1^ Pw "lo Pio r 

+0{D-i), 

(4.23) 

where r and - are determined from 

r 

^ ml + ml - pIq - is ^^^^^^^^^^^^ly ^424) 
m^mi r 

The variable r never crosses the negative real axis even for complex physical masses (m^ 
has a negative imaginary part!). For r < the ie prescription yields Imr = £:sgn(r — ^)- 
Consequently the result (4.23) is valid for arbitrary physical parameters. 

For the field renormalization constants we need the derivative of Bq with respect to 
pIq. This is easily obtained by differentiating the above result 

d ml -ml mi momi 1 

■^-^Bo{pio, mo, mi) = ^ log \ —{- - r) logr 

dpio Pio rno ph r ^^^^^ 

-^(l + 4±|logr)+0(D-4). 
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4.3.3 Scalar three-point function 



The general result for the scalar three-point function valid for all real momenta and 
physical masses was calculated by [28]. It can be brought into the symmetric form 

CQ{pio,P2o,mo,mi,m2) = 
Jo Jo 

+ {ml -ml- pI^)x + [ml - mj + pl^ - ^20)2/ + ml - ie]'^ (4.26) 



a 



1=0 \a=± 



+77(1 - Xia, — ) log 



Vi. 



yoi - 1 / 1 \ , yoi 
- — log — 



yi< 



Vi. 



ri{-Xi+, -Xi-) - r]{yi+,yi-) - 2me{-p]k)0{ - lm{yi+yiJ)) 



log 



1 - Vio 
-Vio 



with (i, j, A; = 0, 1, 2 and cyclic) 

1 



yoi 



2 2 \P%(P% - Pki - pI + - - ml) 



-{Pli - P%){m] - ml) + a{p% - m] + ml)], 



'^Pjk 



and K is the Kallen function 



yi± = voi - Xi±, 

= /«(PlO>P21>P2o)' 

«i = «^(P?fe> m], ml) (1 + iep]k), 



(4.27) 



k(x, y, z) = \J x'^ + y"^ + — 2(xy + yz + zx) . 
The dilogarithm or Spence function hi2{x) is defined as 

"1 dt 



Li2(x) = - / — log(l-a;i), |arg (1 - a;)| < tt. 
Jo t 



(4.28) 



(4.29) 



The 77-function compensates for cut crossings on the Riemann-sheet of the logarithms and 
dilogarithms. For a, 6 on the first Riemann sheet it is defined by 



log(a6) = log(a) + log(6) + r){a, b). 



(4.30) 



All 77-functions in (4.26) vanish if a and all the masses are real. Note that a is real in 
particular for all on-shell decay and scattering processes. 
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4.3.4 Scalar four-point function 

The scalar four-point function Dq{pio,P2o,P3o, mo, mi, 1712,1713) can be expressed in 
terms of 16 dilogarithms [32] . 

Before we give the result we first introduce some useful variables and functions. We 
define 



mj + m] - 



ij 



mimj 



, = 0,1,2,3, 



and rij and fij by 
and 



x"^ + kijX + 1 = (x + rij){x + I /rij). 



(4.31) 
(4.32) 



x^ + {hj - i£)x + 1 = (x + rij){x + I /rij). (4.33) 
Note that for real kij the rjj's lie either on the real axis or on the complex unit circle. 



Furthermore 



P(yo,yi,y2,y3) ^ Yl ^jyiyj + Yy]^ 

0<i<j<3 j=0 



P(l, —, 0, 0) - is] \P(0, 0, — , x) - is 
X I'- ri3 -IL ro2 

X 1 n r 

P(0, — , — ,Q)-is P(l,Q,Q,x)-is 
ri3 ro2 



ax^ + bx + c + isd = a{x — xi){x — X2), 



where 



a = k23/ri3 + ro2koi - A;o3?^02/ri3 - ki2, 

b = (ri3 - l/ri3)(ro2 - l/ro2) + koik23 - /co3^i2, 

c = /coi/ro2 + ri3/c23 - ko3ri3/ro2 - k^, 

d = ki2- ro2koi - ri3k23 + ro2ri3A;o3- 
In addition we introduce 

-fki ^ sgaRe[a{xk - xi)], k,l^ 1,2, 

and 

XkO = Xk, So = fo3, 

Xki = Xk/ri3, Si = foi, 

Xk2 = Xkro2/ri3, S2 = fi2, 
Xk3 = Xkro2, S3 = f23- 



(4.34) 



Qiyo,yi,0,y3) = (l/ro2 - ^02)1/0 + (/^i2 - ^02/^01)2/1 + (/^23 - ^02/^03)2/3, (4.35) 

Q{yo, 0, y2, 2/3) = (1/^3 - ^13)2/3 + (^^12 - ri3k23)y2 + {koi - ri3ko3)yo- (4.36) 
and xi^2 is defined by 

ro2ri3 



(4.37) 



(4.38) 
(4.39) 



(4.40) 
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as well as 
Finally we need 

fj{a,b) 



X 



(0) 
kj 



limxfej as Tij — limf- 



I]- 



ri{a, h) for h not real, 

2m [^(-Ima)^(-Im6) - ^(Ima)^(Im6)] for 6 < 0, 
for 6 > 



with h — linig^o ^• 

Then the result for real ro2 can be written as 

-Do(PlO,P20,P30, mo, mi, 1712,1713) 



mim2m2,m4^a{xi — X2) 



3 2 



I IZ(-l)^^^ Li2(l + SjXkj) + ri{-Xkj, Sj) log(l + SjXkj) 

I j=0 k=l '- 

+ Li2(l + ^) + r;(-a;.„ -) log(l + 



+ E(-i)' 
fe=i 



fc+i 



'7(-a:fe,ro2) 



1 



log(ro2a;fc) + log(Q(-^, 0, 0, 1) - ie 



+ log( 



g(0,0,l,ro24 



(0)^ 



-/ 1 

ri3 



d 



+ i£7fc,3-feSgn(ro2lmri3 



log(^)+log(Q(^, 1,0,0) -is) 

^13 X^ 



, /Q(l,0,0,4°^) . ,^ _ ,V 

+ log( ^^ + z£7fe,3-fcSgn(Imri3)) 



?7(-a;fe, +^(ro2, — ) 

(0)^ 



^^13 



+ log( 



g(0,0,l,ro2a;r) 



+ i£7fe,3_fcSgn(ro2lmfi3)) 



+ Tl{fQ2,^)Ti{-Xk,-^-^) 

ri3 



In the case that |rj,| = 1 for all rj,, the result reads: 



-Do(PlO,P20,P30, "^0,^1,777,2, ma) = 



1 



I j=0 fc=l 



mim2m3m4a{xi — X2) 
Li2(l + SjXfcj) + r]{-Xkj, Sj) log(l + SjXfcj) 

+ Li2(l + ^)+77(-x,,,l)log(l + i^; 



(0) (0) 

log(^P(l, ^, 0, 0) - ^£fo7.,3-0 + log( 

^13 '^13 ^ ^^13 



+ E(-i)'^' 



fe=l 



?^13 
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+ v{-Xk,ro2) 



(0)^ 



77(-Xfc, — ) + r7(ro2, — ) 
ri3 ri3 . 



(0) (0) 



ri3 



+ log( 



(0) 

ri3 



+ (l - 7fe,3-feSgn(6))77(-a:^ifc, -— )?7('^02, — ) [• 



£ is understood as infinitesimally small. 

4.4 Reduction of scalar and tensor integrals for vanishing Gram determinant 

Using the four- dimensionality of space-time the scalar five-point function can be re- 
duced to five scalar four-point functions [29, 31]. Furthermore if the Gram determinant 
of the external momenta of a tensor integral vanishes, 



\X 



N-l\ 



Pi 
P2P1 



P1P2 

pI 



Pn-iPi PN-1P2 



PiPn~i 

P2PN~1 



P%-i 



0, 



(4.44) 



this tensor integral can be expressed in terms of N integrals ^. This is always the 
case for N > 6, because any five momenta are linear dependent in four dimensions. 

4.4.1 Reduction of scalar five-point functions 

Here we assume that the four external momenta appearing in the five-point function 
span the whole four- dimensional spacc^. Then the integration momentum q depends 
linearly on these four external momenta and the following equations holds 



2g2 2gpi 
2p,q 2pl 

2p4q 2p4Pi 



2qpi 
2piP4 



2L>o + ^^00 2gpi 
L>i-L>o + Yio-yoo 2p2 

-D4 - -Do + Yao - Yoo 2pipi 



2qpi 

2piP4 

'2pl 



(4.45) 



with 



Yij = ml + m]- {pi - pjf. 



(4.46) 



and Di as defined in (4.2). Thus we have 



1 



2Do + Yoo 2qp^ 
D,-Do + Y,o-Yoo 2pl 

-D4 - Do I40 - ^00 2p4Pi 



2?P4 

2piP4 



2pl 



(4.47) 



^The exceptional case, when they are hnear dependent will be covered in the next section. 
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Expanding the determinant along the first column we obtain 



2pipi . . . 2pip4 

2T^{Q) + YooT^] ; ••. ; 

2p4Pi . . . 2p4P4 



k=l 



X 



2pi^ . 

2pipi . 

2pk-iPi ■ 

2pk+ipi . 



2piPA 

2pk-iPi 
2pk+iP4 



(4.48) 



2p4pi ... 2pl 

where the arguments of the functions denote again the cancelled propagators. 

The Lorentz decomposition of the vector integrals in (4.48) involves only the momenta 
Pi,...,P4. T^fj,^ ^^{k) does not depend on p^, consequently each term in its Lorentz 
decomposition contains a factor p^^, i ^ k and its contraction with the corresponding 
determinant in (4.48) vanishes. Similarly all terms in the tensor integral decomposition of 
T^(0) + P4^Tq (0) involve a factor — p4^, i = 1, 2, 3, if one performs the shift q ^ q — pi 
to bring the tensor integral to the standard form. Multiplying with the determinants 
and performing the sum in (4.48) these terms drop out. Finally in the term PAjxTQiQi) 
the determinant is nonzero only for k — A where it can be combined with the first term 
in (4.48). Rewriting the resulting equation as a determinant and reinserting the explicit 
form of the tensor integrals we find 



1 



Do^i ■■■Da 



Do + Yoo 2qpi 



'00 



Using 



^40 - ^00 2p4Pl 



2piPj = - Yio - Yoj + Yqo, 
2qpj = Dj -Do + Yoj - Foo, 



2qp4 
2piP4 

2pl 



0. 



(4.49) 



(4.50) 
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adding the first column to each of the other columns and then enlarging the determinant 
by one column and one row this can be written as 



^00 

T,\0) + YooTi 
Yio — Yqo 



Yao — Yc 



00 



^04 

T4(4)+F4oTo^ 

^14 ~ ^04 



I44 - Yc 



04 



(4.51) 



This is equivalent to 



7^5 
-^0 


-To^(O) 


-To\l) 


-To\2) 






1 


^00 


>01 


Y02 


^03 


Yo4 


1 


Yio 


Yn 


Y12 


1^13 


Yu 


1 


Y20 


Y21 


Y22 


Y23 


Y24 


1 


Yso 


Y31 


Y32 


Ys-s 


Yu 


1 


I40 


I41 


Y42 


Y43 


Yu 



0, 



(4.52) 



which can be solved for Tn if the determinant of the matrix l^;,, i,j 



0, 



4 is nonzero. 



Note that in the integral Tq (0) the momenta have not been shifted. In particular (4.52) 
yields the scalar five-point function Tq in terms of five scalar four-point functions. 

4.4.2 Reduction of N-point functions for zero Gram determinant 

For vanishing Gram determinant |Xjv-i| the following relation holds, if the Lorentz 
decomposition of the appearing tensor integrals contains only momenta and no metric 
tensors, which is the case for A?" > 6 or P = (scalar integrals) 



(27r/i) 



A-D 



I 



DoDi ■ ■ ■ Dn-1 



Do + Yoo 
Yio — loo 



2qpi . . . 2qpN-i 
2pj ... 2piPN_i 



Fjv-io - Yqo 2pN-iPi . . . 2p%_^ 



0. 



(4.53) 



Performing the same manipulations of the determinant as in (4.49) to (4.52) above this 
results in 

1 Yqq Yqi 

1 YiQ Yii 



Yqn-i 

YiN-l 



Yj 



N-lN-l 



0, 



(4.54) 



1 ^Af-lO IjV-ll 

valid for |Xjv_i| = and > 6 or P = 0. We stress again that in the tensor integral 
/ip(0) appearing in (4.54) the momenta have not been shifted. Eq. (4.54) determines 



rnN 
Ml . - 

„ in terms of T^-\ 

Ml •••MP Ml •••MP 



0, 



,A^- 1, 



if the determinant of the matrix Yij is 
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nonzero. The vanishing of this determinant corresponds to the leading Landau singularity 
of which is clearly not contained in T^~^. In this case one has to calculate directly 



Eq. (4.54) in particular expresses by T^'K For = 5 and P = (4.54) coincides 
with (4.52), which is thus vahd for arbitrary momenta. For N > 6 one can choose any 
six out of the N denominator factors resulting in different reductions. For < 4, where 
(4.54) is only valid for scalar integrals, the Gram determinant is singular at the edge of 
phase space where some of the momenta pi become collinear, i.e. for forward or backward 
scattering or at the threshold of a certain process. Because in this special situations all 
integrals can be reduced to lower ones one can obtain considerably simpler formulae than 
in the general case (see e.g. [33]). 

With the methods described in this section all tensor integrals with N > 6 can be 
reduced directly to tensor integrals with smaller A^. Note that this may yield tensor 
integrals with P > N because P is not reduced simultaneously as in the reduction method 
described in Sect. 4.2. These tensor integrals are not directly present in renormahzable 
theories. Nevertheless their reduction to scalar integrals can be done with the formulae 
given in Sect. 4.2. 

4.5 UV-divergent parts of tensor integrals 

For practical calculations it is useful to know the UV-divergent parts of the tensor 
integrals explicitly. We give directly the products of D — 4 with all divergent one-loop 
tensor coefficient integrals appearing in renormahzable theories up to terms of the order 



[30]. 



0{D-i) 



{D 



4)^o(m) 

4) BQ{piQ,mQ,mi) 
4) 5i(pio,mo,mi) 



-2m2, 



(D 

(D 



4) -Boo(pio,mo,mi) 

4) -Bii(pio,mo,mi) 

4) Coo{pio,P2o,mo,mi,m2) 

4) Cooi(pio,P2o, mo, mi, ma) 

4) -Doooo (Pio , P20 , P30 , mo , mi , ma , ma) 



1, 

l{plQ-3ml-3ml), 



(4.55) 



2 

3' 



1 

2' 



1 

6' 



12" 



All other scalar coefficients defined in (4.7) and (4.8) are UV-finite. 
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5 Standard matrix elements 



5.1 Definition 

The invariant matrix elements for scattering and decay processes involving external 
fermions and/or vector bosons depend on the polarizations (Tj, a[ and Aj of these particles. 
This dependence is completely contained in the polarization vectors Sfj^X^i^ -^i) ^^'^ spinors 
Vaiip'i, cr^) and Uai{Pi, CTj). ki, p[ and Pi denote the incoming momenta of the vector bosons, 
antifermions and fermions, respectively. For outgoing fermions one has to replace p, pi 
by —Pi —p' and one must use u{—p) — v{p). If we split off the polarization vectors and 
spinors from the invariant matrix element Ai we are left with a tensor involving Lorentz 
and Dirac indices in the general case 

M = (y'l)--- VaAP'n, 0-^ai:.:ar/3i.../3n"/3i(^^l' ^0 " " " ^PniPn, (^n) 

XSfiliki, Ai) . . -S^i^ikm, Am)- 

To be definite we choose m external vector bosons and n external fermion-antifermion 
pairs. The tensor M.^l'^^'p^ can be decomposed into a set of covariant operators together 
with the corresponding scalar formfactors Fj 

MZ::t^^Mt^^-^%F,. (5.2) 

i 

We call the covariant operators J^ia''% multiplied by the corresponding polarization 
vectors and spinors standard matrix elements /Ai 

Mi = Va, (p'l, (t[)... Vo^Sv'n^ OMf^^-;%U/3,(pi, (7i) . . . li/3„(Pn, (7„) 

In this way the invariant amplitude M. is decomposed into polarization independent 
formfactors Fi and the standard matrix elements Mi 

M^Y^MiFi. (5.4) 

i 

The formfactors Fi are complicated model dependent functions involving in general the 
invariant integrals and the counterterms. The standard matrix elements in contrast 
are simple model independent expressions which depend on the external particles only 
but contain the whole information on their polarization. They are purely kinematical 
objects. All of the dynamical information is contained in the formfactors. 

The covariant tensor operators forming the standard matrix elements can be con- 
structed from the external four- momenta p^, p'i and ki, the Lorentz tensors g^^ and £^^'^P'^ 
and the Dirac matrices 7'',75. In general one thus obtains an overcomplete set. Dirac 
algebra and momentum conservation are used to eliminate superfluous operators. Since 
the external particles are on-shell, the Dirac equation for the fermion spinors 

^iu{pi, ai) = miu{pi, ai), 
v{p'i,(y'i)fi = -m'iv{p'i,a'i) 
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and the transversality condition for the polarization vectors 

kt'e^,iki,\i)=0 (5.6) 

reduce the number of independent standard matrix elements further. 

The number of independent standard matrix elements cannot be larger than the num- 
ber of independent polarization combinations of the external particles. In four dimensions 
there are only four linear independent four-vectors. Expressing all four-vectors in a def- 
inite basis allows to derive the missing relations between the remaining standard matrix 
elements. Thus a minimal set of standard matrix elements can be constructed. 

If there are only few external particles there may be less independent standard ma- 
trix elements than different polarization combinations, since there are only few momenta 
available for their construction. In this case some of the polarized amplitudes are related. 

The number of standard matrix elements can be reduced further if the model under 
consideration exhibits certain symmetries. These evidently also apply to the relevant 
standard matrix elements. 

For many applications it is not essential to minimize the number of standard matrix 
elements. All one needs is a complete set. 

Furthermore the choice of the standard matrix elements is not unique. This allows to 
arrange for the most convenient set according to simplicity, the structure of the lowest 
order amplitudes and, if present, symmetries. At least some of the formfactors can be 
chosen as generalizations of the lowest order couplings. This is useful in estabhshing 
improved Born approximations. 

The concept of standard matrix elements is not indispensable for the calculation of 
amplitudes in higher orders. It is, however, extremely helpful in organizing lengthy cal- 
culations, which often are inevitable. All comphcated expressions are cast into the form- 
factors which are polarization independent and thus have to be evaluated only once. 

An alternative method would be to calculate directly the polarized amplitudes. This 
requires a definite representation for the spinors and/or polarization vectors from the 
start. The whole calculation has to be done for each polarization separately. A closer 
look shows that this method can be represented as a particular case of the standard 
matrix element approach. The corresponding covariant operators are constructed from 
the polarization vectors and spinors instead of the momenta, Lorentz tensors and Dirac 
matrices. Their explicit form is 

'■i,ai...anPi--Pn V 2m[ ' ' ' 2m'„ 2mi ' ' ' 2m„ 

where m, m' are the masses of the external spinors. The indices i correspond to different 
polarization combinations. Consequently the number of different standard matrix ele- 
ments equals the number of polarizations of the external particles. For each polarization 
only one standard matrix element is nonzero. In this sense the set of standard matrix 
elements (5.7) is orthogonal. The formfactors equal the polarization amplitudes and are 
directly obtained by inserting explicitly the polarization vectors and spinors in the in- 
variant matrix element. Unlike in the approach outlined above these formfactors are no 
direct generalizations of the lowest order couplings. 

In the following we list complete sets of standard matrix elements relevant for the 
production of bosons in fermion-antifermion annihilation. 
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5.2 Standard matrix elements for processes with two external fermions 



In this section we will give the standard matrix elements for processes involving two 
external fermions (FF) and one [or two] scalar (S) or vector (V) bosons. The momenta 
and spinors of the fermions are denoted by pi, p2 and v{pi) — v{pi, ai), u{p2) — u{p2, (72), 
the momenta and polarization vectors of the bosons by ki, ei = e{ki, Ai) [and k2, £2 ~ 
e{k2, A2)]. The numbers of different polarizations for each scalar, fermion and vector boson 
are 1, 2 and 3, respectively. If we use momentum conservation to eliminate ki [or ki + ^2] 
the standard matrix elements are constructed from the momenta p2 and pi [and ki — /C2] , 
the polarization vectors of the vector bosons, the totally antisymmetric tensor £^^''P'^ and 
Dirac matrices between the spinors. If there are products of £-tensors, pairs of them can 
be eliminated using 



gfMX glll3 gll-y g^^ 

gUa gUP gUy g^ S 

gpa gP/3 gPy g pS 

^cra g(yl^ Q'^^ 



(5.8) 



If any of the left over ^-tensors are contracted with four four-momenta, we write for one 
of these momenta = |{j^,7"} between the spinors. Now all remaining e-tensors are 
contracted with one 7-matrix at least and can be eliminated using the Chisholm identity^ 

e^^vpal" = -ibiilvlp - Qi^vlp + Qi^plv - gvpl^\lb- (5.9) 

All Dirac matrices contracted with pi or p2 can be eliminated using Dirac algebra and 
the Dirac equation. Consequently the only remaining Dirac matrices are contracted with 
polarization vectors [and If^i — lf^2\ and there is at most one of each type. Finally in the 
scalar products involving the polarization vectors only one [or two] independent momenta 
may appear because of transversality and momentum conservation. 

Thus we arrive at the following sets of standard matrix elements (we suppress polar- 
ization indices in the following): 

5.2.1 S ^ FF 

There are 2x2 = 4 different polarization combinations but only two standard matrix 
elements 

= u{pi)oj,v{p2), (5.10) 
where a = ± and lj± = and the fermions are outgoing. 

5.2.2 V ^ FF 

Replacing the scalar by a vector results in 3 x 2 x 2 = 12 different polarizations and 
yet only four standard matrix elements 

Ml = U{pi) ^lUJaV{p2), ^^^^^ 
M2 = u{Pl) (^a V{P2) SiPi. 



"'^Eq. (5.8) and (5.9) can be applied because the standard matrix elements involve only external vectors 
and spinors which remain four-dimensional also in dimensional regularization. 
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5.2.3 FF SS 



Here the number of independent polarizations four equals the number of standard 
matrix elements 

= ^Xpi) |(I^-i-^2) u;^ u{p2), 

M2 = V{pi) Ua U{p2). 



(5.12) 



5.2.4 FF SV 

In this case we find twelve standard matrix elements for 2x2x3 
polarizations 

= v{pi) ^2 u{p2), 

M2 = V{pi) - ^2) UJa u{p2)e2Pl, 

Ml = V{pi) - 1^2) u{p2)e2P2, 

Ml = V{pi) ^2^2 Co-cr m(P2), 



12 different 



(5.13) 



Ml, 

MI2 
Ml, 



5.2.5 FF ^ VV 

There arc 2x2x3x3 = 36 different polarization combinations, however, we can 
construct 40 standard matrix elements 



Ml = v{pi) ^1 (^1 -^1) ^2 u{p2), 

M2 = V{pi) liltl - 1/1:2) {£182) UJa U{p2), 
V{pi) ^1 U{p2) {S2kl), 

v{pi) ^2 u{p2) {eik2), 

V{pi) ^1 UJ^ U{p2) {€2P2), 
V{pi) ^2 U{p2) {eiPl), 

- h) u{p2) (sih) (£2^1), 

- h) ^(P2) (siPi) (£2^2), 

1(^1 - ^2) i^<T u{p2) {eik2) (£2^2), 

K^l - ^2) C^fT U{p2) {SiPi) (£2^1), 
^1 ^2 1i(P2), 

^i(P2) (£1^2), 

^1 C^^ m(p2) (£2^:1), 
1^2 h «(P2) (£1^:2), 
^ll^lUa U{p2) (£2^2), 



-MI2 



— ^VPl 



v{p 
Ml = v{pi 



Mil = v{p 
v{p 
v{p 



M12 

Mil 
Ml 



,12 - v{p 

•^13,1 ^^iV 



Ml,^2 



v{pi 

Ml^^i = vipi 



(5.14) 
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•^14,2 = ^{Pl) h h U{P2) (SlPl), 
•^15 = ^(Pl) '^'T U{p2) (£1^2) (£^2^:1), 

MIq = v{pi) u{p2) {eipi) (£2^2), 

•^17,1 = ^(Pl) U{P2) (£1^2) (£2^2), 
•^17,2 = HPi) Uip2) (SlPl) i£2kl)- 

We have obtained more than 36 standard matrix elements because we have not yet used 
the four dimensionahty of space time, i.e. the fact that the five vectors pi, p2, ki — k2, £1 
and £2 are hnear dependent. The relations between the 40 standard matrix elements can 
be found by representing these vectors in a certain basis using for example Vi = P1+P2, 
V2 — Pi — P2, V3 — ki — k2, V4,/i = £ni'pa'ViV2V^. In this way one can derive the relation 

= 2(pip2 - mima) {M^, + 

— 2 (^2^:2 — ml— 17111112) M'^i — 2 (pi/ci — ml — 11111712) M'^^2 

-2{kj + {kik2)) Ml, - 2 {kl + {kik2)) MI2 - 2MI + 2{Ml, + MI2) 

- 2{mi{ml - P2k2) + m2(m^ - piki)){Ml, - MI2) (5.15) 
+ (mi + m2){piP2 - rnim2) MI2 - 2 777,2 -A^is,! - 2 mi Ali3_2 

+ (mi + m2) (2M^4,i + 2M'[^^2 - ^^is - ^M^^) 
+ (3mi + m2) Ali7,i + (3m2 + mi) M'[-j^2 

and a similar independent one allowing to eliminate four of the 40 standard matrix ele- 
ments (5.14). 

The construction of complete sets of standard matrix elements described above is 
straightforward. The reduction of general structures to these standard matrix elements 
is therefore easy to implement into computer algebra programs. In practical applications 
some of the standard matrix elements may not contribute due to the presence of symme- 
tries and/or the neglection of fermion masses. These aspects will be discussed together 
with the applications in the following chapters. 



5.3 Calculation of standard matrix elements 



For the calculation of the standard matrix elements one has to choose a certain rep- 
resentation for the polarization vectors and spinors. This has to be done only once for 
each process and not in the calculation of individual Feynman diagrams. If there are 
at least four external particles the polarization vectors can be constructed from their 
four-momenta respecting 

ki • £i{ki, Aj) = 0, 



(5.16) 

£i{ki,Xi)ei{ki,\'^) = -SxiX- 



We thus obtain for S2 
e'2{k2A\) 



\l\piP2{2pik2P2k2 - klpip2) +plplkl -pl{pik2)^ - pi{P2k2y 
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X 



X 



P2 +P2) kl - Pik2 (P2k2 +Plk2)) 



- Pi (P2(P1 + P2) kl - P2k2 (P2k2 + P1/C2)) 
+ ^2 (P1P2 (Plk2 - P2k2) +Plpik2 - p\p2k2) 



= (O,cosi?,0, -sini?) 



(5.17) 



\JpiP2{'ipik2P2k2 - klpip2) ^PiPlkl -P2{Pik2y -Pi{P2k2y 
(0,0,1,0) 



(5.18) 



4(^2, L) 



\lkl[{p2k2 +Pik2Y - kl{p2 + piy 
{k, E2 sin 0, cos i?)/ -v/^ . 



/C2 (P2/C2 + P1/C2) - {Pi + P2)''A;2 



(5.19) 



where we have also given the simple expressions in the CMS-system of the fermions and 
bosons. In this system the four-momenta of the external particles read 



P1.2 = (^i,2,0,0,T|p|), 

/ci,2 = (-Ei,2, =F|k| sin'!?, 0, =F|k| cost?). 



(5.20) 



£"12 are the energies and p the three-momentum of the fermions and £^1^2 the energies 
and k the three-momentum of the bosons, 'd is the angle between the spatial vectors p 
and k. 

From the polarization vectors given above the ones for helicity states are obtained as 
£^(A;2, ±) = ^ [4{k2, II) ± ^4(^2, ±)] , 4(A;2, 0) = £^(A;2, L). (5.21) 

The polarization vector ei can be obtained by interchanging 1^2. 

For the case of only three external particles one needs a further independent vector. It 
can be chosen freely but linear independent of the momenta. Using this additional vector 
as one of the polarization vectors the others can be constructed using (5.16). 

Inserting the polarization vectors (5.19) into the standard matrix elements these can 
be reduced to the ones for external scalars, i.e. to (5.10) for the decay V — > FF, and to 
(5.12) for the annihilation processes FF VV. VS. To calculate these remaining Dirac 
matrix elements one either inserts a definite representation for the spinors or evaluates the 
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quantities M.f*M.f via traces and reconstructs M.^ from those if needed. Note that for 
the calculation of \M.\'^ to one-loop order one only has to evalulate the products A^^*Alj' 
for those values of i, where Ff is nonzero in lowest order 

|A<|2 = \Mo + SMi\^ ^ \Mo\^ + 2Rc{M'o6Mi} 

Here Aio: F[o denote the lowest order quantities and 5Aii, SF^^ the one-loop quantities. 

For massless external fermions the Dirac matrix elements (5.10) and (5.12) are equiva- 
lent to the helicity matrix elements. They do not interfere and can thus easily be obtained 
from lAlf P as 



V{pi) U{p2) = \/2PlP2, 

(5.23) 

v{pi) (^1 - ^2) <^(T u{p2) = y4:pi{ki - k2)p2{ki - k2) - 2pip2{ki - k2Y. 

If one is only interested in unpolarized quantities it suffices to calculate Y^poi ■Mi*M.j 
using the polarization sums for vector bosons and spinors. 
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6 Calculation of one-loop amplitudes 

We have described all the ingredients necessary for the calculation of one-loop radiative 
corrections. This chapter shows how one-loop amplitudes are evaluated in practice. 

First one has to specify a Lagrangian and to derive the corresponding Feynman rules. 
Then renormalization has to be carried out and the counterterms have to be determined. 
Both were done in Chap. 2 and 3 for the SM. Once the Feynman rules and the counterterms 
are fixed, the following steps apply to any renormalizable model. 

To calculate the amplitude of a certain process at the one-loop level one has to con- 
struct all tree and one-loop Feynman diagrams with the given external particles allowed 
by the specified model. Next each Feynman diagram has to be reduced algebraically to 
a form suitable for numerical evaluation. This procedure is explained in more detail in 
Sect. 6.1. Finally the expressions for all diagrams have to be put together into a numeri- 
cal program which calculates the amplitude and the corresponding cross section or decay 
rate. 

6.1 Algebraic reduction of Feynman diagrams 

The algorithm for the reduction of one-loop diagrams is the following. The loop inte- 
gral obtained from the Feynman rules contains a product of propagators as denominator 
and a numerator composed of Lorentz vectors and tensors, Dirac matrices and spinors and 
polarization vectors of the external particles. The numerator is simplified using tensor 
and Dirac algebra, the mass shell conditions for the external particles and momentum 
conservation. One can also try to separate terms proportional to one or more of the 
denominators. CanceUing these yields N-point functions of lower degree. Next the loop 
integral is organized into the tensor integrals defined in Sect. 4.1. The Lorentz decompo- 
sition of these integrals is inserted and the whole Dirac and Lorentz structure is separated 
off from the integrals. Using again Dirac algebra and mass shell conditions it can be 
reduced to the appropriate standard matrix elements as discussed in Chap. 5. 

We thus arrive at an expression of the form 

SM^Y^^i^^i (6-1) 

i 

for each one-loop Feynman diagram. The formfactors are linear combinations of the 
invariant coefficient functions of the tensor integrals with coefficients being functions of 
the kinematical invariants. 

The formfactors can be further evaluated by applying the reduction scheme for the 
invariant integrals described in Sects. 4.2 and 4.4. Finally they are obtained as linear 
combinations of the scalar one-loop integrals Aq, Bq, Cq and Dq which are given explicitly 
in Sect. 4.3. This last step may lead to very lengthy expressions. Their algebraic evalu- 
ation needs a lot of time and space. This can be avoided by performing the reduction to 
scalar integrals numerically. 

The evaluation of the counterterm diagrams and the Born diagrams is done in a similar 
way. Since no integrations have to be performed their calculation is much easier. In most 
cases the counterterm diagrams can be obtained from the Born diagrams by replacing the 
lowest order couplings by the corresponding counterterm. 
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Figure 6.1: Box diagram contributing to e'^e — > W'^W' 



As an illustration of the reduction method we present the explicit calculation of a box 
diagram contributing to e+e^ W^W^ (Fig. 6.1). According to the Feynman rules the 
corresponding contribution to the invariant matrix element 6Ai is given by (we include a 
global factor i in the Feynman rules in order to obtain real amplitudes) 



SM = 



p4 2 



9s4 



(6.2) 



with 



{27i)D [g2 _ M|] [{q + hy - M^] {q + h- p^f [{q - - M^] ' 



(6.3) 

^Xai^ = gxai-h - q)u + gav{2k2 - q)\ + Qvxi'lq - k2)a ■ 
Evaluating the numerator and introducing the tensor integrals D and C we arrive at 



2^2 



hrh'^{ki-k2r 



+^i(-2£^ k]^ - 8pt £^ + 2p^ e^) +^iYh 

-^2{-2e^i k\ - 8p^ e\ + 2f^ e\) - ^17''^ 2 8^^ 

+^1 ^e'i e\ + 7''(-16£?^ ^2^2 + 16£^ Pi£i + 2(A;i - p^Y e^e^ 

2{k, - k2r + ^irhiMl - Ak^k2) 
+^i(-£^(3M| + 3M^ + 2piA;i) - 4A;^ £2^)1 + 4p^ £2^)1) 
-^2(£i(3M| + 3M^ + 2p2A;2) + ^k^ e^k2 - 4p^ £1^2) 

+ |^l(2(/Ci -pi)^ £l£2 - £2P2 + 8£^ £iPi) 



(6.4) 
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+^ie2ki{Ap2k2 - 2M^ + 2M|) - ^2 £1^2(4^1^1 - 2M^ + 2M|) 



[A7^^2 - ^2 - ^1 34 + 7'^£i£2 

+Co ^1 2£2^"i - ^2 2£iA;2 + ^1 362^2 - ^2 3eiPi - ^1 4eie2 ^uj_u{p2). 

The three-point integrals arise from terms in the numerator by writing q^={q^ — 
M^)+M'^ and canceUing the first denominator factor. After that the shift q q — ki +pi 
was performed in the three-point integrals (this shift conserves the manifest CP symmetry 
of the diagram) . The arguments of the C and D functions are as follows 

D = D(ki, ki - pi, -k2, Mz, Mw, 0, Mw), 

C ^Cipi,-p2,0,Mw,Mw) . 
Inserting the tensor integral decomposition eqs. (4.7, 4.8) yields the final expression 



(6.5) 



2„2 



5M 



a c 



w 



{M^ [20 Doo + 2(4M^ - 5)^33 + 2{M^ + t)D22 



+ {12M^ + At- 2s)D23 + 2{AMl - s)Drs + {16Ml - 6s + 2M|)L>3 
+ {2t -2s + 6M^ + Ml)D2 + (4M2, - 2s + MI)Dq 



+ 



-4Co - 16 D003 - 8 D002 + 10 L'oo + 2tD22 + 2{M^ + t){D^^ + D13) 



+ 



+ 



+ 2(M2, + 3t)D23 + 2(M2, - 2i + M|)D3 + (M| - t)D2 + (t - 3M|)Do 

■^3:1 + A^3",2] [-3C2 + 2Co - 8L>oo3 - 8L>oo + (4s - IIM^, + ht)D^^ 
- 3{M^ + t)D33 - 2(2M2, + t)D23 + {t - AM^ - SM^D^ + 2(M| - t)^^ 

Mil + Ml2\ [4C2 + 3Co - 8D002 - 26^00 + 2(s - AMl)iDss + D,s) 



-2{t + 2M^)D22 + (4s -2t- 18M^)n 



23 



+ {4s - 8M^ - 2M|)L>3 + {t- 4M2, - 371^1)^2 



16^113 + 8^123 - 8D13 



M; 



8D222 + 16£'223 + 24^22 + 32^23 + 16£'2 
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+ 



•^7,1 + -^7,2] [8^332 + 8L>223 + 8L>123 + 8L>23 + IBi^ia] } , 

where we introduced the standard matrix elements (5.14), the Mandelstam variables 

s^{pi+P2)\ t^{p,-hf, (6.7) 

and put 





^1 ~ ^2 




= P2 = 0. (6.8) 


Furthermore we made use 


of the relat 


ions which follow from the symmetry of the diagram 


under the exchange e"*" ^ 


e-,W+ ^ 


W~ (CP invariance) 






Dii = L>33, 


D12 — -D23, 




= -Do03) 


-D112 = -D332, 


i^l22 = i^223, (6.9) 


Dm 


= -D333, 


-D113 = -D133, 


Ci = C2. 



These reduce the number of independent invariant integrals considerably. Note that not 
all of the 40 standard matrix elements of (5.14) appear in (6.6). This is due to the 

neglection of fcrmion masses and CP invariance of the box diagram. 

This example shows that the reduction method is straightforward and universally 
applicable to one-loop Feynman diagrams, since they all have a similar structure. 

6.2 Generic Feynman diagrams 

The huge number of algebraic calculations makes the evaluation of each Feynman 
diagram very lengthy and tedious. Furthermore there are a large number of diagrams 
contributing to each process. Fortunately many of the diagrams resemble each other 
in their algebraic structure and can be considered as special cases of generic diagrams. 
These are the Feynman diagrams of a theory with only one generic scalar, fcrmion, vector 
boson and Faddeev-Popov ghost each and arbitrary renormalizable couplings between 
those fields (for more details see [34] ) . It suffices to do the algebraic calculations for these 
generic diagrams only. All actual diagrams are obtained from those by substituting the 
actual fields together with their coupling constants and masses. This saves a lot of work 
especially if there are many fields in the theory. 

Clearly the generic diagrams can be calculated with the methods described above. 
The efficiency of generic diagrams is illustrated in the next section using the decay of the 
VF-boson into massless fermions as example. 

6.3 The decay W fifj for massless fermions 

We will now apply the methods described above by calculating the one-loop amplitude 
for the decay of the 1^-boson into massless fermions 

w+{k)^Mp,)f;{P2). (6.10) 

In lowest order there is only one Feynman diagram (Fig. 6. 2) leading to the amplitude 

Mo = -^u{p,)^{k)u;MP2) = -^MT. (6.11) 
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Figure 6.2: Born diagram toW^ fif',. 







Figure 6.3: One-loop diagrams and corresponding counterterm diagram to — > /j/'. 



Neglecting the fermion masses the amplitude (6.11) leads to the following lowest order 
decay width 

aMw ,2 

i — — _ n I Vii 



6 2.2^.''^^' 



(6.12) 



At one-loop order there are six loop diagrams and one counterterm diagram (Fig. 6.3; the 
counterterm is indicated by a cross) . 

The first two loop diagrams correspond to one generic diagram and the other four to 
another generic one. We first calculate the two generic diagrams. The expression for the 
first reads 



^^'^'^ J (2^ 



1 



(27r)^(g2-M2)(g + pi)2(g-p2)2 

(6.13) 
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where denote the generic left- and right-handed fermion-fermion- vector couphngs. Sim- 
phfication and decomposition into tensor integrals yields 



{27r)^ {q^-M^)(q+p,y{q-p2)^ 



{9l 92 93^- + 9l 92 93 ^^+)v{p2)- 

Insertion of the Lorentz decomposition and further simplification gives 

SMi = -Ti^{9i 93 92Mi + gtgt92Mt) 

[(2 - DfCoo - 2k\C^2 + Ci + C2 + Co)]. 

Finally the reduction of the invariant integrals and the use of (4.55) leads to 

SMi = -T^(^r^J^2"-A^r + 9i 93 92 Mt) 

[-2k'CoiO,e,0,M,0,0){l + Mly 

-Bo{k\ 0, 0)(3 + 2^) + 2i?o(0, M, 0)(2 + ^) - 2] 
h{9i 93 92 M^ + 9t9t9tMt)Va{Q, k\ 0, M, 0, 0), 



(6.14) 



(6.15) 



(6.16) 



167r 



where we introduced the generic vertex function Va which is defined in the general case 
in App. C. 

Similarly we obtain for the second generic diagram 



J (27rW 



d^q u{pi)Y{9i ^- + 9i ^+){-i)lp{92 ^- + 92 ^+)v{P2) 



93 



(27r)^ q'[iq+Pir - M^Q - P^)' - Ml] 

fi-p^bi + 2P2 - q)u - gixv{2pi +P2 + q)p + gup{2q + pi- P2)p\£'^ 



= ^93{9i92M^ + 9t9tMt) [^{D - l)Coo - 2k\Ci2 + Ci + C2) 

= ii^93{9i92Ml +gtgtMt) 

[2{Ml + M| + ^)(:7o(0, k\ 0, 0, Ml, M^) - (1 + M+M)5o(P, M^, M^) 
+(2 + ^)5o(0, 0, Ml) + (2 + ^)5o(0, 0, M2)" 

= ^93{9i92M^ + 9t9tMt)V^{Q, P, 0, 0, M^, M^). 

(6.17) 

The general definition of can again be found in App. C. 

Inserting the actual couplings and masses of the six actual diagrams into the generic 
diagrams and adding the counterterm diagram, which can be easily obtained from the 
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Feynman rules or the Born diagram, we find for the virtual one-loop corrections to the 
invariant amplitude for W fifj 

{QfQf'VaiO,M^,0,X, 0,0) 
+gJgJ,Va{0,M^,0,Mz,0,0) 

+QfV^{0,M^,0,0,X,Mw) 

-Qf>V^{0,M^,0,Q,Mw,X) 
+^gjV^{0,M^,0,0,Mz,Mw) 

--^g-V^{0, M^, 0, 0, Mw, Mz) 

^SZli"-^ + + \5Zw + 

The left- and right-handed couplings gfj of the fermions to the Z-boson are defined in 
(A. 14). Note that only one out of the four standard matrix elements (5.11) is contributing 
there and that we need no counterterm to the quark mixing matrix for massless fermions. 
The counterterms are expressed in terms of the self energies in Sect. 3.3. These have to 
be calculated to one-loop order to determine 5M. completely. 

5A4. contains infrared divergencies. These are regularized with a photon mass A. They 
drop out in the decay width if the contribution from the decay W —>■ fif'j'y is added. This 
will be discussed in more detail in Chap. 7. 

The example above was rather simple. If we keep the fermion masses finite or consider 
processes with more external particles the number and complexity of Feynman diagrams 
raises considerably. 



(6.18) 



Ssw "I 
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6.4 Computeralgebraic calculation of one-loop diagrams 

The procedure of generation and algebraic reduction of Feynman diagrams as described 
above is algorithmic and can be implemented in symbolic computation systems. There 
are several attempts to create such systems for high energy physics calculations [35]. In 
addition there exist special packages written in general purpose languages [9, 10, 11, 12]. 
In particular the MATHEMATICA packages FEYN ARTS [11] and FEYN CALC [12] 
have been developed for the automatic calculation of one-loop diagrams following the 
approach outlined in this paper. 

FEYN ARTS generates all graphs to a given process in a specified model together with 
their combinatorial factors (weights). It yields both analytical expressions and drawings 
of the graphs. There is a version under development which uses the concept of generic 
diagrams. It creates all relevant generic graphs together with a list of all possible substi- 
tutions yielding the actual graphs. 

FEYN CALC performs the algebraic evaluation of Feynman diagrams. It starts from 
the output of FEYN ARTS and uses exactly the reduction algorithm described above. It 
can deal with generic diagrams. The FEYN CALC output can easily be translated into 
FORTRAN code. 
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7 Soft photon bremsstrahlung 



As mentioned in the last chapter the virtual one-loop corrections to the decay matrix 
element W — > /j/j are infrared divergent. These divergencies originate from photonic 
corrections and show up in any process with charged external particles. 

However, these processes are not of direct physical relevance since they cannot be 
distinguished experimentally from those involving additional soft external photons. Since 
the photons are massless their energies can be arbitrarily small and thus less than the 
resolution of any detector. Therefore in observable processes in addition to the basic 
process those with arbitrary numbers of soft photons are included. 

For these observable processes one obtains theoretically satisfactory results. Adding 
incoherently the cross sections of all the different processes with arbitrary numbers of 
photons, all infrared divergencies cancel [36]. This cancellation takes place between the 
virtual photonic corrections and the real bremsstrahlung corrections order by order in 
perturbation theory. To one-loop order one only needs to consider single photon radiation. 
For the cancellations only the soft photons, i.e. photons with energy ko < AE, are relevant, 
where AE is a cutoff parameter, which should be small compared to all relevant energy 
scales. Photons with energies fco > are called hard. They can also yield sizable 
contributions especially arising from photon emission collinear to the external charged 
particles. 

In Sect. 7.1 we introduce the soft photon approximation and show that in this ap- 
proximation the bremsstrahlung diagrams are proportional to the Born diagrams. The 
corresponding soft photon cross section for arbitrary Born diagrams is given in Sect. 7.2. 

7.1 Soft photon approximation 

Attaching soft photons to a charged external particle line of an arbitrary Feynman 
diagram yields diagrams which become singular for vanishing momentum of the soft pho- 
ton. This divergence arises from the propagator of the charged particle generated by the 
inclusion of the radiated photon line. In the soft photon approximation the momenta 
of the radiated photons are neglected everywhere but in this singular propagator. This 
approximation is valid if the matrix element of the basic process does not change much 
if a photon with energy AE is emitted, i.e. the basic matrix element is a slowly varying 
function of the photon energy for ko < AE. This is not the case if the basic process 
contains a narrow resonance as e.g. in e+e~ — > i^'^n' ■ Then one must either choose AE 
small compared to the width of the resonance or take into account the strong variation 
exactly [37, 38]. 

We now extract the soft photon matrix elements for external fermions, scalars and vec- 
tor bosons. The general renormalizable couplings of these particles to the photon allowed 
by electromagnetic gauge invariance are (momenta are considered as incoming): 




F 




(7.1) 



F 
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= -ieQs(p-p')n, 



(7.2) 




v:,p' 



Vn,P 



+ieQvgupip - p')^ - ieKv[kpg^^ - Kg^p]. (7.3) 



Quartic boson couplings do not give rise to IR-singularities and are thus not relevant in 
the soft photon approximation. The terms involving the charges Q are obtained directly 
from the covariant derivative with respect to QED. The term proportional to Ky, which 
contributes only to the magnetic moment, is gauge invariant by itself. Further terms 
present in the '-/WW coupling in the SM do not contribute for physical vector bosons. 
Since we will use the unitary gauge in this section they drop out. In renormalizable gauges 
their contributions are cancelled by those of the corresponding unphysical Higgs bosons. 

Consider first radiation from a fermion line. Let the basic matrix element without soft 
photons be 



Mo = A{p)u{p) 




(7.4) 



where u(p) is the fermion spinor with momentum p, p^ — w? and A{p) the remaining part 
of the matrix element. Inserting one photon (polarization vector £, momentum k) into 
the fermion line yields 

£, k 



Ml 



F,p 




(7.5) 



Anticommuting ji — with ^ and using the Dirac equation this can be written as 



Ml = 



-2pk 



A{p - k)[2ep - i£^(7^,r]ii(p), 



(7.6) 



where cr^j, = |[7^,7jy]. The denominator ^ contains the IR-singularity. Neglecting all 
terms proportional k in the numerator we obtain the soft photon approximation 



Mi,s = -eQF-j^A{p)u{p) = -cQf-j^Mq. 



(7.7) 



Note that the contributions of the magnetic moment term, the second term in the square 
bracket in (7.6), are neglected in the soft photon approximation and that the soft photon 



48 



matrix element is proportional to the Born matrix element. For an outgoing fermion 
{u{p)) one finds in the same way 

M^,s = eQp^Mo. (7.8) 

This is equivalent to (7.7) apart from a minus sign originating from the different charge 
flow. 

For an external vector line with polarization vector ev{p) the basic matrix element is 
Mo = A,{p)eUp) = . (7.9) 



The corresponding soft photon contribution is obtained from 

ie[Qvgupi2p - k)^ - nv{kpg^u - k^gi^p)], 

using 

evp = (7.11) 

as 

A^i,. = -eQv^Mo. (7.12) 

pk 

It is proportional to the Born matrix element and independent of the contribution of the 
magnetic moment Ky of the boson V. Again an outgoing vector yields an extra minus 
sign. 

The soft photon matrix element for an external scalar line is derived analogously. 

Radiation from internal charged lines or quartic vertices does not lead to IR- 
singularities and is neglected in the soft photon approximation. 

Summarizing, the 0{a) soft photon matrix element corresponding to an arbitrary 
matrix element Mo can be written as 

MM = -eA^oE(±^30^, (7.13) 

where Pi, Qi are the momentum and the charge of the i-th external particle and k is the 
outgoing photon momentum. The ± sign refers to charges flowing into or out of the 
diagram, respectively. The soft photon matrix element is always proportional to the Born 
matrix element. The proportionality factor depends only on the charges and momenta of 
the external particles. 

7.2 Soft photon cross section 

The soft photon cross section is obtained by squaring the soft photon matrix element 
(7.13), summing over the photon polarizations and integrating over the photon phase 
space with |k| < 

f da\ / da\ r d^k ^ ±piPjQiQ 



- 



dVt ) ^ \dVt y Q {2'kY J\k\</^E 2ujk ~^ Pikpjk 



ZlTu ^ (7-14) 
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where 

= Vk2 + A2 (7.15) 

and ± refers to the relative sign of the i-th and j'-th term in (7.13). As in the virtual 
corrections the IR-singularities are regularized by the photon mass A. Note that these 
integrals are not Lorentz-invariant due to the integration region. The basic integrals 



'|k|<AE 2ukPikpjk 

have been worked out e.g. by 't Hooft and Veltman [28]. 
The result is 



+ 



- log + Ll2 1 — +Ll2 1 



4 Wo + u V V J V V 




with 



{apif - p2 







(7.18) 



and a defined through 



ay,-2ap,pj+p',^0, ^^^_^>o. (7.19) 

Pjo 



For Pi = pj this simplifies to 



lii = 27r <^ log ——- + — log —— \ , (7.20) 



and for = — p ■ = p 



Po + 


IpI 


Po - 


IpI 


qo + 


IpI 


qo - 


IpI 



^^J' = 27ry ^ ^ <^ - log — log -— Li2 —— - 7 log — 7.21 

(Po + go)|p| 12 Po-|p| A2 Vpo + |p|/ 4 - 

^1, go + IpI . 4AE2 / 2^\ 1 . 

+ log ^ log — - 7 log 

2 go - IpI A^ Vgo + |p|/ 4 

Inserting the results for I^j into (7.14) yields the soft photon cross section. Adding it to 
the one- loop corrected cross section for the corresponding basic process the IR-divergencies 
cancel and the limit A — > can be taken. 

Although the inclusion of the real soft photon emission is sufficient to obtain IR-finitc 
results, it is often not adequate for real experiments, because realistic detectors do not 
provide a sufficiently small resolution AE/E necessary for the validity of the soft photon 
approximation. Therefore also hard photons (with ko > AE) are important. Their 
contribution is UV-and IR-finite and can be treated separately. One merely has to make 
sure that the soft and hard part are properly adapted to each other. 

Hard photon corrections are treated with methods different from the ones presented 
in this work. Their contribution depends sensitively on the experimental setup. They are 
usually incorporated by Monte Carlo simulations [15]. 
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8 Input parameters and leading higher order contributions 

In order to complete all ingredients necessary for the calculation of radiative correc- 
tions we have to specify the input parameters. This is done in Sect. 8.1. The leading 
higher order corrections which become important for precision experiments are discussed 
in Sect. 8.2. 

8.1 Input peirameters 

In its original symmetric version the SM depends on the parameters (2.21), which are 
essentially the couplings allowed by the SU{2)w x U{1)y symmetry. These were replaced 

by the physical parameters (2.22), i.e. the particle masses, the electromagnetic coupling 
constant, and the quark mixing matrix. In the on-shell renormalization scheme the renor- 
malized parameters are equal to these physical parameters in all orders of perturbation 
theory. 

The numerical values of the physical parameters must be fixed through experimental 
input. However, this input may not necessarily consist of direct measurements of the 
rcnormalized parameters; it may be obtained from any suitable set of experimental results. 
In practice one uses those experiments which have the highest experimental accuracy and 
theoretical rehabihty. This criterion is certainly fulfilled for the following set of parameters 
whose numerical values are taken from [39]: 

• the fine structure constant 

a = 1/137.0359895(61) 

corresponding to the classical electron charge e — y/Ana, 

• the masses of the charged leptons 

me = 0.51099906(15) MeV, = 105.658387(34) MeV, 

rrir = 1784.1 MeV, 

• the mass of the Z-boson [5] 

Mz = 91.177(21) GeV, 

• and the Fermi constant 

Gf = 1.16637(2)10"^ GeV-^ 
which is directly related to the muon lifetime. 

We do not use the 1^-mass as input parameter because it is experimentally not known 
with comparable accuracy. 

Besides the above listed well known parameters the still unknown masses of the top 
quark and the Higgs scalar are kept as free parameters. If the minimal SM is correct, 
the present experimental data restrict the top quark mass to the region 80 GeV < rrit < 
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200 GeV [5, 8]. For the Higgs mass we use 40GeV < Mh < ITeV, where the lower bound 
is experimental [5] and the upper bound is favored by theoretical consistency arguments. 
If not stated otherwise we will use the values rrit = 140 GeV and Mh = 100 GeV. 

The quark mixing matrix elements Vij are directly taken from experiment. We use the 
parametrization of Harari and Leurer [40] as advocated by the Particle Data Group and 
choose the following numerical values for the parameters in agreement with [39] 

si2 = 0.220, S23 = 0.046, S13 = 0.007 (8.1) 

and 5 — ioT simplicity. This yields approximately the following numbers for the quark 
mixing matrix elements: 

Vud = 0.975, Vus = 0.220, Kb = 0.007, 

= -0.220, = 0.974, Kb = 0.046, (8.2) 

Vtd = 0.003, Vts = -0.046, Vtb = 0.999. 

It remains to discuss the masses rriq of the light quarks (g = d, u, s, c, h). In the 
electrowcak Lagrangian the quarks arc treated as free particles with appropriate masses. 
This is not correct due to the presence of the strong interaction. Therefore the quark 
masses can at best be considered as somewhat effective parameters. Fortunately in typical 
high energy experiments (s ^ m^) theoretical predictions depend on the quark masses 
only through universal quantities such as the hadronic vacuum polarization or the quark 
structure functions. These can be directly determined from experiment. Nonuniversal 
contributions are suppressed as m^/s and thus negligible for sufficiently high energies. 

For processes without external quarks only the hadronic contribution to the vacuum 
polarization 

n^^(.) = 5^ (8.3) 

is relevant. In perturbation theory the contribution of light quarks is given by 

n^^W-3ii Z^g-iog^). (8.4) 

a,u,s,c,o \ Q y 

The large logarithmic terms contained in (8.4) constitute a dominant contribution to 
the radiative corrections. They originate from the charge renormalization in the on-shell 
scheme at zero momentum transfer (see eq. 3.32) involving 



n^^(o) = 



(8.5) 

A:2=0 



In this quantity nonperturbative strong interaction effects cannot be neglected. Since no 
reliable theoretical predictions are available one has to extract 11^^(0) from experimental 
data. Writing 

^tliiO) = Ut^M - Re Ut^,{s) + Re Ut^,{s) 
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the unrenormalized hadronic vacuum polarization n^^(s) can be evaluated perturbatively 
for s ^ and the renormalized one Ren^^(s) is given by the dispersion relation 

Retl^^Js) = —sRe H ds' ^^^^^'^ (8.7) 

with 

^ ^(e^e-^r-^hadrons)^ 

^ ' a{e+e- ^ 7* ^ 

can be taken from experiment up to some scale s, for larger s perturbative QCD 
is used. A recent analysis [41] involving data for the energy range below 40 GeV yields 
for the contribution of the 5 light quarks in the energy region 50 GeV < s < 200 GeV 

Re fi^tP (s) = -0.0288 ± 0.0009 (8.9) 



-0.002980 



log f^^r. .n9 + 0.006307 



(92GeV)2; ■ V(92GeV)2 

For energies around the Z-boson mass this can be approximated by (8.4) using the 
effective quark masses 

m„ = 0.041 GeV, = 1.5 GcV, 

(8.10 

md = 0.041 GeV, = 0.15 GeV, m,, = 4.5 GeV. ^ ^ 

These quark masses, in particular the ones for the three lightest quarks, are effective 
parameters adjusted to fit the dispersion integral and have no further significance. 

In addition to the above parameters we need the strong coupling constant as for the 
QCD corrections. Its value at the scale of the Z-boson mass is given by [5] 

a, = 0.115 ±0.010. (8.11) 

For the numerical evaluation we use in the following 

a, = 0.12. (8.12) 

The PF-mass is determined from the parameters given above through the relation 

Ar summarizes the radiative corrections to muon decay [26] apart from the QED correc- 
tions which coincide with those of the Fermi model. It depends on all parameters of the 
SM and is given by 

M + ^ (q I ^-^ 

'sw M| AttsI, \ 2s^ 



+2—^7^ + jz::r ^ + ^ttt^ i^s 4 ■ (8-14) 



The relation (8.13) can be improved by summing the leading higher order reducible cor- 
rections. This is done in the next section (eq. 8.23). For the set of parameters specified 
above we obtain for the 1^-boson mass 

MvK = 80.23 GeV. (8.15) 
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8.2 Leading higher order contributions 



The natural order of magnitude of one-loop radiative corrections is set by the loop 
expansion parameter - ~ 0.0023. Consequently typical one-loop corrections are of the 
order of one percent. There are, however, two important types of radiative corrections 
which are enhanced by large mass ratios. The first type is associated with light fermions 
the second with a heavy top quark. These corrections can reach several percent. Conse- 
quently the corresponding higher loop corrections may become as big as several permill 
and thus have to be taken into account in predictions for precision experiments. 

8.2.1 Leading logarithms from Ught fermion masses 

The first type of enhanced corrections originates from the renormalization of a at 
zero momentum transfer where the relevant scale is set by the fermion masses. These 
are much smaller than the relevant scales in high energy experiments. The large ratio of 
these different scales leads to large logarithms which can be summarized in the universal 
quantity 

A«(s) = -Reti^\s) = ^ E N[,Q) log M + . . . ^ (Aa)ii + . . . , (8.16) 

where A'^^ is the colour factor of the fermions and the dots indicate nonleading contri- 
butions. Renormalization group arguments can be used to show [42] that the leading 
logarithms (Aq;)^^ are correctly summed to all orders in perturbation theory by the re- 
placement 

1 + (Aa)^^ ^ ] . (8.17) 

Since not only the leading logarithms but the whole fermionic contribution to Act are 
gauge invariant we will sum the latter. The gauge dependent bosonic contribution can 
not be summed, however, because this would violate gauge invariance in higher orders. 
Thus we arrive at 

1 + Aa = 1 + (Aa) /e.m + {^ol^os 1 7-^-^ + (Aq;)6„,. (8.18) 

This corresponds to a resummation of the iterated one- loop fermionic vacuum polarization 
to all orders. 

Since the leading logarithms contained in Aa originate from the charge renormalization 

constant SZ^, they are universal, i.e. they appear everywhere where a appears in lowest 
order. They can be taken into acount by replacing the lowest order o; by a running ol{s) 
defined as 

„ = «(0)^.(.)= ; . (8.19) 

Aa(0) = due to the on-shell renormalization condition for a. Using a(s) instead of a 
effectively corresponds to renormalize a not at zero momentum transfer but at momentum 
transfer s. Then the light fermion masses can be neglected everywhere and no large 
logarithms appear. 

There are similar large logarithms associated with external fermion lines. These are 
related to collinear singularities, arising from the radiation of photons collinear to external 
particles. They can be consistently treated with the structure function method [43]. 
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8.2.2 Leading mf contributions 

The second type of important corrections is also connected to the fermionic sector. 
The top quark gives rise to corrections oc m^/M^, which become large if the top quark 
mass is large compared to the VT-boson mass. For nit = 200 GcV they reach several 
percent. These terms arise from fermion loop contributions to the boson self energies and 
from the Yukawa couplings of the physical and unphysical Higgs fields, which show up in 
vertex and fermionic self energy corrections. The latter are process dependent and are 
therefore not discussed here. The former, however, are universal, they can be traced back 
to the renormalization of sw in the on-shell scheme 

5sw _ 1^^^ (^¥{M^) ^¥{Ml)\ 



Sw '2s^w \ Ml 



lc^r a 3 1 
2X4^4XM;^---" 27 



+ ... = --fAp + ... (8.20) 



where the dots indicate again nonleading contributions. 

There is no general principle that determines the resummation of these corrections. 
The following recipe has been shown to yield the correct leading terms to 0{o?) [44] 

2 2 - -2 V°-21j 



where 



2 



(8.22) 



8V27r2 

incorporates the result [45] from two-loop irreducible diagrams. Note that a has been 
replaced by Gj? in order to obtain the correct leading 0{a^) terms. 

In particular for the relation between and Gp the correct resummation of the 
leading corrections is given by 

TT a 1 1 



Gf^ 



M^V2 SvK 1 - (Aa)/erm 1 + ^Ap 



l + Ar-{Aa)ferm + ^Ap 



(8.23) 



Note that the two types of leading corrections are summed separately. Inserting = 
1 — M^/M| this relation can be used to determine the 1^-boson mass Mw including 
leading higher order contributions. 

Neglecting the nonleading contributions and using (8.19) and (8.21) eq. (8.23) can be 
written as 

^^^^ ^ V2GfM^^. (8.24) 

Sw 

With this relation the appearance of Gi? in Ap can be easily understood. All leading 
universal corrections arise from the renormalization constants of a and sw- Consequently 
they can be absorbed by incorporating the leading finite parts of these renormalization 
constants into the effective parameters a{s) and s^r (including the leading higher order 
corrections). Thus one obtains from the lowest order result directly the corresponding 
result including the leading universal corrections. In particular (8.24) can be obtained in 
this way. Applying this recipe to Ap and using (8.24) introduces naturally Gp. 
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8.2.3 Recipes for leading universal corrections 

Prom the discussion above we infer that the universal leading higher order contribu- 
tions can be taken into account by the following replacements 



a 
















2Tra{s) 






s2 


e2 




7ra(s) 


2 fO. 

w'-w 




s2 ^2 



V2GfMI- — . (8.25) 



4„2 2 2 2 ^2 ^2 ' ^'-z^ _ A - 



Note that this does not include the nonuniversal corrections oc am^/M^ arising from 
enhanced Yukawa couphngs. These have to be evaluated for each process directly. 
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9 The width of the Vt^-boson 



In the coming years the upgrade of the LEP electron positron storage ring to 180 — 
190 GeV CM-energy will allow to study the properties of the H^-boson in detail in a 
model independent way. Besides its mass Mw also its total and partial decay widths 
are of interest. While the leptonic partial widths allow to test lepton universality the 
hadronic partial widths can serve to determine the quark mixing matrix elements [46]. 
The expected accuracy is 6Mw ^ 100 MeV and 6Tw ^ 200 MeV [47]. 

An accurate comparison between these experiments and theoretical predictions of the 
SM requires at least the inclusion of one-loop radiative corrections both for the production 
and the decay of 1^-bosons. The 1^-bosons decay dominantly into fermion-antifermion 
pairs. The total and partial widths of the W^-boson and the corresponding one- loop 
radiative corrections are discussed in this chapter. 

The electroweak and QCD corrections for decays into massless fermions (m/ <C Mw) 
have been calculated in [48, 49, 50]. The hard bremsstrahlung contribution has been 
investigated in [52, 53]. The QCD corrections for the decay into a massive top quark 
and a massless bottom quark were given in [54]. The full one-loop electroweak and 
QCD corrections together with the complete photonic and gluonic bremsstrahlung were 
evaluated for arbitrary finite fermion masses in [55]. 

Since the top quark is probably heavier than the 1^-boson and since all other quark 
masses are small compared to the 1^-boson mass, the fermion mass effects are not of 
great importance for the 1^-decay. However, the matrix element for the W-decay into 
two fermions is directly related via crossing to the one for the decay of the top quark into 
a VT-boson and a bottom quark. In this case the fermion masses are crucial. Since we will 
discuss top decay in Chap. 10 we give the results for 1^-decay for finite fermion masses. 



9.1 Lowest order 

The Born amplitude for the decay — > /j/j was already given in (6.11). For finite 
fermion masses the following result is obtained for the corresponding partial decay width 

To '{Mw,mf,i,mf,j)^Nlj——^\Vij\ —G^, (9.1) 

where k/M^t originates from phase space and 



Gr = J:MT'MT = 2M^ - m% - m?,, - ^^^'^ ^^f''^' (9.2) 

pol 



from the polarization sum of the matrix element squared. The Kallen function k was 
defined in (4.28). The colour factor Nq is given by 

AT/ = / 3 for quarks, 3^ 
I 1 for leptons. 

For leptonic decays we have Vij = Sij. The total width is obtained as a sum over the 
partial fermionic decay widths with ruf^i + mfj < Mw 

= ^rr""' + E^r^'^- (9.4) 

ij i 
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We can write down another possible tree level representation for the partial decay width 
by eliminating a/s^ in favour of 



T7 



127rV2 



w 



(9.5) 



9.2 Electroweak virtual corrections 



In our formulation of the on-shell renormalization scheme the invariant matrix element 
to one-loop order has the following form 

M"^^^^'^^^^ {v,,M-Al + SZ,-^-^]+M-,5V,, 
y Zsw 

k 

+ M:5F:{Mw, mf,i, mf,j)}. (9.6) 

a=l a=± 

The standard matrix elements were defined in (5.11). The functions 5F^ summarize 
the loop corrections to the 1^/j/j-vertex. There are no explicit self energy corrections 
from the external lines. These are all absorbed into the field renormalization constants 
6Zw, SZ^'^ and 6Z-f These and the parameter renormalization constants were given 
in terms of self energies in Sect. 3.3. The exphcit forms of the self energies can be found 
in App. B. 

Fig. 9.1 shows the Feynman diagrams contributing to the vertex corrections for massive 
fermions. They yield the following vertex form factors 

_ 0( 

5F^ {Mw,mi,m2) = 

47r 

I QfQf'Va{ml, M^, ml, X, mi, 7712) + gJgJ,Va{m\, M^, ml, Mz, mi, 7712) 



+ 



E 

<7=± 



QfV;;{mi, M^,mimi,X,Mw) - QrVbi^i, M^,mi,m2, X, Mw) 



+ ^g-V^imj, M^, ml, mi, Mz, Mw) - —g'},V^{ml, M^, m\, m^, Mz, Mw) 

Sw Sw ^ 

+ j^Vciml, M^, ml, Mh, mi, 77^2) 

^Sw 



+ 



+ 



^Sw 



^Sw 



Vdimj, M^, ml, mi, Mh, Mw) + Vdiml, M^, ml, m^, Mh, Mw) 
Ve{ml, M^, ml, mi, Mh, Mw) + Ve{ml, M^, ml, m^, Mh, Mw) 



(9.7) 
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+ Ve(mf , M^, mi mi, M^, Mw) + Ve(m^, M^, mf , ms, M^, M^.) 



o-=± 



^^JV;(m2, M^,, ml mi, M^, M^^) + '-^g-^iml Ml, ml m^, Mz, Mw) 



a 



dF^{Mw,mi,m2) = — mim2 x 
j E fo/Q/' Wa"K, ^^2, A, mi, ms) + M^, m^ M^, mi, ma) 



+ QfW^imi M^, m^, mi. A, Mh.) - Qf'W^{mi M^, mf, ma. A, M^^) 

+ — ^JW^m?, M^,, ml mi, M^, M^) - W.^K, M^,, m?, ma, M^, M^) 

Sw Sw 



+ 



4s2 



^ >V,"(m2,M^,m2,MH,mi,m2) - >V-(m2,M2„m2,Mz,mi,m2) 

■(T=± 



+ 



+ 



0„2 



0„2 



Wd{ml M^, m^, mi, Mj^, Mw) + Wd{ml M^, m?, m2, M^, M^y) 
WeK, M^, ml mi, M^^, M^.) + >Ve(m^, M^, m^, m2, M^, Mw) 



(9.8) 



- WeK, M^, ml mi, Mz, M^y) - We{ml M^, m^, m2, M^, Mw) 



+ 



E 

<T=± 



QfWj{ml M^, m^, mi. A, M^/) - Qf'Wj{mi M'^, m{, m2. A, Mw) 



C('™2 71^2 _2 



^^J>VJ(m2, M^, m^, mi, M^, Mw) + ^^J, WJlm^, M^, m?, m2, M^, M^.) 



(M^,mi,m2) = — miX 

47r 

S E QfQf'^ai'nT'l, Mw, ml A, mi, m2) + gfgJ,X^{ml Mw, ml Mz, mi, m2) 
U=±'- 



+ QfX^iml M^, m^, mi. A, Mw) + —gJX^iml M^, ml mi, Mz, M, 

Sw 



Qf'X,, {ml, M^, m2, m2, Mw, A) - — fi-^/A"^ (m^, M^, m2, m2, M^y, M^) 

Sw 
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+ 



4q2 



J2 {^l-> ^w-> ^l-> ^H, mi, 7712) - {mj, M^, mj, Mz, mi, ma) 

■a=± 



+ Tnr^d{m\, M^, m\, mi, Mn, Mw) 



(9.9) 



0„2 



^ [cr(A';(m?, M^, m^, mi, M^, M^) - X^{ml M^, m^, ma, Mh, Mw 



- X^{mi M^, mi mi, M^, M^y) - X^{mi M^, m^, ma, M^, Mw) 



1 ml 



A;"(mf , M^, m^, mi, Mj^, M^.) - A;"(m^, M^, mf, ma, Mh, Mw) 



- X^{mi M^, mi mi, M^, M^.) - X^{mi M^, m^, ma, M^, Mw 



+ QfXf{mi, Myy, ma, mi. A, Mw) + (^/"^/(ma, M^, mi, ma. A, M, 



Cw 



g}Xf{ml M^, ml mi, Mz, Mw) + c//'A'/(m|, M^, m\, ma, Mz, M^^) 



To obtain these results we had to evaluate six generic diagrams labeled by a, b, c, d, e and 
f in Fig. 9.1. The corresponding invariant functions V, W and X are listed in App. C. The 
formfactor can be obtained from 5Fa~ by the substitutions mi ma, Qf ~Qf'-i 
9} ^ -9f'- 

Squaring the matrix element (9.6), summing over the polarizations of the external 
particles and multiplying the phase space factor yields the one-loop corrected width 



"c 



12 2slr 



{\Vij\^G^[l + 2SZ, - 2^ + 5Zw] 

Sw 



+ mj\'i: E G:5F:{Mw,mf,„mf,j) 



=lo-=± 



(9.10) 



where 



= r™^'(i + CJ, 



(9.11) 
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^ fi 






Figure 9.2: Bremsstrahlung Feynman diagrams for W fifh- 



^2 =E-^r^-^2 

pol 



and we have inserted SVij from (3.23). 
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9.3 Photon bremsstrahlung 

Like any one-loop amplitude with external charged particles (9.6) and consequently 
(9.10) are IR-divergent due to virtual photonic corrections. These singularities are com- 
pensated by the the real bremsstrahlung corrections, i.e. the three-body decay 

W+{k) ^ Mpi)f;{p2h{q). (9.12) 
The corresponding matrix element as given by the Feynman diagrams (Fig. 9.2) is 



-u{Pi) 



-Qf 



(Qf - Qf 



2pir) i 



Qr 



2P2? 



2p2^ i + iii) 



+ 



-2kq 



{qrj - 2kri) ^+2eri ^- 2qe fj \uj_v{p2), 



(9.13) 



where rj denotes the polarization vector of the photon. Performing the polarization sum 
over the square of the amplitude gives 



pol 



a 



QfQ 



2^,^(647r^)|y,,| I ^2p,g) (2p2g) 



(M^ - m% - m}j)G^ 
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f 



(Qf-Qf'T ' 

{-2kqy 

QAQf-Qf) 

—2kq 2piq 

QAQf-Qf) 

—2kq 2p2q 



2 2 

{m% + 2p,q)Gi + (1 + !!^:^i±^)(2pig) {-2kq) + {2p,qf 

2 2 

{m},j + 2p,q)GT + (1 + '"^'2M^^''' )(2P2g) (-2/^?) + (2p2?)^ 



2 2 

(M^ - 2^g)G7 + !!^^|±^(-2^g)2 - 2(2pig)(2p2g) 



(M^ + ml, - mlj)G^ - 2{2p^q){2p2q) 



(M^ - m% + mj,^^.)^^ - 2{2p,q){2p2q) 



(9.14) 



Prom this the complete bremsstrahlung cross section (including soft and hard photons) is 
obtained by integrating over the phase space of the photon and the two fermions as 

1 A^^ f d^q d^pi d^p2 



(27r)5 2Mw 



J 2qo 2pio 2p20 o ^ 



with 



(9.15) 



a 



77 J K{M^,ml,,m}^j) 



-QfQf 



+Q} 



[ml. hi + h) + (1 + 



2M^ 



+Ql 



{m},jl22 + h) + (1 + 



2M^ 



HQf - Qf'f [(M^/oo + lo) + ""^M^'^ i^ - 



+Qf{Qf - Qf) (M^ + m% - m),^^)h, - 2^' 



-QAQf-Qf) 



{Mlr-m% + m),^^)h2-2-^ 



(9.16) 



The bremsstrahlung phase space integrals I,,', = r,',',{M\Y,mf i,mfrj) are given in App. D. 
The IR-singularities contained in the I^i are again regularized by a photon mass A. 



63 



9.4 QCD corrections 



Like the electroweak corrections also the QCD corrections consist of virtual and real 
contributions 

5«^^(M;^, m^,, mr,) = 5Tf + 5?^^, (9.17) 

which are individually IR-divergent. They are obtained from the electroweak results of 
the previous sections by keeping only the terms containing Qj, Q^fi or QfQfi, setting Qf — 
Qfr = 1, replacing a by the strong coupling constant as and multiplying an overall colour 
factor Cp = ^. In particular the virtual QCD-corrections arise only from the diagram of 
Fig. 9.1a with the photon replaced by a gluon and the corresponding corrections to the 
fermion wave function renormalization constants. Since this allows many simplifications 
we give the explicit results 



xQCD 

virt 



+2(M^ - - m}^j){Co + Ci + C2) - 2mlfii - "Im^ f^ (9.18) 



-2 log(^^^^) - 20m,,m,,, \C, + C, 



Q + 

C11 + C12 + C1 +4-^m//j 



C22 + C*12 + C2 



The arguments of the three-point functions are C — C{m?j^^, "^/'j' ^^f',i)- To 

the gluonic bremsstrahlung 5^^^ only the first three terms of (9.16) contribute. For zero 
fermion masses the total QCD-correction reduces to 

(9.19) 



5«^^(Mvi.,0,0) = 



TT 



9.5 Results and Discussion 



For numerical evaluation of the previous results we use the parameters listed in 
Sect. 8.1 including the values for the quark mixing matrix as given by (8.1). The W- 
mass is determined from the relation (8.23). 

In the on-shell scheme the lowest order width is parametrized by a and the particle 
masses (9.1). In this scheme large electroweak corrections arise due to fermion loop 
contributions to the renormalization of a and s\a/. We can improve the results in this 
scheme by resumming the corresponding one-loop contributions to all orders as discussed 
in Sect. 8.2. Thus we obtain for the corrected width 



r^Fo 

= Fo 1 + 5 , 



l+(5i-(AQ;)/erm+^Ap" 



1 - (AQ;)/erm 1 + ^Ap (9.20) 
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Figure 9.3: Electroweak radiative corrections 5^'^ and S'^'^ and QCD-corrections S^'"^ to 
the total ly-width versus the top mass. 



where 



_|_ Aew 
virt ' b 



^QCD 



(9.21) 



is the proper one-loop correction without resummation. As in any charged current process 

we can avoid the large corrections by parametrizing the lowest order decay width with 
Gp and Mw instead of a and (9.5). Using (8.23) we find the relation between the 
decay width in both parametrizations 



1 



1 - {Aa)ferm 1 + ^Ap 



1 + Ar- (Aa)/er-m + T^P 



(9.22) 



1 + Si-Ar\ =fo(l + 5). 

The large fermionic contributions contained in 5i are exactly cancelled by equal contri- 
butions in Ar and consequently the remaining corrections S are small. 

The relative corrections to the total VT-boson decay width are shown in Fig. 9.3. 
They are large and strongly mj-dependent in the on-shell scheme. This behaviour arises 
from the fermionic contributions to the renormalization of the weak mixing angle in (9.6) 
{Ss^^/s'^r) which contain terms oc am^/M^. In contrast to this in the parametrization 
with Gp the corrections depend only weakly on rrit and remain below 0.6%. The QCD 
corrections are practically constant and equal to 2Q;s/(37r), their value for zero fermion 



masses. 
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/;/, 


Mn 




rir 
^ 1) 




fir 


80.0 


79.87 


1.8778 


2.0056 


1.9866 


1.9910 


100.0 


79.99 


1.8996 


2.0147 


2.0002 


2.0063 


120.0 


80.10 


1.9211 


2.0235 


2.0110 


2.0158 


140.0 


80.23 


1.9447 


2.0330 


2.0221 


2.0256 


160.0 


80.37 


1.9713 


2.0435 


2.0340 


2.0363 


180.0 


80.52 


2.0018 


2.0553 


2.0469 


2.0481 


200.0 


80.69 


2.0368 


2.0681 


2.0608 


2.0612 



Table 9.1: W-decay width parametrized by a and parametrized by in lowest 
order and including electroweak corrections. All values are given in GeV. 

Tab. 9.1 shows the lowest order width Fq and the width including electroweak cor- 
rections Fg^ in both paramctrizations for various values of the top mass. The W^-mass 
obtained from (8.23) is also listed there. While the results for the lowest order width in the 
two paramctrizations differ by several percent, the deviation of the first order expressions 
is always less than 0.4%. 

The analytical results were presented above for finite external fermion masses and 
with correct rcnormalization of the quark mixing matrix. However, since the top quark 
is presumably heavier than the l^-boson [57] all relevant actual fermion masses are small 
compared to the VF-boson mass. Therefore in addition to the completely corrected numer- 
ical results r(Mw,Tnf^i,mf'j) for finite fermion masses we also give those for vanishing 
fermion masses F(Mvi/, 0, 0) in Tab. 9.2. The analytical results for vanishing fermion 
masses were listed in Sect. 6.3. Finally we include an improved Born approximation con- 
sisting of the Born widths with zero fermion masses parametrized hj Gp and multiplied 
by the QCD correction factor for zero quark masses 



V 

imp 



6V27r 

GpMh 



rS' = ^|v.l%i + v)' (9.23) 

^ _ 3G^M3. / 2 a, 
'"^^ 2v^7r V ^Stt 

for the leptonic partial widths, the hadronic partial widhts and for the total width, re- 
spectively. 

The numerical values for the partial and total 1^-widths in these different approxi- 
mations are given in Tab. 9.2 assuming a top quark mass of 140 GeV and a Higgs boson 
mass of 100 GeV. The improved Born approximation (9.23) reproduces the exact results 
up to 0.4% (0.6% for the decays into a b-quark). The effects of the fermion masses are 
below 0.3%. They are suppressed by w?/s. There are no mass singularities since the 
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Born 
width 


complete 
one-loop 


m/ = 
in r 


improved 
Born 


Branching 
ratio 




0.2260 


0.2252 


0.2252 


0.2260 


0.1084 




0.2259 


0.2252 


0.2252 


0.2260 


0.1084 




0.2258 


0.2250 


0.2252 


0.2260 


0.1083 


r(W lep.) 


0.6777 


0.6754 


0.6756 


0.6778 


0.3251 


r(w ud) 


0.6450 


0.6672 


0.6672 


0.6696 


0.3211 


r{w us) X 10 


0.3281 


0.3393 


0.3393 


0.3406 


0.0163 


r(W ub) X 10^ 


0.3306 


0.3428 


0.3436 


0.3448 


0.00002 


r(W cd) X 10 


0.3281 


0.3396 


0.3396 


0.3409 


0.0163 


r(w cs) 


0.6432 


0.6656 


0.6657 


0.6682 


0.3204 


r{W cb) X 102 


0.1427 


0.1480 


0.1484 


0.1489 


0.0007 


T(W had.) 


1.3553 


1.4022 


1.4023 


1.4075 


0.6749 


r{w all) 


2.0330 


2.0776 


2.0779 


2.0853 





Table 9.2: Partial and total W-decay widths F in different approximations for 
mt = 140 GeV, Mh = 100 GeV and the corresponding W-mass Mw = 80.23 GeV. 

width is obtained by integrating over the full phase space of the final state particles [56] . 
Consequently the exact numerical values of the masses of the external fermions masses 
are irrelevant^. The branching ratios derived from (9.23) 



BR{W Iv) 
BR{W leptons) 
BR{W Uidj) 



1 



9(l + 2a,/37r)' 
1 

3(l + 2Q;,/37r)' 

\V,,\\l+as/7r) 
3(l + 2a,/37r) ' 



BR{W ^ hadrons) = + "^^Z^) (9.24) 

3(1 + 2as/3n) 

agree numerically within 0.1% with those obtained from the full one loop results. They 

depend only on and Vij. 

The dependence of the M^-width on the unknown top and Higgs masses is shown in 
Tab. 9.3 for Y^^'^ and in Tab. 9.4 for Y^"^". A variation of rrit between 80 and 200 GeV 
affects the partial widths by ~ 4%, a variation of between 50 and 1000 GeV by 

^Thc values given in Sect. 8.1 are not appropriate for the external quarks. We only use them here to 
demonstrate the numerical irrelevance of the fermion mass effects. 
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"'/ 


Mn = 50 


Mil = 100 


.1/// = 30(J 


.1/^/ = 1000 


80.0 


0.2212 


0.2209 


0.2202 


0.2193 


100.0 


0.2227 


0.2224 


0.2217 


0.2208 


120.0 


0.2239 


0.2236 


0.2229 


0.2220 


140.0 


0.2251 


0.2248 


0.2241 


0.2232 


160.0 


0.2264 


0.2261 


0.2254 


0.2245 


180.0 


0.2278 


0.2276 


0.2269 


0.2260 


200.0 


0.2291 


0.2291 


0.2281 


0.2275 



Table 9.3: Partial W-decay width F^^'^ including first order QCD and electroweak cor- 
rections for different values of the top and Higgs masses. All values are given in GeV. 



ITlf 


Mh = 50 


Mh = 100 


Mh = 300 


Mh = 1000 


80.0 


0.6539 


0.6530 


0.6509 


0.6481 


100.0 


0.6585 


0.6575 


0.6554 


0.6526 


120.0 


0.6622 


0.6612 


0.6591 


0.6563 


140.0 


0.6659 


0.6650 


0.6629 


0.6601 


160.0 


0.6700 


0.6691 


0.6670 


0.6642 


180.0 


0.6745 


0.6736 


0.6715 


0.6686 


200.0 


0.6793 


0.6784 


0.6763 


0.6735 



Table 9.4: Partial W-decay width r including first order QCD and electroweak cor- 
rections for different values of the top and Higgs masses. All values are given in GeV. 

~ 1%. This holds as well for the total width as shown in Tab. 9.5. All this is valid 
for constant a, Gp and Mz- In this case the top mass dependence is mainly due to the 
variation of Mw with rrit. Keeping instead Mw, Gp and Mz fixed the dependence on rrit 
is considerably smaller. Remember, however, that the prediction for the decay width has 
the same uncertainty in this parametrization due to the uncertainty of the experimental 
value for the VF-boson mass. 

We have compared our results for the partial leptonic width for zero fermion masses 
to those of Jegerlehner [50] and Bardin et al. [49], who both use the parametrization 
with Gp. Furthermore Jegerlehner includes two-loop QCD corrections into the boson self 
energies. If these are switched off the difference between his and our results is less than 
0.1 MeV. Performing the same comparison with Bardin et al., our values for the partial 
widths are 0.7 MeV to 0.8 MeV larger than theirs. Our results for the QCD corrections 
agree with those obtained by Alvarez et al. [54]. 
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Mh = 50 


Mh = 100 


Mh = 300 


Mh = 1000 


80.0 


2.0375 


2.0347 


2.0283 


2.0196 


100.0 


2.0517 


2.0488 


2.0423 


2.0336 


120.0 


2.0630 


2.0602 


2.0537 


2.0449 


140.0 


2.0747 


2.0719 


2.0654 


2.0566 


160.0 


2.0872 


2.0844 


2.0780 


2.0692 


180.0 


2.1009 


2.0981 


2.0917 


2.0830 


200.0 


2.1157 


2.1129 


2.1066 


2.0980 



Table 9.5: Total W-decay width T including first order QCD and electroweak corrections 
for different values of the top and Higgs masses. All values are given in GeV. 
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10 The top width 



The present lower limit from CDF data indicates that the top mass is at least 89 GeV 
[57]. Moreover, LEP data in combination with radiative correction calculations require 
rrit — 137 ± 40 GeV within the minimal standard model [5, 8] at the la level. Therefore 
the top mass hes presumably above the Wb threshold and the dominant decay of the top 
quark is the one into a VT-boson and a bottom quark {t — > Wb) and the total width of 
the top quark can be well described by the partial width F*'^^ = F(t ^ Wb). 

While the measurement of the top mass will provide a long missing input parameter, 
the measurement of its width will serve as a consistency check on the standard model. 
With the operation of LHC, SSC and/or a high energy e+e~ coUider one expects to obtain 
a sufficiently large number of tops so that both the mass and the width can be measured 
with good accuracy. 

The QCD corrections to the top decay t Wb were already evaluated in [58, 59]. 
The first order electroweak corrections have been calculated by [60, 61]. 

The electroweak corrections to this decay involve particularly interesting contributions 
of 0{am^/M^) which are potentially large for large top masses. Those terms arise not 
only from fcrmion loop contributions to the boson self energies but also from the Yukawa 
couplings of the Higgs fields, which show up in vertex and fermionic self energy corrections. 
As discussed in Sect. 8.2 contributions of the first type can be ehminated if the Born 
approximation is expressed hy Gp and Mw- Surprisingly the effects from strong Yukawa 
couplings turn out to be small, as will be demonstrated in the following. 

We will only consider the decay of free top quarks and sum over the polarizations of 
the ly-bosons. The results are obtained via crossing from the ones for W — > tb. 

10.1 Notation and lowest order decay width 

Because we want to use our results for the decay W~^ tb we consider the decay of 
an anti-top quark. The corresponding decay width is identical to the one of the top quark 
because of the CPT theorem. 

The lowest order decay of an anti-top quark 



i{pi)^W-{k)b{p2) 



(10.1) 



is described by the Feynman diagram of Fig. 10.1 yielding the amplitude 



— e 



Vtb v{pi)^u-v{p2). 



(10.2) 





Figure 10.1: Born diagram for the decay t ^ W b 
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It can be obtained from the Born amplitude (6.11) for the decay W'^ — > tb by crossing. 
This amounts to change the signs of pi and k and use u{—pi) — v(pi) and s{—k) — e{k). 
From (10.2) we get the lowest order width 



a 



8 2s^w 



,K[m 



th\ 



t ' 



Ml 



mt 



ml 



with 



m; + mi- 2M^ + 



[mt 



m: 



2\2 



(10.3) 



(10.4) 



Eq. (10.3) can directly be derived from (9.1) by substituting m^ mt and mj m^ in 
Gi coming from the matrix element squared, exchanging Mw and in the phase space 
factors, changing the spin average from 1/3 to 1/2 and supplying a minus sign originating 
from the different signs of the momenta entering the matrix element squared. This minus 
sign has been incorporated into the definition of Gi which differs from the one in Chap. 9. 
Introducing Gp instead of a the lowest order width reads 



mt,Mw,mb) 



GpM^ 



^^ K{mlM^,ml) ^_ 
\^tb\ — 



mf 



(10.5) 



10.2 Virtual corrections 



With the four Dirac matrix elements 

Mt = v(pi)^u;+v(p2), 

M2 = v{pi)uj-v{p2) t 
Mt = v{pi)uj+v{p2) e 



(10.6) 



Pi, 
Pi 



obtained from (5.11) by setting pi —pi the virtual electroweak one-loop corrections 
take the form of (9.6) with t and b instead of i and j. The corresponding decay width 
Ytwb foiiQ-^s from (9.10) using the substitutions specified after (10.4). 

The QCD corrections can be extracted from the electroweak ones in the same way as 
in Sect. 9.4. 



10.3 Bremsstrahlung 

The real photonic contributions of 0{a) to the top width arise from the radiative 
decay 

i{pi)^W-{k)b{p2h{q). (10.7) 

The corresponding amplitude squared, summed over all polarizations, can be derived from 
(9.14) by replacing the momenta k —k, pi —pi and multiplying an overall factor 
(— 1). From this we get the bremsstrahlung contribution to the top width by integrating 
over the appropriate phase space 

rrV*,M^,m,) = -1^-^ / ^^^^p^5^'\pi-P2-<l-k)W\Mb? 

(27r)5 2mt J 2qo 2ko 2p2o 2 ^ 
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The correction factor reads 



5^'"{mt,Mw,mb) = 



a 



-QtQi 



{M^ - ml - ml)h 



02 



(10.8) 



K%o + /o)-(i + 47^)/° 



ml + ml I 



(10.9) 



+Qt{Qt - Qb) 



-QbiQt - Qb) 



n 



(M^ + ml-ml)Ioi + 2^ 
(M^ - ml + ml)h2 + 2 



r::{mt,Mw,mb). 



The bremsstrahlung phase space integrals carry the arguments /■.'.■ 
They are given in App. D. 

From eq. (10.9) the gluonic bremsstrahlung corrections can be obtained by setting 
Qt = Qb = 1, replacing a by the strong coupling constant as and multiplying with the 
colour factor Cp = ^■ 



10.4 Results and discussion 



Wc again use the parameters listed in Sect. 8.1 as numerical input and calculate the 
Vl/^-mass from the relation (8.23). Unless stated otherwise, we choose for the Higgs mass 
Mh = 100 GeV. 

We perform the same summation of the leading higher order corrections as discussed 
in Sect. 9.5, eq. (9.20), and introduce the parametrization with Gp and as in (9.22). 
In this parametrization the large corrections arising from the renormalization of a and 
arc absorbed into the lowest order expression. This is not the case for large contributions 
proportional to aml/M^^ arising from vertex and fermion self energy diagrams with 
enhanced Yukawa couplings. 

In Tab. 10.1 we give the lowest order width as well as the width including electroweak 
corrections in both parametrizations for various values of the top mass together with 
the ly-mass obtained from (8.23). The results for the first order expressions of both 
parametrizations agree within 0.05%. 

According to (10.3) the width increases with the top mass approximately like mf/M^. 
The corresponding relative corrections arc shown in Fig. 10.2. The QCD corrections yield 
about —10% with only a weak dependence on the top mass. In the on-shell scheme we 
find sizable electroweak corrections which range from +7% at mt = 100 GeV to —13% 
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mt 


Mw 


-ptWb 
^ 


f^tWb 
^ 


ytWb 

ew 


YtWb 
ew 


100.0 


79.99 


0.0887 


0.0940 


0.0951 


0.0951 


120.0 


80.10 


0.3095 


0.3260 


0.3306 


0.3305 


140.0 


80.23 


0.6393 


0.6684 


0.6788 


0.6786 


160.0 


80.37 


1.0850 


1.1218 


1.1135 


1.1133 


180.0 


80.52 


1.6627 


1.7071 


1.7365 


1.7362 


200.0 


80.69 


2.3927 


2.4299 


2.4724 


2.4720 


220.0 


80.88 


3.2989 


3.3083 


3.3665 


3.3661 



Table 10.1: Top decay width F*^^ parametrized by a and F*^^ by Gi? in lowest order and 
including electroweak corrections. All values are given in GeV. 
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Figure 10.2: Electroweak radiative corrections 5^"^ and 5^^ and QCD corrections ^'^c^' iq 
the top decay width versus the top mass. 
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mt/GeV 


c\r\r\ 

200 


o r\r\ 

300 


500 


1 f\f\r\ 

1000 


leading term 


0.023 


0.036 


0.051 


-0.070 


next to 

leading tonii 


-0.043 


-0.065 


-0.108 


-0.217 



Table 10.2: The two leading terms of the expansion in eq. (10.10) for various values of 
the top quark mass. 

at vrit = 350 GeV. The large variation arises from terms oc am^/My^r in the first order 
corrections. Contrarily in the parametrization with Gp the corrections are only ^ +1% 
for rrit = 100 GeV and remain almost constant at ~ +1.7% for rrit > 160 GeV. This 
feature is independent of the Higgs boson mass. However, we expected large corrections 
oc to arise from the vertex diagrams containing large Yukawa couplings; in the similar 
case of hh the corresponding contributions are noticeable. The slow variation of 

the relative correction in this parametrization indicates that the strong Yukawa couplings 
have no sizable effect on the top width. In order to demonstrate the absence of large 
corrections the plot in Fig. 10.2 has been extended up to = 350 GeV, although this is 
well above the present upper limits on the top quark mass. 

Some understanding of this surprising feature can be obtained from the expansion of 
the relative correction factor in the parametrization with Gp and Mw-, 5^"", for large top 
quark masses {rrit 3> Mw, Mz, Mh) 

with Mi — Mw, Mz, Mh- The two leading terms are evaluated in Tab. 10.2 for a wide 
range of values of the top quark mass and Mw = 80 GeV, Mh = 100 GeV. For rrit < 
1 TeV we find the leading term of the expansion to be smaller than the next to leading 
term. Consequently the expansion is only asymptotic and the contributions oc amf/Mw 
are not dominant unless the top mass has a value of several TeV. In the physically 
acceptable range of top quark masses, the quadratic terms are numerically compensated 
by logarithmic contributions. This remains true also for large Higgs masses. The small 
corrections result from intricate cancellations between leading and nonleading terms. 

In Fig. 10.3 we show the lowest order width in both parametrizatons Fq, Fq, the 
electroweak corrections STew, the QCD corrections STqcd as well as the fully corrected 
width F as a function of the top quark mass. 

The dependence of the total width on the Higgs mass is displayed in Tab. 10.3 where 
this parameter is varied from 50 to 1000 GeV. Although the influence of Mh becomes 
stronger for large top masses it never exceeds 1%. 

We have compared the pure QCD corrections to those obtained in [58] and found 
complete agreement. Our results for the electroweak corrections agree with those of [61]. 
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Figure 10.3: Top decay width F*^^ versus the top mass in lowest order Fq, Fq including all 
corrections F and the contribution of the electroweak STew and QCD corrections STqcd- 



rrit 


Mh = 50 


Mh = 100 


Mh = 300 


Mh = 1000 


100.0 


0.0882 


0.0883 


0.0886 


0.0891 


120.0 


0.3022 


0.3021 


0.3020 


0.3022 


140.0 


0.6179 


0.6174 


0.6165 


0.6157 


160.0 


1.0385 


1.0377 


1.0356 


1.0329 


180.0 


1.5742 


1.5734 


1.5697 


1.5641 


200.0 


2.2373 


2.2376 


2.2321 


2.2221 


220.0 


3.0411 


3.0438 


3.0369 


3.0207 



Table 10.3: Top decay width r including first order QCD and electroweak corrections 
for different values of the Higgs mass. All values are given in GeV. 
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11 VF-pair production 



One of the predominant aims of LEP200 is the high precision investigation of the 
properties of the H^-boson, i.e. its mass, its total and partial widths and its couplings. 
Probably the most interesting aspect will be the study of the nonabelian gauge interaction 
which has no direct experimental evidence so far. 

The general properties of VT-pair production at LEP200 have been studied in [62, 63]. 
While the total cross section of e+e" — > W^W~ is extremely sensitive to deviations 
from the triple gauge boson couplings of the SM at energies high above the production 
threshold, the sensitivity to variations of this coupling in the LEP200 energy range is only 
at the percent level. Consequently theoretical predictions should be better than 1% to 
obtain reasonable limits on the structure of the gauge boson self interaction. 

The H^-pair production process allows an independent direct measurement of the W- 
boson mass with an expected accuracy of about 100 MeV [47]. Again this requires the 
knowledge of the total cross section with a precision better than 1%. 

Much effort has been made in recent years to obtain such precise theoretical predic- 
tions for l^-pair production. The virtual electroweak and soft photonic corrections were 
calculated by several authors [64, 65, 66, 67]. The complete analytical results for arbitrary 
polarizations of the external particles were published in [66]. These will be used for our 
evaluations. For the unpolarized case they numerically agree with those of [67] better 
than 0.3% and essentially also with [64]. The hard photon bremsstrahlung corrections 
have been evaluated by [68, 14, 69, 70] for definite initial and final state polarizations. 
The effects arising from the finite width of the W-bosons have been studied in the Born 
approximation in [71] and including the leading weak corrections in [72]. Recently also 
the hard bremsstrahlung to the process e+e" — > VF+l^" — > 4 fermions has been evaluated 
[73]. 

In the following we will review and update existing results for the virtual and real 
electroweak corrections and the finite width effects. We will add some new results on 
approximate formulae for the H^-pair production cross section. 

11.1 Notation and amplitudes 

We discuss the process 

e+(pi, a,) + e-(p2, (T2) ^ W+ih, Ai) + W'ih, A2). (11.1) 

The arguments indicate the momenta and hclicities of the incoming fermions and outgoing 
bosons (cTj = ±|, Aj = 1,0, —1). We introduce the usual Mandelstam variables 

t = {pi- hf = (p2 - k2)^ = M^- 2E^ + 2E^Pcos'&. 

Here E is the beam energy, 'd the scattering angle between the e~ and the and 
P = ^1 — My^/E"^ the velocity of the PF-bosons in the center of mass frame. The electron 
mass has been consistently neglected. 
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In the approximation of zero electron mass the invariant matrix element vanishes due 
to chiral symmetry for equal helicities of the and e~ . Consequently we can write 

A4((Ti, (72, Ai, A2, s, t) = M{a, Ai, A2, s, t) (11.3) 

with (T = o"2 = —<Ji. If wc neglect the CP- violating phase in the quark mixing matrix, CP 
is a symmetry of the process leading to the relation 

M{a, Ai, A2, s, t) = M((7, -A2, -Ai, s, t). (11.4) 

Consequently there are only 12 independent helicity matrix elements instead of 36. 
As discussed in Chap. 5 the general matrix element 

M{(T, Ai, A2, s, t) = v{pi, -a)M'"'u{p2, cr)e^{ki, Xi)e^{k2, A2) (11.5) 

can be decomposed into formfactors and standard matrix elements. Due to the above- 
mentioned symmetries we do not need all of the (overcomplete) 40 standard matrix ele- 
ments given in (5.14) for a general process of vector pair production in fermion-antifermion 
annihilation, but only seven for each fermion helicity 

M{a, Ai, A2, s, t) = E M^K^s, t). (11.6) 

i=l 

The standard matrix elements M.'^ are defined in (5.14) together with 

Ml^Ml, + Ml2, (11-7) 

Only six are linear independent. Using (5.15) we can express by the others 

M^T = --{M'', + M"^) + ^ — Ml + -Ml + Ml (11.8) 

Therefore there are only twelve independent formfactors. The standard matrix elements 
can be calculated with the methods described in Sect. 5.3. The explicit results can be 
found in [66, 74]. 

11.2 Born cross section 

At the Born level three diagrams contribute to 1^-pair production (Fig. 11.1). We 
omitted a Higgs-exchange diagram, which is suppressed by a factor rrie/Mw and thus 
completely negligible. The i-channel i^e-exchange diagram contributes only for left-handed 
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Figure 11.1: Born diagrams for e^e — > W^W . 



electrons whereas the s-channel diagrams containing the nonabehan gauge couphng con- 
tribute also for right handed electrons. The analytical expressions read 



^- + 2{M-,-M-,)e' 



Sw 



s-Ml 



9^2 



+e 



(M3 - M-,) 



s s- Ml 



;ii.9) 



Mo{+, Ai, A2, s, t) = 2{Mt - Mt)e^ 



Cw I 
- —9e- 



1 



s s-Ml 



s Sw"" s- M| 



2{Mt-Mt). 



where we have inserted the explicit form of the Z-boson fermion couplings g~ , (A. 14). 
The corresponding cross section for arbitrary longitudinal polarizations of the leptons and 
bosons is given by 



da 



^ E \{l-2aP+){l + 2aP-)\Mo{a,X,,X2,s,t)\\ 

Ai,A2 



;ii.io) 



and arc the polarization degrees of the leptons {P~ = ±1 corresponds to purely right- 
and left-handed electrons, respectively). 

The Born cross section determines the main features of ly-pair production. We first 
study the threshold behaviour [75, 76]. For small P the matrix elements behave as 



:il.ll) 



Consequently the s-channel diagrams vanish at threshold and the i-channel graph domi- 
nates in the threshold region. For P <^ 1 the total cross section is given by 



(7o(s) 



2 1 



S As^r 



4/5 + (/?'). 



;il.l2) 



All terms oc (3'^ which are present in the differential cross section drop out in the total cross 
section, s-channel diagrams yield contributions oc In the SM the coefficient of the 
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term in (11.12) is roughly equal to the one of the leading /3 term. As long as that coefficient 
is not enhanced drastically the f3^ term is negligible in the threshold region, i.e. for 2M\.y < 
y/s < 2Mw + 10 GeV (for = 2Mw + 10 GeV we have /5 = 0.33 and (3^ = 0.04). 
Consequently the shape of the total cross section close to threshold is completely governed 
by the linear rise in (3 and hence by kinematics alone. Any change in the couplings will 
only affect the coefficient of the /3 term and thus the normalization of the cross section. 
Moreover many new physics effects such as anomalous gauge couplings contribute to the 
s-channel only and thus do not affect the leading term. The inclusion of the finite width 
of the ly-boson smears the threshold considerably (see Sect. 11.5). However, since the 
next to leading term becomes only sizable several above threshold only the leading 
term is relevant for the cross section in the region of the nominal threshold. 

This fact allows a model independent determination of the W-mass from the Im- 
pair production threshold [75]. The measured cross section up to about 10 GeV above 
threshold is fitted with a three-parameter curve 

a{s)^^ + basMiMw,s) (11.13) 

where a/s accounts for the background, 6 is a model dependent normalization factor and 
asM the W^-pair production cross section in the SM depending on the PF-mass. Eq. (11.13) 
is valid including radiative corrections and finite width effects. 

At high energies the VF-pair production cross section is subject to large gauge cancel- 
lations arising from the contributions of longitudinally polarized 1^-bosons. For s ^ 
the matrix elements behave as 



— rsj 



t ' s M^' 



[U.U) 



1 . 1 , . .„ . M? 



but 

1 

IT I O/A/fcr A /<0■^ _. W 

1 1 



-Ml + ^-^2{Ml -Ml)^^-r 0(1), 

(11.15) 



2{Ml-Ml) -0(1). 



Consequently the Born matrix elements (11.9) have a good high energy behaviour. While 
the cross sections corresponding to the s-channel or t-channel diagrams alone violate 
unitarity at high energies 

T^O? 1 



the SM cross section respects it 

-KO? 1 



w 



s 4s^24M^' 



s N 1 1 5 



2Km:-i)-o- 



+ 



(11.17) 



The gauge cancellations arc illustrated in Fig. 11.2. They reach one order of magnitude 
at 400 GeV and two orders at 1 TeV. They only occur for longitudinal ly-bosons. After 
the gauge cancellations the i-channel again dominates the SM cross section. Compared to 
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Figure 11.2: Gauge cancellations in the total cross section for W-pair production. Shown 
are the Born cross sections arising from the s-channel cro,s and t-channel (Jo,t diagrams 
alone and the SM cross section cto. 

the ^-channel contribution all other contributions to the total cross section are suppressed 
by ~ 501og(s/M^). For example the cross section for right handed electrons is 



We now consider the complete Born expressions for the VF-pair production cross sec- 
tion. The differential cross section for the unpolarized case and for longitudinally polarized 
Vl/^-bosons is shown in Fig. 11.3 and 11.4, respectively. Due to the t-channel pole the un- 
polarized cross section is strongly peaked in forward direction at high energies and drops 
smoothly with increasing scattering angle. In contrast the differential cross section for 
longitudinal VF-bosons has a minimum for a certain energy dependent finite scattering 
angle 7^ vr. 

In order to show the importance of the separate contributions we give in Tab. 11.1 the 
integrated cross section for different center of mass energies and different polarizations of 
the leptons and 1^-bosons. 

The gauge cancellations depend crucially on the values of the SM couplings. Any 
deviations from these values can lead to sizable effects at higher energies since they are 
enhanced by a factor /5s/M^. This fact concerns especially anomalous three gauge boson 





and the cross section for longitudinal I^-bosons 




(11.19) 
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Figure 11.4: Lowest order differential cross section for the production of longitudinal 
1^-pairs at different center of mass energies. 
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165.0 


10.761 


21.466 


0.056 


5.254 


1.349 


4.158 


180.0 


16.928 


33.600 


0.256 


9.397 


1.743 


5.788 


200.0 


17.709 


35.070 


0.348 


11.189 


1.428 


5.092 


250.0 


15.011 


29.745 


0.277 


11.487 


0.731 


2.793 


500.0 


6.534 


13.019 


0.050 


6.178 


0.113 


0.244 


1000.0 


2.139 


1.867 


0.010 


2.395 


0.027 


0.017 



Table 11.1: Integrated lowest order cross section for different polarizations and different 
center of mass energies and M^r — 80.23 GeV. 

couplings which have been studied by many authors [77] . The sensitivity to these effects is 
best at high energies and large scattering angles where the t-channel pole is not dominant. 
Nevertheless one hopes to determine the anomalous couplings up to 20% at LEP200 [6]. 

Using right-handed electrons one could study a pure triple gauge coupling process, 
but this would require longitudinally polarized electron beams. Furthermore the right- 
handed cross section is suppressed by two orders of magnitude compared to the dominant 
left-handed mode, mainly because there is no t-channel contribution. On the other hand, 
nonstandard couplings or other new physics can enhance it drastically exactly for this 
reason. 

11.3 Virtual and soft photonic corrections 

The radiative corrections can be naturally divided into three classes, the virtual 
corrections, the soft photonic, and the hard photonic corrections. Since the process 
e^e~ W^W~ involves the charged current, the radiative corrections cannot be sepa- 
rated into electromagnetic and weak ones in a gauge invariant way. We first discuss the 
virtual and soft photonic corrections. 

11.3.1 Virtual corrections 

The virtual corrections get contributions from the i/g-, 7- and Z-self energies, the 
7-Z-mixing energy, from the vertex corrections to the 667-, eeZ-, eVgW-, WW^- and 
WW Z -vertices and from box diagrams. The necessary counterterms involve in addition 
the e- and W^-self energies. Altogether one has to calculate more than 200 individual 
diagrams. These can be treated using the methods described in the first part of this 
review. The number of generic diagrams to be evaluated is about 30. The results can be 
expressed in terms of the formf actors defined in (11.6) 

5Mi{a, Ai, A2, s, t) = E M^SF^r (11.20) 

i=l 
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The formf actors 5Ff can be evaluated for every CP-invariant set of diagrams separately. 
For CP-violating diagrams wc need in addition the standard matrix elements i— A^f 2! 

— A^4 2 and Aij I — ■M.7 2- These drop out in CP-invariant combinations. The 
explicit analytical results for the formfactors are given in terms of the scalar coefficients 
of tensor integrals in [66, 74]. The reduction to scalar integrals and their evaluation is 
done numerically using the formulae given in Chap. 4. The contribution of the virtual 
corrections to the cross section is given by 



Jda\ P 



J2 t(1 - 2aP+)(l + 2c7P-)2Re {M*o5Mi 



(11.21) 



Ai,A2 



The cancellations already present at the Born level occur as well at the level of radiative 
corrections. These cancellations only work for gauge invariant quantities. Consequently 
the inclusion of the leading higher order contributions must be done such that gauge 
invariance is respected. Otherwise one may introduce sizable unphysical corrections. This 
will be discussed in more detail in Sect. 11.3.3. 

The presence of these cancellations enforces very careful tests of the numerical stability 
of the computer programs. The reliability of the results is founded on agreement between 
independent calculations [66, 67]. 



11.3.2 Soft photonic corrections 

The soft photonic corrections can be easily obtained using the results of Chap. 7. The 
soft photon matrix element reads {k is the photon momentum) 



Ms = eMo 



ep2 epi ^ eki ek2 
kp2 kpi kki kk2 



This yields the soft photon cross section as 



with 
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(11.23) 
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Adding the soft photon cross section to the contribution of the virtual corrections 
(11.21) the IR- singularities cancel. Moreover also the large Sudakov double logarithms 
log^ (m^/s) drop out. 

11.3.3 Leading weak corrections 

In order to set up improved Born approximations which are often very handy the 
first step is to extract the leading corrections [78]. The universal corrections involving 
Aa and Ap can be easily obtained from (8.25) including the leading O(a^) contributions. 
There are no nommivcrsal corrections oc amf/M^/ to the W^-pair production cross section 
for not too high energies, i.e. as long as the unitarity cancellations are not sizeable. In 
the LEP200 energy region also terms involving logm^ or logM^ may become important. 
These have been evaluated in the hmit M^, ^ s. In addition close to threshold apart 
from the large bremsstrahlung corrections which will be discussed in the next section there 
is a sizable effect of the Coulomb singularity. This can be simply obtained from general 
considerations or to 0{a) directly from the loop diagrams involving photons exchanged 
between the final state 1^-bosons. Altogether this yields the following approximation 



m: 



9„2 



1 - Aa 1 + ^Ap 



a 1 / 1 1 "^f a 11 1 Mff aTT 



+ e 
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(11.25) 
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4:ns- Ml 
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All terms in (11.25) respect the high energy cancellations apart from those involving 
logm^ and logM^. However, these were obtained for M^, ml 3> s, whereas the unitarity 
cancellations work for s ^ Mfj, ml. In this limit the terms containing logm^ and log Mjj 

are absent. These may, however, cause large effects for small energies and large top quark 
or Higgs boson masses. This phenomenon was called delayed unitarity cancellation in 
[79]. Introducing Gp instead of e^/s^ and the running a{s) we obtain 
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We have included four functions Ci{s,t), i = 1,2, o" = ± in this approximation. These 
are necessary to describe the angular dependence of the differential cross section. The 
complete one-loop invariant matrix element for ly-pair production involves 12 formfactors 
F^. It turns out, however, that only four of them namely the Cf are relevant in the 
LEP200 energy region. For higher energies even Ci can be omitted. The functions 
Cf have been determined such that they reproduce the corresponding exact one-loop 
formfactors sufficiently well in the LEP200 energy region [78]. 

We want to stress that the naive summation of the Dyson series of the self energies 
may lead to incorrect results, i.e. a wrong high energy behaviour. This happens because 
the leading corrections are not only contained in the self energies but also in the vertex 
corrections. The actual place of their appearance depends on the choice of the field 
renormalization . 



11.3.4 Numerical results for the virtual and soft photonic corrections 

We now present some numerical results for the radiative corrections in the soft photon 
approximation. The numerical input parameters are defined in Sect. 8.1. The soft photon 
cutoff is chosen as AE/E = 0.1. Different choices of AE/E uniformly shift the absolute 
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Figure 11.5: Radiative corrections to the differential cross section relative to lowest order 
for the unpolarized case at various center of mass energies. 

value of the corrections but do not influence their angular dependence very much. Fig. 11.5 
shows the relative correction factor S defined through 

for the unpolarized case. While the variation with the scattering angle is relatively flat 
for LEP200 energies it becomes stronger with increasing energy. In the forward direction 
where the Born cross section is dominated by the t-channel pole the energy dependence is 
very weak. In the backward direction, however, the percentage correction varies strongly 
with energy and reaches large negative values up to —50% at 1 TeV. Nevertheless since 
the absolute cross section is small for large scattering angles (see Fig. 11.3), the relative 
corrections to the integrated cross section stay below 20% up to 1 TeV. Note that the 
one-loop corrections are large exactly in that region where the sensitivity to new physics 
is highest. 

The behaviour of the corrections for purely transverse VT-bosons is similar to the unpo- 
larized case. In contrast to this the corrections for purely longitudinal bosons (Fig. 11.6) 
exhibit a strong angular dependence arising from the minima in the lowest order cross 
section (Fig. 11.4). 

The sensitivity of the total unpolarized cross section on the unknown masses of the 
Higgs boson and top quark is illustrated in Tab. 11.2 and 11.3. A change of rrit from 

80 to 200 GeV affects the cross section by less than 3% apart from the region very close 
to threshold. The large effect close to threshold is due to the variation of M]y and thus 
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Figure 11.6: Radiative corrections to the differential cross section relative to lowest order 
for purely longitudinal ly-bosons at various center of mass energies. 

the variation of the threshold with m^. A variation of Mh between 50 and 1000 GeV 
influences the total cross section by less than 0.5%, again with the exception of the 
threshold region. Note that this is valid for constant a, Gp and Mz- Fixing instead Gf, 
M\Y and Mz the dependence on rrit is much weaker. This allows to determine My/ from 
the cross section practically independently on rrit and Mh as pointed out by Jegerlehner 
[80]. These results for the top and Higgs mass dependence remain valid if we include hard 
photonic corrections. 

Using the functions Cf given in [78] the improved Born approximation (11.26) agrees 
with the full one-loop order result within 0.5% for ^/s < 220 GcV and within 1% for 
\/s < 270 GeV in the case of the total cross section. For the differential cross section 
the deviation is at most 1% for a/s < 210 GeV. The largest difference occurs for large 
scattering angles, i.e. where the cross section is small. 

11.4 Hard photon bremsstrahlung 

11.4.1 Complete calculations 

The complete hard photonic bremsstrahlung to e'^e" — > W~^W~ was determined in 
[68, 70]. The polarized amplitudes for the process e+e~ — > W'^W~^ were calculated 
using three different methods. The first one, described in detail in [69] uses the Weyl 
representation for Dirac matrices and spinors and results in expressions for the amplitudes 
in terms of the components of momentum and polarization vectors in the center of mass 
frame of the incoming leptons. The second method used in [14] is based on the Weyl- van 
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80 


120 


160 


200 




79.87 


80.10 


80.37 


80.69 




cr/pb 


165.0 


10.120 


9.743 


9.218 


8.461 


180.0 


15.521 


15.620 


15.646 


15.626 


200.0 


15.944 


16.112 


16.220 


16.301 


500.0 


5.689 


5.760 


5.807 


5.847 


1000.0 


2.064 


2.088 


2.103 


2.113 



Table 11.2: Total unpolarized cross section for e^e — > W^W including virtual and soft 
photonic corrections for different top quark masses at various center of mass energies. 



Mn/GeV = 


50 


100 


300 


1000 


Mw/GeV = 


80.26 


80.23 


80.16 


80.06 


v/i/GeV 


cr/pb 


165.0 


9.488 


9.503 


9.614 


9.793 


180.0 


15.654 


15.638 


15.612 


15.598 


200.0 


16.168 


16.168 


16.113 


16.105 


500.0 


5.785 


5.785 


5.776 


5.764 


1000.0 


2.097 


2.096 


2.093 


2.088 



Table 11.3: Total unpolarized cross section for e+e — > W'^W including virtual and soft 
photonic corrections for different Higgs boson masses at various center of mass energies. 
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der Waerden formalism. It yields concise analytical formulae for the amplitude which 
are manifestly Lorcntz invariant. The relative numerical difference between both results 
was found to be less than 10~^ for the amplitude squared. In [70] the amplitudes were 
calculated numerically. 

Prom this the total cross section is obtained as 



a{s) 



r9 ^5, / C^^^iEl-Mr<^^'nPi+P2-^i-^2-^) 

(27r)5 2s J 2kio 2/c2o 2fco 4 ^ 



1 1 
2^8(27r) 



d cos ■j?2 d cos d4> dko 



ko\k2\ 



|k2|(\/s - ko) + kok2oC2o 



pol 



(11.28) 

where ■j?2, '& are the polar angles of the Vr~-boson and the photon, is the azimuthal 
angle of the photon with respect to the incoming electron and 



C20 = sin 'i?2 sin r) cos cp + cos '62 cos {}. 



;il.29) 



The nontrivial phase space integrations are performed using Monte Carlo routines. Thus 
experimental cuts can be easily implemented. 

Eq. (11.28) contains the soft photon poles. These are eliminated by a cut ko > AE on 
the photon energy. After combining soft and hard photonic corrections the cut dependence 
drops out. This has been checked numerically for AE/E between 10"^ and 10~^. 

11.4.2 Leading logarithmic approximation 

The leading logarithmic (LL) QED corrections to the W^-pair production cross section 
were already calculated in [43]. The resulting cross section is given by [68] 



o-ll(s) 



4M2 /s 



dz4>{z)aQ{zs), 



;il.30) 



where ao{zs) denotes the (improved) Born cross section at the reduced CMS energy 
squared zs. The flux ^(2;) reads 



<P{z) = 6{l-z) 
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TT 
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(11.31) 



+ (^l)' {5(1 - z) (^2V £ + 31og£ + ^ - ^) 



+ 0{l-e-z) 



ii^(2 log(l -z)-logz+^) + ^{l + z) logz -{1-z) 



where L = \og{Q^/ml) is the leading logarithm and e = AE/E the soft photon cutoff. 
(j){z) is given including 0(q;^) LL-contributions. Furthermore some nonleading terms 
are incorporated taking into account the fact that the residue of the soft photon pole 
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is proportional to L — 1 rather than L for the initial state radiation. The scale is a 
free parameter. It can only be determined through higher order calculations. In [68] the 
integral in (11.30) was performed numerically. Neglecting the 0{a'^) leading logarithms 
and the nonlogarithmic terms it has been evaluated analytically [78]. 

11.4.3 Numerical results 

The results presented in this section were calculated by [68, 81]. The parameters are 
the same as in Sect. 8.1. In the Born cross section (Tq, entering the LL approximation, a 
has been replaced by everywhere. Thus the large fermionic corrections are absorbed 
at least in the dominant t-channcl contributions. 

The total cross section is plotted in Fig. 11.7 in the LEP200 energy range. Shown 
are the Born cross section with a replaced by G^, the cross section including the full 
0{a) corrections and including in addition the 0{a^) LL corrections. The corresponding 
numbers are given in Tab. 11.4 for various CMS energies. The uncertainty of the full 0{q) 
result is due to the Monte Carlo integration of the hard bremsstrahlung corrections. This 
error refers also to the last column of Tab. 11.4. The 0{q.) LL results were evaluated for 
two scale choices 

The second one is motivated by the fact that the total cross section is dominated by the 
^-channel pole. It reproduces the exact 0{a) results better. The difference is found to be 
less than 2% for > 165 GcV. Choosing the scale = s the deviation from the exact 
0{a) result is about 5% at 165 GeV. Also at higher energies the scale choice = —tmm 
turns out to be preferable. It reproduces the complete 0{a) result including hard photon 
bremsstrahlung within 1% for 170 GeV < < 500 GeV. The effect of the ©(a^) LL 
contribution is demonstrated in the last column of Tab. 11.4. It reaches about 5% at 
165 GeV, decreases with increasing energy and is small for y/s > 190 GeV. A practically 
identical result is obtained by soft photon exponentiation. 

The large deviation between 0{a) LL and the exact result close to threshold is due 
to the Coulomb singularity (see Sect. 11.3.3). It amounts to about 10% at 1 GeV above 
threshold and is not included in the 0{a) LL result. Note also the large 0{a'^) correction 
close to threshold (28% at 161 GeV). 

11.5 Finite widths effects 

Reahstic calculations for 1^-pair production must include the decays of the 1^-bosons 
into fermions. These are especially important around threshold. 
In real experiments one observes the reaction 

e+e- W+W~ final states. (11.33) 

The W^-bosons give rise to peaks in the invariant mass distributions of the final state 
particles. Therefore one has to calculate the cross section for e~^e~ fihfsfAil, g)- 
This task has been attacked but not completed so far [83]. Above the VF-pair production 
threshold the dominant contributions come from Feynman diagrams containing resonant 



90 



20 



ff 




21 — 
ISO 



— — — Bom 

— — " incl. fuU 0(a) 

incl. m 0(a) mi LL 



170 



17S 



180 les 

V?/GeV 



190 



195 



200 



Figure 11.7: Total cross section for 1^-pair production including hard photonic correc- 
tions. 
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161.0 


4.411 


2.003 


2.158 


2.556 ± 0.002 


3.255 


165.0 


11.761 


8.141 


8.429 


8.553 ± 0.006 


9.089 


170.0 


15.465 


11.967 


12.285 


12.264 ± 0.010 


12.606 


175.0 


17.413 


14.264 


14.578 


14.484 ± 0.013 


14.690 


180.0 


18.501 


15.730 


16.028 


15.920 ± 0.014 


16.033 


190.0 


19.361 


17.272 


17.525 


17.373 ± 0.016 


17.375 


200.0 


19.351 


17.810 


18.015 


17.796 ± 0.017 


17.712 


250.0 


16.106 


16.223 


16.257 


16.033 ± 0.023 


15.937 


300.0 


13.473 


13.734 


13.682 


13.543 ± 0.026 


13.470 


500.0 


7.142 


7.664 


7.519 


7.449 ± 0.019 


7.430 



Table 11.4: Total cross section for e'^e 
tions. 



Vl^"*"!^ in pb including hard photonic correc- 
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Figure 11.8: General structure of diagrams containing two resonant VF-propagators. 



H^-propagators. The cross section for W^-pair production is obtained by calculating all 
diagrams containing two resonant W^-propagators. Diagrams contributing to the same 
final state but without two resonant propagators are considered as background. They are 
suppressed by a factor M^^/r^ 40 if the full range of invariant masses ^/sl of the final 
state particles is included. If ^/s^ is restricted by a cut A 



Miy - A < < Mv^^ + A, (11.34) 

the suppression is even M^/(r^A) ^ 300 for A ^ 10 GeV. Explicit calculations show 
that the background contributions are below 1% for ^/s > 2My\r [83]. It becomes, however, 
more relevant below the nominal threshold. But even above threshold nonresonant Born 
contributions to e'^e~ — > f\f2hfi{l-i g) must be taken into account to obtain an accuracy 
of better than 1%. 

There are three types of diagrams which may give resonant contributions. The most 
important ones are factorizablc diagrams with the structure shown in Fig. 11.8 which 
evidently contain two resonant H^-propagators. The corresponding cross section is given 

by 

a{s) = dsids2a*{s,si,S2)p{si)p{s2)9{y/s - ^/s^ - ^/s^), (11.35) 

J(Mvv^-A)2 

where a*{s, Si, S2) is the 'cross section' for the production of two off-shell 1^-bosons and 

with the 'decay width' T^{s) for an off-shell PF-boson. Note that 

p(s) 5{s - M^) for ^ 0. (11.37) 

The off-shell quantities a* {s, 81,82) and r^(s) are not gauge invariant. However, the 
leading resonant contributions to a{s) are. Eq. (11.35) closely resembles a Breit-Wigner 
approximation for the unstable VF-bosons. In the threshold region a* {8, 81, 82) depends 
strongly on Si and 82- Consequently (j(s) deviates considerably from (7*{s, M^^, M^), 
the cross section for on-shell stable W^s. Fig. 11.9 [81] shows this effect in lowest order 
and with the full 0{a) corrections to a* and r^(s) included. This dependence is mainly 
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Figure 11.9: Total cross section for 1^-pair production in lowest order and including the 
full 0{a) corrections with and without finite width effects. 



due to the threshold factor fi;^/^(s, si, S2) contained in a*{s, si, 52)- Extracting this factor, 
a*/K}^^ depends only weakly on si and S2- This is also the case for V^{s)/^/s with 
respect to s. Replacing these quantities by their on-shell values we find the following 
approximation 

with 



(11.38) 

*) = I .^yi'SL.,,.,. . (11.39) 



This approximation is gauge invariant because a* {s, Ml,, M^) and r'^(M^) are physical 
on-shell quantities. It is particulary useful above the nominal threshold, whereas it gets 
worse below threshold because there at least one of the 1^-bosons has to be off-shell. 

For high energies (s ^ ^w) also k{s, Si, S2) varies only weakly with si and S2 in the 
resonance region S2 ~ M^,. Replacing it by its on-shell value we can perform the 
integrations and obtain 

(7(s) Ri(7*(s,M^,M^), for s>Mw. and A > F^. (11.40) 

Eq. (11.35) incorporates all resonant lowest order contributions and all one-loop correc- 
tions associated either with the production of l^-pairs or the decay of the l^-bosons (given 
in Chap. 9). This includes in particular all self energy and vertex corrections and thus all 
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Figure 11.10: Examples for additional diagrams leading to resonant contributions. 

leading corrections. A more thorough analysis of this kind has been carried through for 
the case of Z-pair production [82]. 

Feynman diagrams which do not fit into the structure shown in Fig. 11.8 give nonres- 
onant contributions and can thus be neglected with the exception of two classes. Both 
originate from photonic corrections. The first type results from virtual photons exchanged 
between the external hues connected to different blobs in Fig. 11.8. An example is shown 
in Fig. 11.10a. These diagrams give rise to resonant contributions coming from photons 
which are nearly on-shell. From similar cases in ^-pair production wc know that these are 
cancelled by the corresponding bremsstrahlung diagrams if one integrates over the whole 
photon phase space. For stringent cuts, however, resonant contributions survive. 

The second type of diagrams consists of those where a real photon is emitted from 
the internal W^-boson line (Fig. 11.10b). There are three Vt^-propagators in the diagram. 
If the photon is hard the corresponding resonances appear in three different regions of 
phase space. Therefore these diagrams seem not to fit into the simple Breit-Wigner-like 
picture discussed above. 

In order to take into account these resonant contributions properly one has to calculate 
the virtual and real photonic corrections to e^e~ — > 4 fermions. This has been done for 
real photon radiation [73]. Evaluation of the virtual photonic corrections is under way 
[83]. 
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12 Conclusion 



With the electroweak standard model (SM) we have a theory that describes all known 
experimental facts about the electroweak interaction. It has succesfuUy survived all preci- 
sion experiments at low energies and at LEPIOO. The upcoming experiments at LEP200, 
HERA and the planned hadron colliders will allow to investigate sectors of the SM which 
were not directly accessible so far. For a conclusive confrontation of future experimental 
results with the SM precise predictions are mandatory. 

For the adequate theoretical description of the experiments at LEPIOO the calculation 
of radiative corrections was inevitable. Although experiments outside the Z-region will 
not profit from the presence of a resonance, the expected experimental accuracy will be 
such that radiative corrections will be indispensable. Moreover radiative corrections allow 
to obtain information on otherwise not accessible quantities such as the mass of the top 
quark or the Higgs boson. 

One of the next important classes of experiments will be the investigation of the W- 
boson and its nonabelian couphngs at LEP200. We have presented the relevant formulae 
necessary for the corresponding higher order calculations. Together with [66, 14] these 
cover the complete analytical expressions for processes with on-shell 1^-bosons. The one- 
loop virtual corrections are settled for the polarized differential and total cross section. 
Also hard bremsstrahlung has been calculated by several authors. We have given an 
improved Born approximation including the leading two-loop contributions for the to- 
tal and differential cross section. The effects of the finite width of the VF-bosons have 
been discussed for the lowest order cross section and the cross section including radia- 
tive corrections. While the inclusion of the non-photonic corrections is simple the correct 
simultaneous treatment of photonic corrections and finite width effects involves nonfac- 
torizable box diagrams. These contributions are under consideration. 

We have discussed the total and partial VT-decay widths including all one-loop and 
leading two-loop corrections. Because the W^-boson decays only into light fermions the 
widths can be described by a very simple expression with an accuracy better than 0.6%. 

Furthermore we have given results on the i-quark decay width. Also in this case the 
electroweak corrections can be incorporated into a simple approximation valid for a top 
mass below 250 GcV with an accuracy of about 1.7%. 

Apart from giving these explicit results we have discussed many techniques needed for 
the calculation of one-loop corrections. We have compiled a comprehensive set of formulae 
which are relevant for the calculation of one-loop radiative corrections within and outside 
the standard model. We have listed the complete set of Feynman rules for the electroweak 
standard model including the counter terms. These were expressed by the self energies of 
the physical particles in terms of two-point functions. We have given explicit results for 
the scalar N-point functions for = 1, . . . , 4 and the relevant formulae for the reduction 
of the higher scalar functions and the tensor functions to those. Furthermore we have 
outlined a general strategy for the calculation of one-loop diagrams. Finally we have given 
the general expressions for the soft photonic corrections. 

If the SM will prove to describe the upcoming experimental results succesfully, further 
precision checks will become necessary and a lot more calculations of radiative corrections 
will be required. These calculations will be even more involved than the existing ones 
because the structure of the corresponding physical processes will in general be more 
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complicated. The techniques and formulae presented in this review are general enough 
to serve as a basis for the evaluation of radiative corrections to reactions which will be 
studied at future colliders such as gauge boson scattering processes (W^"*"!^" W~^W~), 
electron photon reactions {e'j — > I'eW), reactions with three or more final state particles 
and so on. 

The methods described here have been implemented in the computer algebra package 
FEYN CALC . Many of the quoted formulae are included in this package. We hope that 
this compilation together with the packages FEYN CALC and FEYN ARTS can serve 
as a useful tool, facilitating future calculations of radiative corrections. 
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A Feynman rules 



In this appendix we hst the Feynman rules of the SM in the 't Ho oft- Feynman gauge 
including the counter terms in a way appropriate for the concept of generic diagrams. I.e. 
we write down generic Feynman rules and give the possible actual insertions. We omit 
any field renormalization constants for the unphysical fields. For brevity we introduce the 
shorthand notation 

C = C^y, S ^ Sw (A.l) 

In the vertices all momenta are considered as incoming. 
Propagators: 

for gauge bosons V — ^, Z, W in the 't Hooft Feynman gauge (^^ = 1) 





k 


for Faddeev-Popov ghosts G = , , 






for scalar fields S = H, x, 4> 




k 

s • — — •s 


and for fermion fields F 


= h 




F. ? .F 













i 


k^ 






i 




-M|' 




+ mp) 


p2 





In the 't Hooft-Feynman gauge we have the following relations 

M„7 = 0, M„z = M^ = Mz, 
Tadpole: 



w- 



W-counterterm: 



V.,,k 



Cik^ - C2 



with the actual values of Vi, V2 and Ci, C2 



w+w- 


: Ci 


= SZw, 


C2 


- M^5Zw + 5M^, 


zz 


: C^ 


= SZzZ: 


C2 


= Ml5Zzz + 5Ml, 


AZ 


■ Ci 


— I^Zaz + ^SZzA, 


C2 


= Ml\5ZzA, 


AA 


: Ci 


= SZaa, 


C2 


= 0. 



(A.2) 



(A.3) 
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SS-counterterm: 



Si, k 



-X- 



^2 



CiA;2 - C2 



with the actual values of Si, S2 and Ci, C2 



HH 


: Ci 


^5Zh, 


C2 


= M]j6Zh + 5Ml, 


XX 


: Ci 


= 0, 


C2 


— e 5* 1 X/\f2 


0+0- 


: Ci 


= 0, 


C2 


— e 5t 1 X!i#2 



(A.4) 



FF-counterterm: 



-X 



Cij^a;_ + Cr^uj+ -CgUJ_- C~^uJ+ 



fjft 



with the actual values of Fi, F2 and Cl, Cj?, , Cg 

{C, = \ {6Zif + 6Zlf^ , Cn = \ {SZif + 6Zif) , 

VWV-couphng: 



> 


< 







with the actual values of T^, V2, V3, V4 and C 



W^+W^+W^-W^- : C 

W+W-ZZ : 

ly+T^-AZ : 

W+W-AA : C = 



1 + 26Ze - 2 



Ss 



26Z, 



w 



s 



1 + 25Ze - + SZw + SZzz\ + ISZaz, 
l + 26Ze-^'-^ + 6Zw + ISZzz + \5Zaa 



s 

— \5Zaz — h^^ZzAi 



2"^A^ 2s^ 

1 + 25Z, + 5Zw + SZaa] + -SZzA- 



(A.5) 



(A.6) 
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VW-coupling: 



-ieC 



gnu{k2 - ki)p + Qupik'i - k2)i, + gp^,{kl - ks)^ 



with the actual values of Vi, V2, V3 and C 
AW+W~ : C 
ZW+W- : C^-l(l + 5Ze-y-^ + 5Zw + l5Zzz) + l5ZAz. 



1 + SZ, + SZw + 'jSZaa - kl^ZzA, 



SSSS-coupling: 




with the actual values of 5*1, 5*2, S3, S4 and C 
HHHH 



HHxx 
HH(f)(f) 



c = 

c = 
c- 



3 Mj 



3 Ml 
2^2 



1 + 25 Z^ 



^"7 + -TIT- -r 



2s MwM'jr 



1 + 2(5Ze - 


25s 
s 


+ 




+ 


e 5t 


M'i 


2s MwM'jj 


1 + 25Ze - 


2^ 

s 


+ 


5M% 


+ 


e 5t 




2s MwM'jj 


1 + 25Ze - 


25s 
s 


+ 


5M% 


+ 


e 5t 




2s MwM'jj 


1 + 25Ze - 


25s 
s 


+ 


5M% 


+ 


e St 




2s Mvt^Ml^ 



SSS-coupling: 




with the actual values of 5"!, S2, S3 and C 
HHH 



HXX 
H4)(t) 



C = 



3 Mj, 



2s Mvv' 



1 Ml 
2s Mvv^ 



1 + 



5s 
s 



5s 
s 



5A4 _e 



e 5t 
2s MvkM^ 

5f 

2s MhTmJ 



p h -^OZjh 



2 M,i 



2 
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WSS-coupling: 



T A 

> 


c 




^2 



with the actual values 

ZZ(t)+(t)- 
ZA(l)+(j)- 



AA(f)+(f)- 
ZZHH 

ZZxx 
ZAHH 1 
ZAxx i 

W^A(I)'^H 

W^Z(j)^X 
W^A(t)^X 



of Fi, 1^2, -^i, ,^2 and C 

1 + 25Ze - 2^ + + 5Z 

1 + - 2^ + 5Zy 



- 2s2 



H 



c 



2^2 c2 



1 + 25Ze + (^2.^2)^2 ■ 



+ 5Z 



zz 



1 + 25Ze + 



c 

c 
c 

c 

c 
c 



1 + 25Z, + 5Zaa 



5Z 



ZAi 



1 



1 

2s2c2 

1 \5Z 



1 + 2bZ, + 2^^ + bZzz + 



1 + 2(5Ze + 2^ 



5s 



+ 5Z 



zz 



2i2^2"^2^' 

1 + 2(5Ze - f + i^Zi^ + \8Zh + 
1 + 25Ze - f + + \bZn + I^Zaa - Yc¥ZzA, 



"2c 



"2s 



C = 1 + 2<5Ze - ^ + i^Z^ + UZzz T ^5^^^z, 



2s 



1 + 25Ze 



5s 
s 



+ \^Zw + \SZaa 



2s 2" 

±1 
2c 2 



(A.IO) 
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VSS-coupling: 





Si,ki 














= ieC{ki - /C2)^ 




S2,k2 





with the actual values of V, Si , S2 and C 



AxH 


C 


ZxH 


c 


A(f)+(f)- 


c 


Z(f)+(f)- 


c 




c 




c 



2cs 



1 + SZ, + + ^6Zh + UZi 



l + SZe + ISZaa + 



1 + 5Z,+ 



-C' 1 



2sc 2 



SZ 



ZA 



2sc 



5s 



+ '^SZzz 



— ^SZaz, 



1 + - ^ + UZw + \5Zh 



1 + SZe- 



5s 



+ 16Z 



w 



SW-couphng: 







s 






V2,u 



with the actual values of 5", Vi, V2 and C 





: C = 


Mw\ 


1 


+ SZ^-'-^ + 1'-^ + 16Zh + 6Z^' , 






HZZ 


: C = 


Mw^ 






HZA 


: C = 


Mw^ 


I^SZzA, 






(t)^W^Z 


: C = 


-Mw 


s 

c 


l + 5Ze+^'f + 1^ + ISZw + ISZzz 




- Mw 


(t)^W^A 


: C = 


-Mw 


1 + + + ^SZw + ^SZaa] - Mwl\ 


SZzA 
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VFF-coupling: 




with the actual values of V, Fi , F2 and C"*", C 



ififj 



Zfifj 



W djUi 



where 



+ dijgjldZzA, 



C+ = 0, 
C+ = 0, 



c- 



1 

\/2s 



'i^i, (1 + - f + + 

C- = ^5,, [1 + - f + i^Zn. + t + 6Z\f^) 



9f 

9J 



5Z, + ^'-f 



/3 



(A.13) 
(A. 14) 



5Z,+ 



+ 5gj. 

The vector and axial vector couplings of the Z-boson are given by 



2sc 



(A.15) 
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SFF-coupling: 




ie{C-uj^ + C+UJ+) 



with the actual values of S, Fi , F2 and C+, C 



Xfifj 



(f) djUi 



Iji^i 



c- = 

(:;+ = 
c- = 

c+ = 
c- = 

c+ = 
c- 
c+ 
c- 



1 1 

"2s Mw 



1 1 

2s 



<5m# 



^2s^-'W^,/MvK 



Sm 



+ \5Zh) 



5ijmf^i(l + 5Ze 



Ss 



SM: 



S ™-f.i 

+ l{mf,,6Z(f + 5Zij'^^mfA 



^2s^^WJ Mw 



1 1 



t/2s Mw 



5Mw 
Mw 



rud,: 



+1 Eki^Zi^^^Vkjrudj + Vikmd,kSZ^f] 



1 1 



V2s Mw 



1 1 



'mu,iVij(l + SZe 



5s_ _|_ Smu,i 



SM: 



Mw 



+muASVii 



d,L\ 
kj ) 



V2s Mw 



5Z^-^ + 



SMw 
Mw 



+5Vjimu,i 



iT.ki^Z'^k^Vlm. 



-/2s Mw 



y]k^u,k5Zl]^) 



jk 



1 "ij.i X 
■ a/2s Mvi, 



+1 EA:(<^^^^m,,,4 + m,,,-KUz,"f ) 



0, 
0, 



1 m, , 



^s M^y ^ij 



l + 6Z^-^ + 

s mi,, 



5Mw 
Mw 



ii 



(A.16) 
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VGG-coupling: 




with the actual values of V, Gi , G2 and C 

± 





C 




. c 




: C 




C 




C 




. c 



1 + 5Z, + \5Zaa\t-A5Zza, 
l + 5Z,-y-f + \5Zzz\±\5ZAz 



s 



l + 5Ze-y-^ + l5Zw 



T 

T 
± 



l + 5Z,-^^ + ^5Zw 



l + 6Ze + UZ. 



w 



1 + SZ,+ \5Zw 



SGG-couphng: 




ieC 



with the actual values ol S, Gi , 6*2 and C 

-^.Mw\l + 5Ze + + + },5Zh 





C 




c 




c 




c 




c 




c 



2 M,i 



2 



2sc 



1 + 6Z, 



s 



2 



2sc 



l + 5Z,+ 



2 

5s , iSMf 



(s'^—c'^)c'^ s 



+ 



2 



1 + 



2 A^i^ 
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B Self energies 



In this appendix we list all self energies of the physical fields. 
The gauge boson self energies read 

St^(^') = -^{lT.NL2Q}[-{e + 2m%)Bo{e,mf,,mf,) 

^ f,i 



(B.l) 



+ 



Se + 4M^ Bo{e, Mw, Mw) - 4M2,5o(0, Mw, Mw) L 



1 



3swCw 



+ -)e + (12c^ + A)M^]Bo{e, Mw, Mw) 



- {12clr - 2)M^Bo(0, Mw, Mw) + -k^ 



(B.2) 



^r'ik') = -^{^E^^{((5;)' + (97)') [-{k' + 2m%)B,{k^mj,,m,,] 

1 1 3 1 

+ 2mliBo{0, mf,i, ruf^i) + -k'^ + mliBo{k'^, mf,i, mf,i) \ 



+ 



1 



(184 + 2c^ - 2)^' + (244 + ISc^VF - 10)M^j5o(A;2, Mw, Mw) 



(244 - 84 + 2)M^Bo{0, Mw, Mw) + (4^ - l)-k^ 



+ 



J-^ j (2M| - 10M| - e)Bo{e, Mz, Mh) 

- 2M|So(0, Mz, Mz) - 2M|So(0, Mh, Mh) (B.3) 

- \2 " {Mk^ Mz, Mh) - So(0, Mz, Mh)) - -k' 
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+ n^liBoiO, mi^i, m^) + (5o(A;^ 0, m^) - Bo(0, 0, mi,i)) 



+ 



2A;2 



, rudj) - Bq{0, mu,i, maj)) 



-\ (2M^ + 5k')Bo{k', Mw, A) - 2M^5o(0, M^., M^.) 



(Bo{k^, Mw, A) - Bo{0, Mw, A)) + 



+ 



(40c^ - l)k^ + (16c^ + 54 - 10c^^2)M^]So(A;2, M^, Mz) 

- (16c^ + 2) [m^So(0, Mi^, Mw.) + M|5o(0, M^, Mz)l + (4^ - 



- {^cl + 1) 



A;2 



2\2 



(5o(A;2, Mw', Mz) - 5o(0, M^, Mz)) | 



+ (2M^ - lOM^, - A;2)5o(A;^ Mw, Mh) 

- 2M^Bo{0, Mw, Mw) - 2M^Bo(0, Mh, M„) 

- ~ {Bo{k', Mw, Mh) - Bo{0, Mw, M^)) - 
For the self energy of the physical Higgs boson wc obtain 

{2 
E ^c^T?&T l^Mmu) + {^m% - k')B,{k\ m/,,, m/,,)' 



(B.4) 



2s2 



(6M^ - 2e + ^)^o(A;^ M^, Mi^)+(3 + £^)MMw) - 6M, 



M? 



2 



4c2 „2 



(6M| - 2A;^ + ^)Bo{e, Mz, Mz)+{3 + ^) A(Mz) - 6M 



Mi 



2M| 



2M|- 
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w 



/\/f4 J\f2 



The fermion self energies are given by 



(B.5) 



a 
4n 



+ 5^. 



1 Imj^i 



Bi{p'',mf,i,Mz)+Bi{p^ Mh) 



(B.6) 



a 



2Bi{p\mf,i,X) + l 



+ Sij{g+f 2B^{p\mf,,Mz) + l 



2 r 



1 lm'f^ 
^^'^2^2 Ml 



Bi{p^,mf,i,Mz)+Bi{p^ Mh) 



+ 



2s^r My^r ^ 



(B.7) 



4So(p^m/,,,A)-2 



4Bo{p\mf,i,Mz)-2 



1 irn^,^^ 



Bo{p^mf,i,Mz)-Bo{p^ 



+ ^ E VikV^j'^Mp', mf,,k, Mw) 



w 



(B.8) 



/' is the isospin partner of the fermion / and Nq the colour factor, i, j, k run over the 
fermion generations. For down- type quarks VikV^j has to be replaced by V^^Vkj. 
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The two-point function Bq was given in Sect. 4.3. For Bi we find 



2 2 

Bi{p'^,mo,mi) = ^ ° (.So(p^,mo,mi) - ■Bo(0,mo,mi)) 

V ^ ^ (B.9) 

For the field renormahzation constants one needs in addition the derivatives of the 
self energies with respect to k'^ or p^, respectively. These are easily obtained from the 
expressions above. was given in Sect. 4.3, can be calculated from (B.9) as 



— 5i(p^mo,mi) = i^-^(5o(p^"^o,"^l) - ^o(0,mo,mi 

2 2 2 o 



r5o(p ,mo,mi). 



These derivatives become IR-singular for = and = or vice versa. This leads 
to IR-singular contributions in the field renormahzation constants of charged particles 
arising from photonic corrections to the corresponding self energies. Because these reduce 
to very simple expressions we give the photonic contributions to the field renormahzation 
constants of the PF-boson and the charged fermions explicitly 



«w-|photo„ic = log + ^ (i + 5(A + 1 - log ) . (B.11) 



^ a I photonic ~ a I photonic 



.2 



A - log + 4 + 41og 



(B.12) 
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C Vertex formfactors 



The vertex formfactors V, W, X, can be expressed by the scalar one-loop integrals 
BqIttiI, Ml, M2), Co(mf , ml, m\, Mq, Mi, M2) and the scalar coefficients of the vector and 
tensor integrals Bi{ml, Mi, M2), Cit^j){m\, ml, ml, Mq, Mi, M2) 



Ve 



V 



ml, ml, ml, Mq, Mi, M2) = Bo{ml, Mi, M2) - 2 - (M, 



M^)Ci (C.l) 



-{M^ -ml- Mi)C2 - 2{ml -ml- ml){Ci + C2 + Co), 

ml, ml, ml, Mq, Mi, M2) = 3Bo{ml, Mi, M2) + 4.MlCo 
+ (4m^ + 2ml - 2^o + ^0 - ^i)Ci 
+(4m^ + 2ml - 2^o + ^0 - ^1)^2, 



ml. 


ml. 


ml 


ml. 


ml. 


ml 


ml. 


ml. 


ml 
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mo, 


2 

m2 


ml. 


ml. 


ml 


/ 2 


2 

mo. 
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m2 


m^, 


m^. 


m2 


ml. 


ml. 


ml 


ml, 


ml, 


ml 


ml. 


ml. 


ml 


ml. 


ml. 


ml 


ml. 


ml. 


ml 


ml. 


ml. 


ml 



Mo, 


Ml 


M2) 


= SmlCo, 




M 


M \ 
IVI2) 


2 2 


Mo, 


Ml 


M2) 


= ml{Ci-Co), 




Ml 


1V12) 


ml 


Mo, 


Ml 


M2) 


m^Co + m2C2, 


Mo, 


Ml 


M2) 


= rriiCi, 


Mo, 


Ml 


M2) 


= 2(Ci + C2 + Co), 


Mo, 


Ml 


M2) 


= -2Co, 


Mo, 


Ml 


M2) 


= 3{Ci + C2), 


Mo, 


Ml 


M2) 


— 3Co, 


Mo, 


Ml 


M2) 


= ^j^2^[Boiml,Mi,M2) 


Mo, 


Ml 


M2) 


2 r 1 

= M2 + "^2*^2 , 


Mo, 


Ml 


M2) 


= -C2, 



l-M^{Ci + C2) 



(C.2) 

(C.3) 

(C.4) 

(C.5) 

(C.6) 

(C.7) 
(C.8) 
(C.9) 
(C.IO) 
(C.ll) 
(C.12) 

(C.13) 

(C.14) 
(C.15) 
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m 




ml, 


ml. 


Mo, 




[ml 


ml, 


ml. 


Mo, 




[ml 


ml. 


ml. 


Mo, 




[ml 


ml. 


ml. 


Mo, 


K 


[ml 


ml. 


ml. 


Mo, 




[ml 


ml. 


ml. 


Mo, 


K 


[ml 


ml. 


ml. 


Mo, 


K 


( 2 


2 

mo, 


2 

m^. 


Mo, 




[ml 


ml. 


ml. 


Mo, 




[ml 


ml. 


ml. 


Mo, 


K 


[ml 


ml. 


ml. 


Mo, 


xl 


[ml 


ml. 


ml. 


Mo, 




[ml 


ml. 


ml. 


Mo, 




[ml 


ml. 


ml. 


Mo, 



Ml 
Ml 
Ml 
Ml 
Ml 
Ml 

Ml 

Ml 
Ml 
Ml 

Ml 
Ml 
Ml 



M2 
M2 
M2 
M2 
M2 
M2 

M2 

M2 
M2 
M2 

M2 
M2 
M2 



= 4 
= 2 
= 6 



M'w 
-Ci, 

— C2 — Co, 

-4 [Cii + C12 + 2Ci + C2 + Co] , 
Ci + C2 + Co 



2Cii + 2C12 — C2 
Ci + C2 



-2 



mo 



Ml 



C22 + C12 



2C2, 



mi 
M^ 



Cii + C12 + Ci 



C22 + C12 + C2 , 



mf 
M^ 



Ci + C2 + Co 



-2Ci. 



(C.16) 

(C.17) 
(C.18) 
(C.19) 
(C.20) 
(C.21) 
(C.22) 

(C.23) 

(C.24) 
(C.25) 
(C.26) 

(C.27) 
(C.28) 
(C.29) 



Using the reduction methods decribed in Chap. 4 the vector and tensor coefficients can 
be expressed by scalar integrals. For illustration we give the explicit reduction formulae. 
The vertex function is defined as 



C... = C...{pi,p2, Mo, Ml, M2) ^ C...{ml,ml,ml Mo, Ml, M2) 



(C.30) 



[^2 _ Mi][{q+pif - Ml][{q + p2? - M|] 



with 



2222 -^/22 2\ 

= P2 = ^2> P1P2 = -2(^0 -^2) 



(C.31) 
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For the three-point vector functions (4.18) yields {P — 1, M — N — 1 — 2) 

r* T^3 D3,fc' 

^k — J-k — 1^2 jkk'ti 

with /c, /c' = 1, 2 and 



\{rn\ + ml — ml) 



l{mj + m2 - mg) 



ml 



Evaluating X2 ^ this gives 



2\ d3,2 



where 



K = K{ml, m^, ml), 
from (4.28). The i?'s are obtained from (4.19) as 



R 



3,2 



Bo{ml Mo, M2) - {m\ - Ml + Mo')Co - B^{ml M2, Mi) 
So(m?, Mo, Ml) - (m^ - M| + M2)Co - So(m^, M2, Mi) 



The tensor coefficients are evaluated analogously as (P = 2, M = 2) 

1 



00 



^3,00 _ ^3,1 _ ^; 



3,2 



L)-2 

Cki = = (^2 ^)kk'[R^''' — Ci 



ooj 



or more explicitly 

Coo 



1 r 



4 



Bo{ml, M2, Ml) + (M^ - M^ + m^)Ci 



+ (M2 - M^ + ml)C2 + 1 + 2M2Co 



G 



21 



C'll — 

4 

<-^12 — 5- 

^ _ 4 

Ooo — ;t 



m2(-R?'^ - Coo) + :^{'ml-ml- ml)Rl'^ 



m^-mj- ml){Rl'^ - Coo) + mfRi 



2d3,2 



"^2-^2'^ + |("^o - "^1 - mi){R2'^ - Coo) 



,2wp3,2 



-{ml -ml- ml)R^2:^ + mi(i?2'^ - Cqq) 
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with 



= M^Co + 5oK, M2, Ml), 

1 

2 
1 



R 



,3,1 



?3,2 _ -L 

'^^ "2 

,3,2 _ 1 
~ 2 







- Ml + M2)Ci - B^{ml M2, Ml)' 






fii (ml, Mn, Mo) 


V 1 


- M^ + Ml)C2 + (So + Si)K, M2, Ml)' 




.BiK,Mo,Mi) 


- (m^ 


- M| + M2)Ci - B^{ml M2, Ml) 








- (m^ 


- Ml + M2)C2 + (So + Si)(m2, M2, Ml)" 


. (C.39) 



Note that C12 = C21 can be calculated in two different ways. In the evaluation of Coo we 
used (4.55) 

Bi was given in (B.9). The results for the scalar integrals can again be found in 
Sect. 4.3. 
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D Bremsstrahlung integrals 



For the decay width of a massive particle with momentum po and mass tuq into two 
massive particles with momenta pi, p2 and masses mi, 1712 and a photon with momentum 
q and mass A we need the following phase space integrals 



1 



(mo, nil, 1712) = -ij 



d^pid^P2d^q __^{±2qpjJ---{±2qpjJ, 



2pio 2p2o 2go 



5(Po-Pi-P2-?) 



{±2qp,,)---{±2qpJ 



■ (D-1) 



Here j^, ii = 0, 1, 2 and the plus signs belong to Pi,P2, the minus signs to po- 
Introducing the abbreviations 



K — n{m^, ml, ml). 



(D.2) 



as defined in (4.28) and 
with 



ttIq — ml — ml + K 
2mim2 ' 
ml — ml + m| — k 
2mom2 ' 



/32 



m5 + mT — mo — n 



2momi 



(D.3) 



(D.4) 



= 1, 

we get compact expressions for the final results. From (D.l) it is evident that the integrals 
with the indices 1 and 2 interchanged are obtained by interchanging mi and m2- We list 
only the independent integrals. The IR-singular ones are given by 



-^00 — 



4mo 



1 



/tlog( 



K 



'11 



hi = 



Amlml 
1 



Amomim2 
2 



log( 



) -K - (ml - ml) log(^) - ml log(/3o) 
Xmomi^^y^ - - "^^^ log(|)-m?log(A) 



(D.5) 
(D.6) 



Ami 



21og( 



Amomim2 



)log(/?2) + 21og2(/32) - log2(/?o) - log'(/?i) 



+2Sp{l - Pi) - Sp{l - p'o) - Sp{l - p!) 



(D.7) 



(12 



'01 ~ -'02 



4mo 



-21og( 



Amomim2 



)log(/3o) + 2V(/3o) - V(/3i) - V(/32) 



+25p(l - p',) - Sp{l - P!) - Sp{l - Pi) 



(D.8) 



For the IR finite integrals we obtain 



1 

4mo 



mo + m^ + ma) + 2momi log(/32) + 2mom2 log(/3i) + 2m-jn2 log(/3o) 
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4ml 



Ami 



Ami 



-2m2log(/?2)-2m2log(A)-«:J, 
-2m2log(/?2) - 2m2log(/?o) - « , 
m{ log(/32) - ml{2ml - 2ml + mj) log(/3i) 



K, 



4 



Ami 



ml log(/52) - m2(2mi - 2ml + "^2) log(/3o) 



K 



— — {m^ — SttIq + hm^] 



-I - II 



Ami 



m\ log(/?o) - ml{2mi - 2m{ + ml) log(/?i) 



— '-^{ml — Sml + SttIq) 



-'nn — 



'00 



7-02 
-'11 



7-00 

-'11 



Ami 
1 



m^ log(/32) + log(/3i) + — ^ + -(3m? + Smj - ml) 



Ami 

7-0 7-02 
-'1 -'11 



6mo 4 



log(/52) + log(/3o) + ^ + ylSm^ + Sm^ - mj) 



6mi A 



Ami 



2ml{m\ + ml- ml) log(/3o) + ^ + 2Kml 



-'oo 



-"o -'00 



Ami 



2m2(mo + m2 - ) log(Pi) + - — ^ + 2fi;m2 



7-22 tZ tUZ 

-'11 — -'1 -'11 — 



2 7-02 



Ami 



K 



2mJmQ + - 777.2) log(/^2) + 7; — 9 + 2/tmo 

bmf 



Note the symmetries in 1 and 2. 
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